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ABSTRACT 

The  Numerical  Electromagnetics  Code  (NEC)  was  used  to 
evaluate  the  admittance,  average  power  gain,  and  the  electric 
near  and  far  field  of  a  monopole  antenna  mounted  on  a  cubical 
box  over  a  perfectly  conducting  ground  plane.  Two  models  of 
the  box,  employing  surface  patches  and  wire  grids,  were 
evaluated.  The  monopole  was  positioned  at  the  center,  the 
edge,  and  at  a  corner  of  the  box's  top  surface.  Admittance 
and  average  power  gain  of  the  antenna  were  calculated.  NEC 
results  were  examined  and  compared  with  experimental  data  and 
with  results  from  "PATCH",  another  independent 
electromagnetic  modeling  code.  The  near  electric  field  was 
calculated  for  both  models.  Computer  graphics  techniques 
were  presented  for  plotting  NEC  near  field  results  using 
DISSPLA  (Display  Integrated  Software  System  and  Plotting 
LAnguage) ,  a  commercial  graphics  package.  Contour  and  3-D 
amplitude,  and  phase  plots  of  the  near  electric  fields  were 
presented.  Radiation  patterns  were  calculated  to  relate  far 
field  and  near  field  behavior  of  the  antenna.  Surface  patch 
and  wire  grid  models  are  compared  and  conclusions  were 
presented. 
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I.   INTRODUCTION 

A.   GENERAL  BACKGROUND 

The  arrangement  of  HF  antennas  aboard  ships  requires  the 
consideration  of  many  factors.  Among  these  are:  input 
impedance,  gain,  blockage  angles,  patterns,  weight,  firing 
arcs  of  weapons,  and  transmission  line  length  restrictions. 
Additional  factors  that  greatly  influence  naval  operations 
are  RF  radiation  hazards  (RADHAZ)  and  electromagnetic 
interference  (EMI) . 

Shipboard  operations  are  carried  out  within  fixed  (small) 
distances  from  HF  transmitting  antennas,  so  the  Navy  has  a 
unique  and  long-standing  operational  problem--the  radiation 
from  these  antennas  can  be  hazardous  to  personnel,  ordnance, 
fuel,  and  electronic  equipment  because  of  the  intensity  of 
the  fields  close  to  the  radiating  elements  (near  fields) . 
The  goal  of  the  antenna  designer  is  to  select  a  location  for 
the  antenna  that  permits  it  to  perform  with  a  minimum  of 
restrictions  placed  on  personnel  activity,  electronic 
communication  equipment,  fuels  and  ordnance. 

In  the  selection  of  a  site  for  an  antenna,  it  becomes 
evident  that  all  of  the  desired  requirements  cannot  be  met, 
therefore,  trade-offs  have  to  be  made  based  on  a  priority 
system  determined  by  the  ship's  primary  mission.  For 
instance,  on  a  combatant  ship,  placing  antennas  relative  to 
gun  mounts,  missile  launchers,  and  magazines  are   of  major 
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importance.  On  an  auxiliary,  however,  the  underway 
replenishment  gear  might  receive  the  greatest  attention. 

For  the  above  reasons,  the  Navy  has  recognized  the  need 
and  focused  its  interest  on  the  investigation  of  the  near 
fields  of  antennas  for  a  number  of  years. 

Near  field  structure  is  very  complex,  and  previous 
theoretical  analysis  is  practical  for  only  very  simple 
antennas  in  uncomplicated  geometrical  settings.  With  the 
advent  of  high-speed  computers,  approximate  solution 
technigues  such  as  the  Method  of  Moments  became  practical. 
This  technigue  is  suitable  for  the  calculation  of 
electromagnetic  fields  anywhere,  including  in  close  proximity 
to  the  radiating  element.  The  degree  of  accuracy  of  the 
solution  is  a  function  of  the  number  of  computations  reguired 
to  obtain  that  solution  (computer-time) . 

The  Method  of  Moments  technigue  is  the  theoretical  basis 
for  the  Numerical  Electromagnetics  Code  (NEC) ,  which  is  a 
code  for  the  simulation  and  analysis  of  the  electromagnetic 
response  of  antennas  and  other  metal  structures  [Ref.  1]  . 
NEC  is  the  computer  simulation  tool  that  was  used  in  this 
investigation  of  near  fields. 

B.   REASONS  FOR  THIS  STUDY 

1.   Hazards  of  Electromagnetic  Radiation  (RADHAZ) 

Since  near  fields  of  Navy  HF  antennas  are  primarily 
responsible  for  radiation  hazards  (RADHAZ) ,  the  accuracy  of 
near   field  predictions  must  be  adeguate  to  ensure  that 
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personnel  and  equipment  can  continue  to  function  safely  in 
modern  shipboard  environments.  There  are  three  types  of 
radiation  hazards: 

a.  Hazards  of  Electromagnetic  Radiation  to  Personnel 

(HERP) . 

b.  Hazards  of  Electromagnetic  Radiation  to  Ordnance 

(HERO) . 

c.  Hazards  of  Electromagnetic  Radiation  to  Fuels 

(HERF) . 

These  hazards  and  their  general  near  field  criteria  are 
discussed  in  detail  in  Reference  2 . 

2 .   Electromagnetic  Interference  (EMI) 

Near  fields  of  shipboard  antennas  can  also  be  the 
cause  of  radio  frequency  interference.  On  board  ships  RF 
interference  problems  can  be  quite  complex.  In  land 
installations,  there  are  several  preventive  measures  for 
reducing  interference,  such  as  isolation  of  communication 
related  equipment  and  separation  of  transmitting  and 
receiving  antennas  to  avoid  near  field  interaction.  On  a 
ship,  it  is  almost  impossible  to  do  this.  The  design  is 
constrained  by  the  nature  of  the  ship  and  care  has  to  be 
taken  in  operating  radiating  equipment  in  order  to  minimize 
interferences . 

Because  shipboard  transmitting  and  receiving  antennas 
are  close  together  and  the  ship's  superstructure  provides  a 
very  large  nonlinear  surface  to  interact  with  communications 
signals   (near  field  interaction) ,   intermodulation  problems 


are  likely  because  of  currents  induced  by  nonlinear  devices 
in  the  near  field  region  of  an  antenna. 

Intermodulation  sources  are  nonlinear  items  present 
in  the  immediate  surroundings  of  the  receive  or  transmit 
antennas.  Metallic  life  lines  are  major  sources  of 
interference;  expansion  joints  and  mooring  or  anchor  chains 
used  to  secure  the  ship  at  the  pier  are  found  to  introduce 
significant  nonlinearities.  Other  nonlinear  items  include 
ladders,  guard  rails,  antenna  guying  wires,  and  booms.  [Ref. 
3] 

These  environmental  nonlinear  sources  can  be 
classified  as  two  major  categories,  (1)  nonlinear  resistive 
"junctions"  created  by  metallic  mating  surfaces,  and  (2) 
ferromagnetic  materials.  Ferromagnetic  materials,  such  as 
steel,  generate  interfering  frequencies  when  driven  by  strong 
fundamental  currents  caused  by  the  near  fields  of  an  antenna. 

Table  1  from  Reference  4  is  a  listing  of  a  number  of 
significant  losses  in  lives,  fleet  assets,  and  dollars  due  to 
RADHAZ  and  due  to  Electromagnetic  Interference  (EMI) . 


TABLE  1  [Ref.  4] 
MAJOR  LOSSES  DUE  TO  EMI  OR  EMI  "FIXES" 

During  Time  of  Military  Active  Environment: 

—  32  aircraft/134  crew  lost,  $172M  damage,  EMI  triggered  A/C 
rocket  detonation  on  USS  Forrestal  flight  deck  -  Late  60 's. 

—  Most  of  the  topside  destroyed,  crew  injuries,  $100+M 
damage,  anti-radiation  missile  detected  friendly  cruiser  - 
Late  60 's. 

—  All  eguipment  lost,  crew  fatalities,  $10+M  damage,  anti- 
radiation  missile  detected  USMC  radar  -  Early  70' s. 

—  Frigate  lost,  crew  fatalities,  $200+M  replacement,  HMS 
Sheffield  secured  ASMD  radar  to  use  SATCOM  to  avoid  EMI-1984. 

During  Peacetime  Environment  -  Late  70' s  thru  1986: 

—  Planes  Crashed:   NATO  aircraft  (A/C)  overseas  -  A/C  lost, 
crew  fatalities,  $10+M  replacement.   USN  aircraft  in 
California  -  A/C  lost,  crew  fatalities,  $10+M  replacement. 

—  FFG  ASW  helo  while  hovering  over  FFG,  ship  radiation 
caused  blade  servo  fault,  helo  damage,  crew  injuries,  $750K 
damage. 

Ship  Missiles  Lost  and  Near  Miss  Hits  of  USN  Ships: 

—  Cruiser  destroyed  missile  to  avoid  another  cruiser,  1 
missile  and  1  target  drone  lost,  $10+M. 

Ship  Damage  Due  to  EMI: 

—  Terrier  launcher  rail  fire,  target  drone  lost,  $100K+ 
replacement. 

—  Craft  dropped  from  foilborne  operation  due  to  ship  control 
system  EMI  -  Bow  damaged,  crew  injuries  unknown,  $10+M. 

—  Fires  in  propulsion  room  due  to  EMI  "fixes"  installed  to 
prevent  total  shutdown  of  plant  in  presence  of  RF,  minor 
damage,  no  injuries,  crew  repaired. 

Equipment  Damage  Due  to  EM  Environment: 

—  ECM  eguipment  burned  out  by  radiator  on  neighbor  ship, 
equipment  lost,  $1M+  replacement. 

Sailors  RF  Burned: 

Common  occurrence  for  crew.   Increases  noted  from  1982-1985: 

—  Personnel  RADHAZ      -  ten-fold  increase 

—  Personnel  RF  burns    -  six-fold  increase 

—  Hazards  to  ordnance  -  five-fold  increase 


C.   SCOPE  OF  THE  STUDY 

The  Numerical  Electromagnetics  Code   (NEC)   is  used  to 
evaluate  the  effects  of  the  shipboard  topside  environment  on 
the  admittance  and  the  electric-near  field  structure  of  a 
6  cm  monopole  antenna  mounted  on  a  cubical  box. 

A  ship's  superstructure  has  been  approximated  by  two 
theoretical  scaled  models  of  cube-shaped  boxes  of  0.1  m  sides 
(  ^/3  at  a  frequency  of  1  GHZ)  over  a  perfectly  conducting 
ground  plane.  For  the  first  model,  the  surface  patch 
modeling  technique  is  used  and  for  the  second,  the  wire  grid 
modeling  technique. 

The  antenna  is  /y5  at  a  frequency  of  1  GHZ  and  is  mounted 
at  the  center,  edge,  and  corner  of  the  top  surface.  NEC 
results  for  admittance  are  obtained  for  these  three 
positions. 

For  the  surface  patch  model,  the  results  are  compared  to 
measured  experimental  data  from  References  5  and  6.  For  the 
same  model,  NEC  results  are  compared  with  "PATCH"  code 
results.  "PATCH"  is  another  independent  code  recently 
developed  using  a  much  different  theoretical  and  numerical 
algorithm  implementation  [Ref.  7].  The  computational 
accuracy  of  NEC  is  then  examined  and  discussed.  For  the  wire 
grid  model,  the  admittance  characteristic  results  from 
Reference  8  are  compared  with  the  experimental  data  mentioned 
above.  Comparison  between  surface  patch  and  wire  grid 
numerical  modeling  techniques  are  made. 


NEC  is  then  used  for  calculation  of  near  electric  fields 
for  surface  patch  and  wire  grid  models.  Fortran  programs 
were  constructed  for  plotting  NEC  results  of  near  electric 
fields  in  conjunction  with  DISSPLA  (Display  Integrated 
Software  System  and  Plotting  LAnguage) ,  [Ref.  9],  which  is  a 
commercial  proprietary  plotting  package. 

Contour  plots  and  3-D  plots  of  near  electric  field  are 
presented.  The  phase  variation  of  the  near  electric  field  is 
examined  and  plotted  also.  The  effect  of  the  shipboard 
topside  environment  on  the  electric  near  field  of  the 
monopole  antenna  is  discussed. 

For  the  surface  patch  model  radiation  patterns  were 
calculated  in  order  to  relate  the  far  field  and  near  field 
behavior  of  the  antenna. 

D.   THEORETICAL  BACKGROUND 
1.   Near  Field  Theory 

The  fields  around  an  antenna  may  be  divided  into  two 
regions,  one  near  the  antenna  called  the  near  field  or 
Fresnel  zone  and  one  at  a  large  distance  called  the  far- 
field  or  Fraunhofer  zone.  Figure  la  indicates  these  two 
regions. 
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Figure  la.  Fields  of  a  Radiating  Antenna  [Ref.  8] 
The  boundary  between  the  far  field  and  near  field  is 
the  distance  r  =  2D2/^  ,  where  D  is  the  maximum  length  of  the 
antenna  in  meters  and  2  is  the  wavelength  in  meters.  The 
distance  from  the  surface  of  the  antenna  to  this  boundary  is 
called  the  near  field  region,  while  beyond  this  boundary  the 
region  is  called  far  field.  The  near  field  region  is  divided 
into  two  subregions  —  the  reactive  and  radiating  near  field. 
The  reactive  near  field  is  extended  to  2/2ir  from  the 
antenna's  surface,  while  in  practice  a  distance  of  /\ 
represents  this  boundary.  The  phase  of  the  magnetic  and 
electric  field  is  almost  quadrature  in  regions  within  a 
wavelength  of  the  antenna  (reactive  near  field) .  Beyond  the 
distance  of  a  wavelength,  the  electric  and  magnetic  fields 
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are  propagating  in  phase  (radiating  near  field)  until  the  far 
field  is  reached.  In  the  far  field,  the  shape  of  the  field 
pattern  is  independent  of  the  distance,  while  in  the  near 
field  the  shape  depends  on  it. 

Consider  a  dipole  antenna  which  is  enclosed  in  an 
imaginary  sphere.  Figure  lb  shows  the  energy  flow  near  the 
dipole. 
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Figure  lb.   Energy  Flow  Near  a  Dipole  Antenna  [Ref.  8] 


The  waves  are  reflected  near  the  poles  while  in  the  eguator 
are  diffusing  perpendicular  to  the  dipole.  Conseguently , 
energy  is  oscillating  near  the  antenna  while  part  of  it  flows 
outward  in  the  region  of  the  equator.  The  oscillating  energy 
represents  reactive  power,  while  the  outward  flow  represents 
radiated  power.  The  near  field  is  related  to  charge  and 
current  density  while  the  far  field  is  related  only  to 
current. 


2 .   Numerical  Electromagnetics  Code  (NEC)  -  Method  of 
Moments 

NEC  is  a  computer  code  which  analyzes  the 
electromagnetic  reponse  of  antennas.  The  program  uses  an 
electric  field  integral  equation  (EFIE)  and  a  magnetic  field 
integral  equation  (MFIE)  to  do  this  analysis.  The  EFIE  and 
MFIE  integral  equations  are  solved  in  NEC,  for  current 
distribution  using  the  method  of  moments  approach.  This 
approach  reduces  the  integral  equation  to  a  system  of  linear 
algebraic  equations  in  the  form  of  a  matrix,  where  the 
unknowns  are  coefficients  in  an  expansion  of  the  current. 
The  matrix  equations  can  be  solved  for  the  current  in  a  high- 
speed computer.  There  are  two  limitations  to  the  moment 
method:  (1)  the  amount  of  computer  storage  necessary  for  the 
elements  of  the  impedance  matrix  and  (2)  the  time  required  to 
calculate  those  elements  and  solve  the  system  of  equations. 
When  the  current  distribution  is  computed,  all  effects  of  the 
antenna  structure  can  be  determined.  [Ref.  1] 

During  this  investigation,  two  modeling  approaches 
are  used,  surface  patch  modeling  and  wire  grid  modeling.  In 
surface  patch  modeling  a  two-dimensional  field  equation  and 
appropriate  surface  boundary  conditions  are  applied  to  solve 
for  two  orthogonal  components  of  current  on  small  patches 
representing  the  actual  surface.  Either  an  electric  field 
integral  equation  (EFIE)  or  a  magnetic  field  integral 
equation  (MFIE)  may  be  used  theoretically.  NEC  employs  the 
MFIE  method  while  the  code  mentioned  previously,  "PATCH", 
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uses  an  EFIE  method.  The  wire  grid  modeling  approach 
approximates  the  antennas  and  the  surfaces  on  the  ship  with 
thin  wires  using  the  EFIE  method.  Surfaces  are  represented 
as  wire  grids  in  which  the  grid  cells  are  much  smaller  than 
the  wavelength.  Either  modeling  approach  can  provide  data 
upon  which  engineering  judgments  can  be  made.  The 
investigations  of  the  two  modeling  approaches  in  the  past 
shows  that  wire  gridding  provides  slightly  more  accurate  data 
in  comparison  to  surface  patch  modeling.  This  is  another 
point  of  investigation  of  this  thesis.  However,  surface 
patch  modeling  can  apparently  offer  some  amount  of  savings  in 
computer  time  cost.  [Ref.  2] 
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II.   SURFACE  PATCH  MODELING 

A.   GUIDELINES  FOR  MODELING  USING  NEC 

A  conducting  surface  is  modeled  by  means  of  multiple, 
small  flat  surface  patches,  corresponding  to  segments  used  to 
model  wires.  The  patches  are  chosen  to  completely  cover  the 
surface  to  be  modeled.  The  parameters  defining  the  surface 
patch  are  the  cartesian  coordinates  of  the  patch  center,  the 
components  of  the  outward-directed,  unit  normal  vector  and 
the  patch  area.  The  program  computes  the  surface  current  on 
each  patch  along  the  orthogonal  unit  vectors  which  are 
tangent  to  the  surface.  Four  patch  shape  options  exist: 
arbitrary  shape,  rectangular,  triangular,  and  quadrilateral. 
For  the  rectangular,  triangular,  and  quadrilateral  patches, 
the  outward  normal  vector  fi  is  specified  by  ordering  of 
corners  1,  2,  and  3  using  the  right-hand  rule.   The  vectors 

A  * 

t-L  are  parallel  to  the  side  from  corner  1  to  corner  2,  and  t2 
=  n  x  t]_.  The  coordinates  of  corners  1,  2,  and  3  have  to  be 
specified  in  the  program.  Figure  2  shows  how  the  multiple 
patches  are  defined  for  a  rectangular  surface.  Rectangular 
patches  were  used  for  the  modeling  purposes  of  this  thesis. 
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Figure  2 .   Rectangular  Surface  Covered  by  Multiple  Patches 
[Ref.  1] 


When  a  wire  is  connected  to  a  surface,  the  wire  must  end 
at  the  center  of  a  patch  with  identical  coordinates  used  for 
the  wire  end  and  the  patch  center.  The  program  then  divides 
the  patch  into  four  equal  patches  about  the  wire  end.  Only  a 
single  wire  may  connect  to  a  given  patch  and  a  segment  may 
have  a  patch  connection  on  only  one  of  its  ends.  A  wire  may 
never  connect  to  a  patch  formed  by  subdividing  another  patch 
for  a  previous  connection.  The  use  of  surface  patches  is 
restricted  to  modeling  voluminous  bodies.  Very  long  narrow 
patches  should  be  avoided  when  subdividing  the  surface.  [Ref. 
1] 


B.   LIMITATIONS  AND  OPTIMIZATION  CRITERIA 

The  accuracy  of  the  results  depends  on  the  patch  size 
measured  in  wavelengths.  A  minimum  of  about  25  patches 
should  be  used  per  square  wavelength  of  surface  area,  with 
the  maximum  size  for  an  individual  patch  less  than  or  equal 
to  0.04  square  wavelengths. 
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In  NPGNEC  (an  NPS  customized  version  of  NEC-3)  which  was 
used  for  simulation  purposes  in  this  thesis,  the  number  of 
patches  plus  the  number  of  segments  must  be  less  than  3  00. 
When  the  total  number  is  greater  than  3  00  and  less  than 
1,000,  the  NPG1000  version  of  NEC-3  must  be  used;  this  is 
done  for  the  wire  grid  modeling  in  the  second  part  of  this 
thesis.  Since  the  program  does  not  integrate  over  patches 
except  at  a  wire  connection  (active  patch) ,  the  patch  shape 
does  not  affect  the  results.  The  parameters  affecting  the 
results  are  the  location  of  the  patch  centroid,  the  patch 
area,  and  the  outward  unit  normal  vector. 

The  optimum  model  for  complex  structures  which  the  first 
part  of  the  thesis  investigated  can  be  estimated  by  varying 
the  segment  and  patch  density,  observing  the  results  and  the 
convergence  of  the  solution.  In  the  case  of  an  edge- 
mounted  antenna,  the  accuracy  of  the  results  are  expected  to 
depend  upon  the  size  of  the  patches  near  the  edge.  Smaller 
patches  are  suggested  at  edge  areas  since  the  current 
magnitude  may  vary  rapidly  in  this  region. 

C.   PREVIOUS  STUDIES 

In  References  5  and  6,  experimental  and  computational 
investigations  were  performed  to  determine  the  admittance 
characteristics  of  a  monopole  antenna  mounted  on  a  cubical 
conducting  box  of  0.1  m  sides  (3/3  at  frequency  1  GHz)  over  a 
ground  plane.  The  antenna,  a  6  cm  monopole  (  ^/5  at  frequency 
1  GHz) ,  was  tested  for  three  different  mounting  positions  on 
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the  top  surface.  Experimental  admittance  data  for  the  6  cm 
monopole  antenna  were  plotted  vs.  frequency  [Ref.  5]. 
Numerically  calculated  data  using  the  "PATCH"  computer  code 
were  also  tabulated  [Ref.  5],  "PATCH"  is  a  frequency  domain 
electromagnetic  analysis  code  based  on  a  method  of  moments 
solution  to  the  Electric  Field  Integral  Equation  (EFIE)  [Ref. 
7].  In  this  code  objects  are  modeled  by  planar  triangular 
patches  which  easily  conform  to  surfaces  and  boundaries  of 
general  shape  and  allow  variable  patch  densities  over  the 
surface  of  the  object. 

D.   MODEL  OF  6  CM  MONOPOLE  ANTENNA  MOUNTED  ON  A  CUBICAL 

CONDUCTING  BOX — ADMITTANCE  AND  AVERAGE  GAIN  INVESTIGATION 

Monopole  antennas  mounted  on  arbitrary  shaped  objects 
have  received  much  attention  because  they  represent  practical 
situations  where  antennas  are  mounted  on  aircraft,  ships, 
vehicles,  and  buildings.  The  antennas  are  not  always 
centrally  mounted  and  are  often  located  near  edges  and 
corners. 

Super-structures  of  Navy  ships  can  be  modeled  as 
conducting  boxes  which  form  a  "ship-like"  structure.  Navy 
shipboard  topside  antennas  are  often  monopoles  and  are 
attached  to  the  top  of  these  boxes,  in  different  locations. 

The  numerical  model  of  a  cubical  five-sided  box  of 
0.1  m  per  side  was  constructed  using  NEC  (the  bottom  was  not 
included  as  a  surface  since  the  box  was  placed  on  a  perfectly 
conducting  ground  plane) .   A  6  cm  monopole  antenna  was  placed 
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at  the  center,  at  the  edge  (3.63  cm  from  center)  and  at  a 
corner  (5.14  cm  on  a  diagonal  from  center),  as  shown  in 
Figure  3 .  This  model  configuration  was  selected  since 
experimental  data  was  available  for  it  for  comparison  from 
Reference  5. 


Figure  3 .   Geometry  of  Monopole  Attached  to  the 
Conducting  Box  Over  a  Ground  Plane 


The  first  part  of  the  investigation  checked  the  input 
impedance  and  average  gain  of  the  model  using  NEC-3  (NPGNEC) . 
Comparisons  of  calculated  values  for  input  impedance  and 
average  gain  were  made  to  experimental  values  from  Reference 
5  in  order  to  validate  the  computer  model  before  proceeding 
with  an  investigation  of  near  fields. 

This  modeling  exercise  compared  driving  point  admittance 
from  NEC  as  patch  density  is  varied,  to  both  experimental 
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values  and  PATCH  code  results  from  Reference  5.  The  number 
of  patches  on  the  top  of  the  box  was  varied  in  search  of  an 
optimum  value  of  surface  samples,  which  would  later  be  used 
for  near  field  calculations.  The  top  was  subdivided  to 
retain  symmetry,  as  much  as  possible,  and  to  closely  match 
positions  of  the  antenna  on  the  experimental  model.  The 
monopole  was  positioned  at  the  center,  at  an  edge  (3.63  cm 
from  the  center)  and  at  a  corner  (5.14  cm  on  a  diagonal  from 
the   center),  very  close  to  the  experimental  model  positions. 

Nine  frequencies  from  1  to  1.4  GHZ  were  selected.  NEC 
input  data  files  for  the  three  different  locations  of  the 
monopole  are  presented  in  Appendix  A. 

The  Numerical  Green's  Function  option  (NGF)  was  used  to 
exploit  partial  symmetry,  reducing  solution  time.  The 
nonsymmetric  part,  the  6  cm  monopole  antenna,  was  added  to 
the  NGF  file  of  the  box  for  impedance  and  average  gain 
calculations. 

1 .   Monopole  at  the  Center 

The  monopole  was  divided  into  five  segments  and 
placed  at  center  of  the  top  surface.  The  top  surface  was 
divided  into  25  (5x5),  49  (7x7),  81  (9x9),  and  121  (11x11) 
patches.  The  impedance  (in  OHMS)  and  admittance  (in  mMHOS) 
are  listed  in  Tables  4,  5,  6,  and  7  for  the  nine  frequencies. 
Experimental  data  from  Reference  5  are  shown  in  Table  2 .  NEC 
statistics  for  the  four  top  surface  configurations  are  in 
Table  3.   The  number  of  patches  used  for  each  case,  the  patch 
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area,  and  the  variation  of  run  time  are  shown.  Double 
precision  NEC  (DNPGNEC)  was  used  for  the  5x5  and  7x7  patch 
case  to  investigate  the  possibility  of  better  accuracy  than 
single  precision  NEC.   The  results  are  in  Tables  8  and  9. 

To  gain  insight  into  the  accuracy  of  the  results, 
the  average  power  gain  was  calculated.  Correct  average  power 
gain  does  not  necessarily  insure  accurate  results. 
Therefore,  input  impedance  was  also  checked.  The  average 
power  gain  provides  a  good  check  of  the  solution  accuracy  if 
it  is  close  to  two  for  radiating  systems  located  over 
perfectly  conducting  ground  which  is  our  case,  or  close  to 
one  for  free  space  conditions.  Average  power  gain  within 
10%  of  theoretical  is  adequate  for  engineering  purposes.  For 
all  cases,  the  calculated  average  power  gain  is  between  1.96 
(minimum)  to  2.02  (maximum)  which  is  in  good  agreement  with 
theory. 

NEC  results  and  measurements,  for  the  monopole  at  the 
center,  while  the  patch  density  varies  on  the  top  surface 
(Case  1:  5x5,  Case  2:  7x7,  Case  3:  9x9,  and  Case  4:  11x11), 
are  presented  in  Figures  4,  6,  8,  and  10.  In  Figures  5,  7, 
9,  and  11,  the  PATCH  code  results  have  been  superimposed  on 
those  of  NEC  and  the  measurements  for  comparison.  The  PATCH 
code  is  compared  to  experimental  results  in  Figure  12. 

Single  precision  (NPGNEC)  calculations  vs  double 
precision  (DNPGNEC)  results  (Tables  4,  6,  8,  and  9)  for 
conductance  and  susceptance  vary  only  in  the  second  and 
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third  decimal  places.    This  difference  is  negligible  and 
cannot  be  observed  in  a  plot.   For  all  cases,  the  antenna  was 
fed  at  the  base  segment. 
CASE  1:   5x5  patches  on  top  surface  (Figures  4  and  5) 

The  correlation  between  measurements  and  NEC 
calculations  for  admittance  in  Figure  4,  was  generally  quite 
poor.  Apparently,  the  solution  has  not  converged.  There  is 
good  correlation  for  susceptance  between  1.2  and  1.4  GHz. 
PATCH  results  of  Figure  5  correlate  well  with  measurements  in 
comparison  to  NEC. 
CASE  2:   7X7  patches  on  top  surface  (Figures  6  and  7) 

NEC  values  of  admittance  correlate  well  with 
measurements  (Figure  6) ,  but  the  run  time  observed  in  Table  3 
is  almost  double  compared  to  Case  1.  The  solution  has 
apparently  converged.  The  increased  accuracy  is  probably  due 
to  the  top  of  the  surface  patch  box  having  been  more  finely 
divided  into  smaller  patches  which  provides  finer  details  of 
surface  currents.  Observing  Figure  7,  we  see  that  NEC  has 
somewhat  better  correlation  with  measurements  than  PATCH 
does. 
CASE  3:   9x9  patches  on  top  surface  (Figures  8  and  9) 

There  is  good  correlation  of  NEC  results  with 
measurements  in  Figure  8,  especially  in  conductance,  but  not 
quite  as  good  as  in  Case  2 . 

NEC  performs  better  than  PATCH  for  both  conductance 
and  susceptance.    As  the  number  of  patches  increases  (81 
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patches  on  the  top  surface  vs  49  for  Case  2  and  25  for  Case 
1) ,  the  run  time  is  doubled  and  tripled  correspondingly. 
CASE  4:  11x11  patches  on  the  top  surface  (Figures  10  and  11) 
The  correlation  of  NEC  results  is  good  especially  in 
conductance  (Figure  10)  and  is  better  than  PATCH  (Figure  11) . 
Run  time  is  doubled  compared  to  Case  3 ,  as  can  be  observed 
from  Table  3 . 

GENERAL  RESULTS:  The  solution  seems  to  converge  at  Case  2 
with  7x7  =  49  patches  on  the  top  surface,  patch  area  0.0002  0 
m2  (0.0022  square  wavelengths  at  1  GHz)  and  0.00005  m2 
(0.00055  square  wavelengths  at  1  GHz)  for  each  of  the  four 
central  "subpatches"  where  NEC  subdivides  a  normal  size  patch 
at  the  connection  point  to  the  monopole.  NEC  gives  excellent 
results  compared  to  measurements  for  Cases  2,  3,  and  4  and 
has  better  performance  than  PATCH  for  these  three  cases. 
PATCH  results  are  better  than  NEC  results  for  Case  1  (5x5  = 
25  patches  on  the  top) .  For  easier  comparison  between  the 
two  codes,  NEC  vs.  PATCH  admittance  values  for  the  Cases  1, 
2,  3  and  4  are  plotted  separately  in  Figures  13,  14,  15,  and 
16.  In  Figures  17  and  18,  NEC's  conductance  and  susceptance 
vs  measurements  as  the  top  surface  patch  density  increases 
(5x5,  7x7,  9x9,  and  11x11  patches)  are  available  for 
comparison.  Consequently,  the  optimum  model  appears  to  be 
Case  2  (7x7  =  49  patches  on  the  top)  ,  the  one  that  is 
selected  for  the  near  field  investigation  for  surface  patch 
numerical  model. 


20 


TABLE  2 
MEASURED  DATA  FOR  A  6CM  MONOPOLE  ON  A  CUBICAL  BOX 


FREQUENCY 
(GHZ) 

6CM  MONOPOLE  POSITION 

ADMITTANCE  (mMHOS) 

AT  CENTER 

AT  EDGE  3.5cm 
FROM  CENTER 

AT  CORNER  5.15cm 
FROM  CENTER 

1.00 

5.0 

+  J20.0 

7.0   +  J17.0 

8.0   +  J12.0 

1.05 

14.0 

+  J28.0 

13.0   +  J19.5 

12.0   +  J13.0 

1.1 

37.0 

+  J24.0 

24.0   +  J17.0 

18.0   +  jlO.O 

1.15 

41.0 

-  j  8.0 

32.0   +  j  3.0 

21.0   +  j  4.0 

1.2 

26.0 

-  J16.0 

26.0   -  j  7.5 

20.0   -  j  2.0 

1.225 

19.0 

-  J15.0 

22.0   -  jlO.O 

17.0   -  j  4.5 

1.3 

11.0 

-  jio.o 

12.0   -  j  9.0 

12.0   -  j  5.0 

1.4 

7.0 

-  j  6.0 

7.0   -  j  6.0 

8.0   -  j  3.0 
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TABLE  3 

NEC  STATISTICS  FOR  MONOPOLE  AT  CENTER  OF  PATCH  BOX 

AS  THE  NUMBER  OF  PATCHES  ON  TOP  INCREASES 

(NPS  MAINFRAME  COMPUTER  -  IBM  SYSTEM  370) 


CASE 

NUMBER  OF 
PATCHES 

PATCH 
AREA  (m2) 

RUN  TIME   (sec) 
SINGLE  PRECISION  NEC 

1 

surface  on  sides 
18  x  4  =         72 

0.00056 

13.937 

surface  on  top 
5x5=         25 

0.00040 

monopole  attach- 
ment patch       4 

0.00010 

6.747 

TOTAL  # PATCHES  101 

2 

surface  on  sides 
18  x  4  =         72 

0.00056 

25.883 

surface  on  top 
7x7=        49 

0.00020 

monopole  attach- 
ment patch       4 

0.00005 

8.940 

TOTAL  # PATCHES  12  5 

3 

surface  on  sides 
18  x  4  =         72 

0.00056 

54.170 

surface  on  top 
9x9=         81 

0.00012 

monopole  attach- 
ment patch       4 

0.00003 

13.107 

TOTAL  # PATCHES  157 

4 

surface  on  sides 
18  x  4  =         72 

0.00056 

107.00 

surface  on  top 

11  x  11  =      121 

0.00008 

monopole  attach- 
ment patch       4 

0.00002 

17.767 

TOTAL  # PATCHES  197 
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TABLE  4 
CALCULATED  DATA 
MONOPOLE  AT  CENTER  OF  SURFACE  PATCH  BOX 
(5X5  PATCHES  ON  THE  TOP  SURFACE  -  BASE  SEGMENT  FEED) 


FREQUENCY 
(GHZ) 

AVERAGE 
POWER 
GAIN 

IMPEDANCE  (OHMS) 
ADMITTANCE  (mMHOS) 

1.00 

2.02 

10.82   -  J20.21 
20.58   +  J38.45 

1.05 

2.02 

14.38   -  j  5.51 
60.64   +  J23.24 

1.10 

2.01 

18.86   +  j  8.57 
43.94   -  J19.97 

1.15 

2.01 

24.29   +  J21.89 
22.72   -  J20.47 

1.175 

2.01 

27.34   +  J28.20 
17.72   -  J18.28 

1.2 

2.01 

30.60   +  J34.25 
14.50   -  J16.24 

1.225 

2.01 

34.03   +  J39.99 
12.34   -  J14.50 

1.3 

2.01 

44.86   +  J55.32 
8.84   -  J10.90 

1.4 

2.02 

58.47   +  J71.83 
6.81   -  j  8.37 
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TABLE  5 
CALCULATED  DATA 
MONOPOLE  AT  CENTER  OF  SURFACE  PATCH  BOX 
(7X7  PATCHES  ON  THE  TOP  SURFACE  -  BASE  SEGMENT  FEED) 


FREQUENCY 
(GHZ) 

AVERAGE 
POWER 
GAIN 

IMPEDANCE  (OHMS) 
ADMITTANCE  (mMHOS) 

1.00 

1.99 

10.99   -  J37.53 
7.19   +  J24.54 

1.05 

1.99 

14.60   -  J21.82 
21.18   +  J31.65 

1.10 

1.99 

19.14   -  j  6.81 
46.36   +  J16.50 

1.15 

1.99 

24.64   +  j  7.36 
37.26   -  jll.12 

1.175 

1.99 

27.74   +  J14.08 
28.66   -  J14.55 

1.2 

1.99 

31.04   +  J20.50 
22.43   -  J14.81 

1.225 

1.99 

34.52   +  J26.62 
18.17   -  J14.01 

1.3 

1.99 

45.53   +  J42.98 
11.61   -  J10.97 

1.4 

1.99 

59.41   +  J60.67 
8.24   -  j  8.41 
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TABLE  6 
CALCULATED  DATA 
MONOPOLE  AT  CENTER  OF  SURFACE  PATCH  BOX 
(9X9  PATCHES  ON  THE  TOP  SURFACE  -  BASE  SEGMENT  FEED) 


FREQUENCY 
(GHZ) 

AVERAGE 
POWER 
GAIN 

IMPEDANCE  (OHMS) 
ADMITTANCE  (mMHOS) 

1.00 

1.98 

11.11   -  J48.08 
4.56   +  J19.74 

1.05 

1.98 

14.75   -  J31.73 
12.05   +  J25.92 

1.10 

1.97 

19.33   -  J16.14 
30.48   +  J25.45 

1.15 

1.97 

24.88   -  j  1.42 
40.06   +  j  2.29 

1.175 

1.97 

28.00   +  j  5.54 
34.36   -  j  6.80 

1.2 

1.97 

31.34   +  J12.22 
27.70   -  J10.80 

1.225 

1.97 

34.84   +  J18.57 
22.35   -  jll.91 

1.3 

1.98 

45.96   +  J35.57 
13.60   -  J10.53 

1.4 

1.98 

59.99   +  J54.01 
9.20   -  j  8.29 
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TABLE  7 
CALCULATED  DATA 
MONOPOLE  AT  CENTER  OF  SURFACE  PATCH  BOX 
(11  X  11  PATCHES  ON  THE  TOP  SURFACE  -  BASE  SEGMENT  FEED) 


FREQUENCY 
(GHZ) 

AVERAGE 
POWER 
GAIN 

IMPEDANCE  (OHMS) 
ADMITTANCE  (mMHOS) 

1.00 

1.97 

11.19   -  J50.03 
4.26   +  J19.04 

1.05 

1.96 

14.85   -  J33.49 
11.06   +  J24.95 

1.10 

1.96 

19.46   -  J17.73 
28.08   +  J25.58 

1.15 

1.96 

25.05   -  j  2.86 
39.41   +  j  4.49 

1.175 

1.96 

28.20   +  j  4.19 
34.69   -  j  5.16 

1.2 

1.96 

31.55   +  J10.94 
28.29   -  j  9.80 

1.225 

1.96 

35.04   +  J17.36 
22.90   -  jll.33 

1.3 

1.97 

46.27   +  J34.56 
13.87   -  J10.36 

1.4 

1.97 

60.42   +  J53.21 
9.32   -  j  8.20 
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TABLE  8 

CALCULATED  DATA 

MONOPOLE  AT  CENTER  OF  SURFACE  PATCH  BOX 

(5X5  PATCHES  ON  THE  TOP  SURFACE  -  BASE  SEGMENT  FEED) 

(DOUBLE  PRECISION  CALCULATION) 


FREQUENCY 
(GHZ) 

AVERAGE 
POWER 
GAIN 

IMPEDANCE  (OHMS) 
ADMITTANCE  (mMHOS) 

1.00 

2.02 

10.82   -  J20.20 
20.60   +  J38.46 

1.05 

2.01 

14.38   -  j  5.49 
60.67   +  J23.19 

1.10 

2.01 

18.86   +  j  8.58 
43.91   -  J19.99 

1.15 

2.01 

24.29   +  J21.90 
22.70   -  J20.47 

1.175 

2.01 

27.35   +  J28.22 
17.71   -  J18.27 

1.2 

2.01 

30.61   +  J34.26 
14.50   -  J16.23 

1.225 

2.01 

34.03   +  J40.00 
12.34   -  J14.50 

1.3 

2.01 

44.87   +  J55.33 
8.84   -  J10.90 

1.4 

2.02 

58.47   +  J71.84 
6.81   -  j  8.37 
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TABLE  9 
CALCULATED  DATA 
MONOPOLE  AT  CENTER  OF  SURFACE  PATCH  BOX 
(7X7  PATCHES  ON  THE  TOP  SURFACE  -  BASE  SEGMENT  FEED) 
(DOUBLE  PRECISION  CALCULATION) 


FREQUENCY 
(GHZ) 

AVERAGE 
POWER 
GAIN 

IMPEDANCE  (OHMS) 
ADMITTANCE  (mMHOS) 

1.00 

1.99 

10.99   -  J37.52 
7.19   +  J24.54 

1.05 

1.99 

14.60   -  J21.81 
21.20   +  J31.66 

1.10 

1.99 

19.14   -  j  6.79 
46.39   +  J16.47 

1.15 

1.99 

24.65   +  j  7.37 
37.23   -  jll.14 

1.175 

1.99 

27.75   +  J14.09 
28.65   -  J14.55 

1.2 

1.99 

31.05   +  J20.51 
22.42   -  J14.81 

1.225 

1.99 

34.52   +  J26.63 
18.16   -  J14.00 

1.3 

1.99 

45.53   +  J42.99 
11.61   -  J10.96 

1.4 

1.99 

59.42   +  J60.67 
8.24   -  j  8.41 
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Figure  4.   Monopole  at  Center  of  Patch  Box 

(5x5  Patches  on  Top) 

NEC  Admittance  vs  Measurements 
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Figure  5.   Monopole  at  Center  of  Patch  Box 

(5x5  Patches  on  Top) 
NEC  Admittance,  vs  Measurements  and  "PATCH"  Code 
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Figure  6.   Monopole  at  Center  of  Patch  Box 
(7x7  Patches  on  Top) 
NEC  Admittance  vs  Measurements 
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Figure    7.       Monopole    at    Center   of   Patch   Box 
(7x7   Patches   on   Top) 
NEC   Admittance    vs   Measurements    and    "PATCH"    Code 
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Figure  8.   Monopole  at  Center  of  Patch  Box 

(9x9  Patches  on  Top) 

NEC  Admittance  vs  Measurements 
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Figure  9.   Monopole  at  Center  of  Patch  Box 
(9x9  Patches  on  Top) 
NEC  Admittance  vs  Measurements  and  "PATCH"  Code 
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Figure  10.   Monopole  at  Center  of  Patch  Box 

(11x11  Patches  on  Top) 

NEC  Admittance  vs  Measurements 
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Figure  11.   Monopole  at  Center  of  Patch  Box 

(11x11  Patches  on  Top) 

NEC  Admittance  vs  Measurements  and  "PATCH"  Code 
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Figure  12.   Monopole  at  Center  of  Patch  Box 
"PATCH"  Code  Admittance  vs  Measurements 
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Figure  13.   Monopole  at  Center  of  Patch  Box 

(5x5  Patches  on  Top) 

NEC  Admittance  vs  "PATCH"  Code 
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Figure  14.   Monopole  at  Center  of  Patch  Box 

(7x7  Patches  on  Top) 

NEC  Admittance  vs  "PATCH"  Code 
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Figure  15.   Monopole  at  Center  of  Patch  Box 

(9x9  Patches  on  Top) 

NEC  Admittance  vs  "PATCH"  Code 
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Figure  16.   Monopole  at  Center  of  Patch  Box 

(11x11  Patches  on  Top) 

NEC  Admittance  vs  "PATCH"  Code 
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Figure  17.   Monopole  at  Center  of  Patch  Box 
NEC  Conductance  vs  Patch  Density 


42 


e 

-2     >■ 

4-i-isi 

Do 

z 

Ed 

X 

o 

=i==C| 

I 

C3              ^ 

"  1 

PQ 

Wiftlr^!C5;x 

CO 

g 

, 

n        > 

T 

u 

oboiul< 

1 

U 

1 

o 

Ed.Sd  Ed'&£,»d 

\ 

2. 

ZjZZZ^ 

\ 

«;==,=;  = 

n  p!<i+  x 

k 

5  ^ 

1 

H 

| 

-  z 

£ 

>h 

Ed 

K 

J3 

U 

yv 

7\ 

Ed 

\ 

W 

yC- 

hJ 

/' ' 

) 

c 

/' 

'^s^ 

c 

, 

• 

- 

.xV 

z 

/i  .-■' 

^,- 

c 

1  <t 

■ 

2 

\'\ 

(O 

e 

A 

e 
© 

Oi          09           OS           Of-           OC           02           OT            0          OT-        02-      OC 

- 

(SOHWITIIW)  33NVld3DSnS 

Figure  18.   Monopole  at  Center  of  Patch  Box 
NEC  Susceptance  vs  Patch  Density 
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2 .   Monopole  at  Edge 

The  6  cm  monopole  was  divided  into  five  segments  and 
placed  at  an  edge  (Fig.  3)  .  The  measured  data  [Ref .  5]  are 
for  a  monopole  attached  to  an  edge  at  a  distance  3.5  cm  from 
the  center.  Two  computer  model  configurations  of  the  top 
surface  were  selected  in  an  attempt  to  match  the  geometry  of 
the  experimental  model:  81  (9x9)  patches  which  corresponds  to 
a  monopole  position  of  3.33  cm  from  the  center  and  121 
(11x11)  patches  which  corresponds  to  3.63  cm.  The  difference 
in  distance  for  the  position  of  the  monopole  in  the  NEC  model 
compared  to  the  actual  physical  geometry,  is  0.17  cm  less  and 
0.13  cm  more,  correspondingly. 

For  each  configuration,  the  freguency  was  varied  from 
1  to  1.4  GHz  in  nine  computer  runs.  The  impedance  (in  OHMS), 
admittance  (in  mMHOS)  and  calculated  values  of  the  average 
gain  are  shown  in  Tables  10  and  11.  DNPGNEC  was  used  for 
the  case  of  9x9  patches  on  top,  in  order  to  investigate  the 
possibility  of  higher  accuracy  than  single  precision  NEC;  the 
results  are  shown  in  Table  12.  From  Tables  10,  11,  and  12, 
the  calculated  average  gain  varies  from  1.89  (minimum)  to 
1.92  (maximum),  which  is  in  good  agreement  with  the 
theoretical  value  of  2.0. 

Results  of  NEC  predictions  vs.  measurements  for  the 
monopole  at  the  edge,  with  top  surface  patch  density  varying 
(Case  1:  9x9,  Case  2:  11x11),  are  presented  in  Figures  19  and 
21.   Figures  20  and  22,  show  PATCH  code  results  of  Reference 
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5  superimposed  with  NEC  calculations  and  include  measurement 
plots  for  comparison.  PATCH  code  and  experimental  results 
are  presented  in  Figure  23.  Experimental  data  for  the 
monopole  located  at  the  edge  are  shown  in  Table  2 .  The  total 
number  of  patches  and  the  patch  area  are  listed  in  Table  3. 
To  compare  performance  between  NEC  and  PATCH  for  each  case, 
see  Figures  24  and  25. 

Comparing  NPGNEC  vs.  DNPGNEC  results  (Tables  11  and 
12)  ,  we  observe  that  the  values  of  conductance  and 
susceptance  differ  only  in  the  second  and  third  decimal 
place.  This  difference  is  negligible  and  cannot  be  observed 
in  plots  of  NPGNEC  vs  DNPGNEC  results.  For  all  cases,  the 
antenna  was  fed  at  the  base  segment. 
CASE  1:  9x9  patches  on  the  top  surface  (Figures  19  and  20) 

The  agreement  is  very  good  for  conductance  (G)  and 
susceptance  (B)  ,  when  NEC  results  are  correlated  with 
measurements  (Figure  19) . 

In  Figure  20,  a  conductance  comparison  between  the 
performance   of   NEC   and   PATCH   vs.   measurements   shows 
performance  to  be  almost  the  same  between  freguencies  of  1  to 
1.15  GHz. 
CASE  2:  11x11  patches  on  the  top  surface  (Figures  21  and  22) 

For  this  configuration,  the  position  of  the  monopole 
in  the  NEC  model  is  closer  to  the  experimental  position  (3.63 
cm  vs  3.5  cm)  than  in  the  9x9  patch  case.  NEC  conductance 
appears  to  follow  the  same  general  trend  as  measurements. 
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PATCH'S   conductance   values   are   17%   closer   to 
measurements  than  NEC.   NEC  and  PATCH  have  almost  identical 
performance  for  susceptance  and  both  have  the  same  resonant 
frequency  of  1.20  GHz  vs  1.16  GHz  for  themeasurement . 
GENERAL  RESULTS: 

In  both  cases,  NEC  results  are  reasonably  close  to 
experimental  data.  For  clarity,  results  for  NEC  and  PATCH 
for  Case  1  are  plotted  in  Figure  24  and  those  for  Case  2  are 
plotted  in  Figure  25. 

NEC  performs  better  than  PATCH  in  Case  1  and  vice 
versa  for  Case  2 .  The  solution  seems  to  converge  at  Case  1 
(9x9  patches  on  top)  .  However,  the  model  of  Case  2  is 
selected  for  the  near  field  calculations  because  the  position 
of  the  monopole  is  closer  to  the  experimental  model  position. 
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TABLE  10 
CALCULATED  DATA 
MONOPOLE  AT  EDGE  OF  SURFACE  PATCH  BOX  (3. 3 3 CM  FROM  CENTER) 
(9X9  PATCHES  ON  THE  TOP  SURFACE  -  BASE  SEGMENT  FEED) 


FREQUENCY 
(GHZ) 

AVERAGE 
POWER 
GAIN 

IMPEDANCE  (OHMS) 
ADMITTANCE  (mMHOS) 

1.00 

1.90 

21.62   -  J52.92 
6.62   +  J16.19 

1.05 

1.90 

24.44   -  J38.36 
11.81   +  J18.54 

1.10 

1.91 

27.79   -  J24.12 
20.52   +  J17.81 

1.15 

1.91 

31.89   -  J10.25 
28.42   +  j  9.13 

1.175 

1.91 

34.26   -  j  3.49 
28.88   +  j  2.94 

1.2 

1.91 

36.85   +  j  3.08 
26.95   -  j  2.26 

1.225 

1.91 

39.64   +  j  9.49 
23.86   -  j  5.71 

1.3 

1.92 

49.09   +  J27.47 
15.51   -  j  8.68 

1.4 

1.92 

63.24   +  J48.51 
9.95    -  j  7.63 
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TABLE  11 
CALCULATED  DATA 
MONOPOLE  AT  EDGE  OF  SURFACE  PATCH  BOX  (3.63CM  FROM  CENTER) 
(11  X  11  PATCHES  ON  THE  TOP  SURFACE  -  BASE  SEGMENT  FEED) 


FREQUENCY 
(GHZ) 

AVERAGE 
POWER 
GAIN 

IMPEDANCE  (OHMS) 
ADMITTANCE  (mMHOS) 

1.00 

1.89 

23.67   -  J57.39 
6.14   +  J14.89 

1.05 

1.89 

26.31   -  J42.92 
10.38   +  J16.93 

1.10 

1.89 

29.41  -  J28.69 

17.42  +  J16.99 

1.15 

1.90 

33.24   -  J14.74 
25.14   +  jll.15 

1.  175 

1.90 

35.47   -  j  7.92 
26.85   +  j  5.99 

1.2 

1.90 

37.92   -  j  1.24 
26.34   +  j  0.86 

1.225 

1.90 

40.59   +  j  5.27 
24.22   -  j  3.15 

1.3 

1.90 

49.75   +  J23.72 
16.38   -  j  7.80 

1.4 

1.91 

63.89   +  J45.58 
10.37   -  j  7.39 
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TABLE  12 

CALCULATED  DATA 

MONOPOLE  AT  EDGE  OF  SURFACE  PATCH  BOX  (3.63CM  FROM  CENTER) 

(11  X  11  PATCHES  ON  THE  TOP  SURFACE  -  BASE  SEGMENT  FEED) 

(DOUBLE  PRECISION  CALCULATION) 


FREQUENCY 
(GHZ) 

AVERAGE 
POWER 
GAIN 

IMPEDANCE  (OHMS) 
ADMITTANCE  (mMHOS) 

1.00 

1.89 

23.68   -  J57.38 
6.14   +  J14.89 

1.05 

1.89 

26.31   -  J42.91 
10.39   +  J16.93 

1.10 

1.89 

29.42  -  J28.67 

17.43  +  J16.99 

1.15 

1.90 

33.24   -  J14.73 
25.14   +  jll.14 

1.175 

1.90 

35.48   -  j  7.90 
26.85   +  j  5.98 

1.2 

1.90 

37.93   -  j  1.23 
26.34   +  j  0.85 

1.225 

1.90 

40.60   +  j  5.29 
24.22   -  j  3.15 

1.3 

1.90 

49.76   +  J23.73 
16.37   -  j  7.80 

1.4 

1.91 

63.90   +  J45.58 
10.37   -  j  7.39 
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Figure  19.   Monopole  at  Edge  of  Patch  Box 

(9x9  Patches  on  Top) 

NEC  Admittance  vs  Measurements 
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Figure  20.   .Monopole  at  Edge  of  Patch  Box 

(9x9  Patches  on  Top) 

NEC  Admittance  vs  Measurements  and  "PATCH"  Code 
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Figure  21.   Monopole  at  Edge  of  Patch  Box 

(11x11  Patches  on  Top) 

NEC  Admittance  vs  Measurements 
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Figure  22.   Monopole  at  Edge  of  Patch  Box 

(11x11  Patches  on  Top) 

NEC  Admittance  vs  Measurements  and  "PATCH"  Code 
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Figure  23.   Monopole  at  Edge  of  Patch  Box 
"PATCH"  Code  Admittance  vs  Measurements 
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Figure  24.   Monopole  at  Edge  of  Patch  Box 

(9x9  Patches  on  Top) 

NEC  Admittance  vs  "PATCH"  Code 
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Figure  25.   Monopole  at  Edge  of  Patch  Box 

(11x11  Patches  on  Top) 

NEC  Admittance  vs  "PATCH"  Code 
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3 .   Monopole  at  Corner 

The  6  cm  monopole  in  the  NEC  model  was  placed  at  the 
corner,  5.14  cm  on  the  diagonal  from  the  center  and  fed  at 
the  base.  The  position  of  the  monopole  for  the  experimental 
model  is  5.15  cm.  The  top  surface  was  divided  into  121 
(11x11)  patches.  The  calculated  values  for  the  impedance, 
admittance  and  average  power  gain  are  shown  in  Table  13. 
Experimental  data  are  presented  in  Table  2  and  the  total 
number  of  patches  as  well  as  the  patch  area  in  Table  3. 
Average  power  gain  is  between  1.8  0  (minimum)  and  1.8  3 
(maximum)  and  is  within  the  10%  tolerance  of  the  theoretical 
value. 

NEC  results  of  admittance  vs.  measurements  are 
presented  in  Figure  26.  PATCH  code  results  from  the 
calculated  data  of  Reference  5  have  been  superimposed  on 
NEC's  data,  for  comparison,  in  Figure  27.  PATCH  code  results 
are  compared  to  experimental  results  Figure  28.  NEC  vs. 
PATCH  is  compared  in  Figure  29. 

NEC  is  in  excellent  agreement  with  measurements  in 
both  conductance  and  susceptance  (Figure  26) .  PATCH  seems  to 
perform  slightly  better  than  NEC  (Figure  27) ,  however, 
conductance  and  susceptance  in  the  range  of  1.15  to  1.4  0  GHZ 
are  virtually  identical  for  both  codes. 
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TABLE  13 
CALCULATED  DATA 
MONOPOLE  AT  CORNER  OF  SURFACE  PATCH  BOX 
(5.14CM  ON  THE  DIAGONAL  FROM  CENTER) 
(11  X  11  PATCHES  ON  THE  TOP  SURFACE  -  BASE  SEGMENT  FEED) 


FREQUENCY 
(GHZ) 


AVERAGE 
POWER 
GAIN 


IMPEDANCE  (OHMS) 
ADMITTANCE  (mMHOS) 


1.00 


1.80 


36.48 
8.13 


-    J56.15 
+    J12.52 


1.05 


1.81 


35.53 
11.60 


-    J42.85 
+    J12.90 


1.10 


1.82 


40.93 
16.17 


+    D 


J29.24 
11.55 


1.15 


1.82 


44.23 
20.16 


+    D 


J15.43 
7.03 


1.175 


1.83 


46.32 
20.88 


j    8.53 
3.85 


1.2 


1.83 


48.74 
20.49 


j    1.71 
j    0.72 


1.225 


1.83 


51.50 
19.24 


+    J    4.99 
-    j    1.86 


1.3 


1.82 


61.73 
14.08 


+    D 


23.94 
j    5.46 


1.4 


1.81 


78.84 
9.52 


+    J45.42 
-    j    5.48 
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Figure  26.   Monopole  at  Corner  of  Patch  Box 

(11x11  Patches  on  Top) 

NEC  Admittance  vs  Measurements 
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Figure  27.   Monopole  at  Corner  of  Patch  Box 

(11x11  Patches  on  Top) 

NEC  Admittance  vs  Measurements  and  "PATCH"  Code 
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Figure  28.   Monopole  at  Corner  of  Patch  Box 
"PATCH"  Code  Admittance  vs  Measurements 
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Figure  29.   Monopole  at  Corner  of  Patch  Box 

(11x11  Patches  on  Top) 

NEC  Admittance  vs  "PATCH"  Code 
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E.   NEAR  ELECTRIC  FIELD 

Near  fields  are  more  difficult  to  calculate  than  far 
fields.  When  calculating  radiation  in  close  proximity  to  an 
antenna,  the  antenna  no  longer  appears  as  an  infinitesimal 
point  (as  is  the  assumption  for  far  field)  ,  and  the  energy 
radiating  from  it  comes  from  separated  locations  insofar  as 
the  strength  of  radiation  from  each  sub-section  of  the 
radiator  is  concerned.  This  means  that  terms  in  the  field 
expressions  with  powers  of  1/r2  (r  is  the  distance  from  the 
origin  of  the  antenna  to  the  field  point)  are  appreciable  in 
magnitude  compared  to  the  1/r  dependent  terms  which  are 
dominant  in  the  far  field.  The  complex  Poynting  vector 
1/2  (E  x  H*)  thus  contain  terms  with  higher  powers  of  (1/r), 
in  addition  to  the  radiation  field  term.  The  near  field  is 
thus  very  dependent  on  the  charge  density  and  the  current 
while  the  far  field  is  affected  only  by  the  current. 

For  near  field  calculations,  NPGNEC  computes  the 
magnitude  of  the  peak  electric  field  (E-Total)  in  (V/m) , 
which  is  the  vector  sum  of  the  three  components  Ex,  Ey,  and 
Ez  (the  calculation  of  Epea]<;  is  a  special  feature  that  has 
been  added  to  the  NPGNEC  version  of  NEC  which  is  in  use  at 
NPS)  .  NEC  outputs  the  magnitude  and  phase  of  each  component 
Ex,  Ey,  and  Ez  separately. 
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1.  Near  Electric  Field  in  2-Dimensions  For  the  6  cm 
Monopole  Antenna  on  the  Cubical  Conducting  Box 

Using  the  optimum  models  in  NEC  for  the  monopole  at 

the  center  (7x7  patches) ,  at  the  edge  (11x11  patches)  and  the 

corner   (11x11   patches) ,   we  proceed  with  the  near  field 

investigation.   The  NE  card  is  used  in  conjunction  with  the 

PL  card;  both  are  described  in  detail  in  Reference  1.   The 

input  data  files  are  presented  in  Appendix  A  for  the  near 

field   calculations   and   for  the  three  positions   of  the 

monopole  on  the  patch  box.   For  all  cases,  NEC  was  tasked  to 

produce  two  hundred  points  in  x,  y,  or  z  directions,  with  a 

spatial  variation  of  0.01  m.   The  starting  point  for  the  near 

field  calculation  is  0.05/}  away  from  the  wire  and  box  metal 

structure.   The  calculation  of  the  near  field  is  extended  to 

a  distance  of  at  least  6^  in  the  x,  y,  or  z  directions, 

respectively.     In  that  manner,  we  ensure  that  the  field 

points  are  far  enough  away  from  any  metallic  structure  and  at 

a  distance  where  the  near  field  variation  can  be  observed  in 

detail.   The  plots  of  the  near  field  are  presented  in  Figures 

31  and  32   (Monopole  at  center  —  variation  in  x  and  z)  , 

Figures  33,  34,  and  35  (Monopole  at  edge  —  variation  in  x, 

y,  and  z)  ,  and  Figures  36,  37,  and  38  (Monopole  at  corner — 

variation  in  x,  y,  and  z) . 

In  order  to  compare  NEC  near  fields  for  the  monopole 

on  the  box  with  theoretical  values,  we  consider  a  linear 

current  element  I=I0eJwt  of  length   Az  oriented  in  the  z 

direction  and  located  at  the  origin  as  in  Figure  30  [Ref. 
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10]  .  For  convenience,  assume  I0  is  the  amplitude  of  the 
current.  This  antenna  is  a  simple  radiating  structure  but  it 
will  demonstrate  basic  properties  of  the  near  electric  field 
for  all  small  linear  antennas. 


is: 


Figure  30.   A  Linear  Current  Radiator 
The  complete  electric-field  intensity  of  the  antenna 


-"    Vo"     V  r        r 


-kjL+J22-  +  ±-  )e-Jk^l         (2-1) 


sin# 


65 


The  only  part  of  the  field  entering  into  the  expression  for 
the  radiated  power  is  that  part  consisting  of  the  terms 
varying  as  r0_1,  that  is; 

E$=  iMgL  sin9e-;v 

This  is  called  the  radiation  field.  The  parts  of  the  field 
varying  as  r~2  and  r~3  are  called  the  induction  field.  The 
induction  field  does  not  represent  an  outward  flow  of  power, 
but  instead  depicts  a  storage  of  reactive  energy  in  the 
vicinity  of  the  radiating  current  element  (Reference  10)  . 
Consequently,  the  terms  that  are  functions  of  distance  r  of 
the  E-Field  in  Eq.  2.1  are: 


All  the  other  terms  are  phase  terms  or  constants.  Generally, 
the  magnitude  of  electric  field  has  an  r  dependence  which  can 
be  expressed  as: 


i-T-r-r 
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For  the  monopole  at  the  box  center,  the  near  electric  field 
is  presented  in  Figure  31.  From  this  figure  and  NEC  output, 
we  observe  for  a  distance  close  to  zero,  r  =  0.0166  m,  the 
value  of  the  field  is  |E(r)|  =  45.9  (V/m) .  Using  this  value 
in  Eq.  2.2,  we  can  find  the  constant  factor  C  =  0.000222065 
(normalization  factor) .  Using  the  calculated  value  of  the 
constant  C  in  Eq.  2.2  and  varying  the  distance  r  from  the 
linear  current  element,  we  calculate  the  corresponding  |E(r)| 
for  each  distance;  the  r-dependence  is  plotted  in  Figure  39. 
Comparing  Figure  31,  for  the  near  electric  field  of  the 
monopole  on  the  box  center  with  the  near  field  of  the  current 
element  in  Figure  39,  both  near  electric  fields  follow  the 
same  trend  which  is:  a  maximum  in  close  proximity  to  the 
antenna  with  a  gradual  reduction  as  the  observation  point 
moves  away.  A  knee  of  the  curve  appears  approximately  at 
0.125  m  (0.413  away  from  the  antenna).  The  results  for  the 
monopole  at  the  edge  and  corner  are  almost  the  same.  The 
"spikes"  observed  in  Figures  31,  33,  and  34  are  from  the  edge 
of  the  box. 
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Figure  31.   Monopole  at  Center  of  Patch  Box 

(7x7  Patches  on  Top) 

Near  Electric  Field  on  X  Axis 
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MAGNITUDE  VS  DISTANCE  freq=1qhz 
monopole  at  center  of  patch  box  7x7  patcnes  on  top 
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Figure  32.   Monopole  at  Center  of  Patch  Box 
Near  Electric  Field  on  Z  Axis 
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MAGNITUDE  VS  DISTANCE  freq  Ighz 
near  e-field  6  cm  monopoie  at  edge  3.63  cm  from  center 
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Figure  33.   Monopoie  at  Edge  of  Patch  Box 
Near  Electric  on  X  Axis 
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MAGNITUDE  VS  DISTANCE  freq  1  ghz 
near  e-field  6  cm  monopole  of  edge  3.63  cm  from  center 
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Figure  34.   Monopole  at  Edge  of  Patch  Box 
Near  Electric  Field  on  Y  Axis 
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MAGNITUDE  VS  DISTANCE  freq  1  ghz 
near  e-field  6  cm  monopole  at  edge  3.63  cm  from  center 
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Figure  35.   Monopole  at  Edge  of  Patch  Box 
Near  Electric  Field  on  Z  Axis 
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MAGNITUDE  VS  DISTANCE  frea  Ighi 
near  ©-field  6  cm  monopole  at  corner  5.14  cm  on  diagonal 
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Figure  36.   Monopole  at  Corner  of  Patch  Box 
Near  Electric  Field  on  X  Axis 


73 


MAGNITUDE  VS  DISTANCE  frea  Ighz 
near  e-field  6  cm  monopole  at  corner  5.M  cm  on  diagonal 
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Figure  37.   Monopole  at  Corner  of  Patch  Box 
Near  Electric  Field  on  Y  Axis 
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MAGNITUDE  VS  DISTANCE  freq  Ighz 
near  e-field  6  cm  monopole  at  corner  5.1-4  cm  on  diagonal 
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Figure  38.   Monopole  at  Corner  of  Patch  Box 
Near  Electric  Field  on  Z  Axis 
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Figure  39.   Near  El'ectric  Field  of  a  Current  Element 
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F.   NEAR  ELECTRIC  FIELD  MAGNITUDE  AND  PHASE  CONTOUR  PLOTS 

Fortran  programs  were  developed  to  plot  the  magnitude  and 
phase  of  the  near  electric  field  using  NEC  output  data.  The 
programs  use  subroutines  of  the  DISSPLA  (Display  Integrated 
Software  System  and  Plotting  LAnguage)  graphics  software 
package  [Ref.  9].  Appendix  B  lists  the  Fortran  source  codes 
and  describes  their  function.  PROGRAM  I :  Contour  plot  data 
sets  of  the  magnitude  Epeak  as  well  as  the  magnitude  of  the 
individual  components  Ex  and  Ez  of  the  electric  field  (a 
section  of  a  plane  in  3-dimensional  space) ,  composed  of  a 
30  X  30  point  window  in  space,  are  produced.  PROGRAM  II: 
This  is  similar  to  Program  I  but  for  the  phase  of  x  and  z 
components  of  the  electric  field,  for  100X100  points. 
PROGRAM  III:  Three-dimensional  plot  data  sets  of  Epea^  are 
produced.  For  these  plots,  a  3-D  map  of  the  electric  field 
is  shown  on  the  top  and  the  2-D  contours  of  the  field  are 
presented  on  the  bottom.  CONTOR  FORTRAN :  This  is  a 
subroutine  which  is  called  by  the  above  three  programs  and 
creates  the  plots.  Positive  value  contours  are  presented  by 
solid  lines  and  negative  contours  by  dashed  lines.  In  these 
programs,  NEC  output  for  the  near  electric  field  is 
interfaced  with  the  plotting  facilities  of  DISSPLA. 

For  reference,  two  simple  antennas  with  known  results  are 
examined:  (1)  a  0.15  m  (  3/2)  dipole  antenna  in  free  space, 
and  (2)  a  0.06  m  (2/5)  monopole  antenna  over  a  perfect  ground 
plane.    NEC  input  data  files  which  produce  magnitude  and 
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phase  of  near  electric  field  for  the  two  antennas  are  listed 

in  Appendix  B. 

1.   Half  Wave  Dipole  In  Free  Space  vs  6  cm  Monopole  Over 
Ground  Plane 

Near  field  plots  for  magnitude  and  phase  for  the  two 
antennas  are  presented  in  Figures  4  0  through  45.  The  near 
electric  field  is  displayed  in  a  section  of  a  plane  in  3- 
dimensional  space.  For  comparison  purposes,  the  same  view 
plane  in  the  x  and  z  directions  are  used  for  all  such  plots 
in  this  study.  The  origin  of  the  view  space  plane  is  very 
close  to  the  wire  antenna  and  box  surface  (just  0.053  away) 
and  extends  to  a  distance  of  10^  in  the  x  and  z  directions. 
For  these  two  simple  antennas,  additional  view  planes  were 
selected  in  order  to  observe  how  the  E-field  varies  away  from 
the  antenna.  For  this  case,  the  results  for  peak  E-field  for 
the  dipole  are  presented  in  Figures  40(a)  and  40(b)  (0.5^  and 
0.8^  away)  and  for  the  monopole  in  Figures  43(a),  43(b),  and 
43(c)  (1^  ,  2~X  ,  and  8^\  away).  The  electric  near  fields  are 
referenced  to  one  volt  per  meter  for  all  plots. 

It  is  well  known  that  the  A/2  dipole  has  the  same 
behavior  as  a  /1/4  monopole  (note  that  we  test  a  A/5  monopole 
in  this  study) . 

Before  comparing  the  dipole  and  monopole  field  plots, 
consider  the  following: 

ZD 
It  is  known  that  Z^  =  .   Assume  resonance  for  the  two 

2 

antennas,  then: 
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Rm  = 


Rd 


(2.3) 


The  gain  of  an  antenna  over  an  isotropic  radiator  can  be 
expressed  as: 


G  = 


;/7? 


Pr/4ttR2     Pr/47TR2 


where  S  =  power  density  (watts/m2) 
and  PR  =  power  radiated  (watts) 

Then  for  the  ^/2  dipole: 

and  for  the  monopole: 

It  is  well  known  that: 

(In  practice     GD  =  2.15  dBi  and  GM  =  5.15  dBi) 
Forming  the  ratio  of  Equations  2.4  and  2.5  and  using  Equation 
2.6,  we  get: 


ED2/7? 

GD    -                     , 
PR/47TR2 

(2.4) 

%2/7/ 

Gjyj    — 

PR/47rR2 

(2.5) 

GM    =    2GD 

(2.6) 

GM     EM' 


Thus 


Ep     1 

EM    J2 


(2.7) 


However      PD  =  ID2  Rd 
PM  =  IM2  RM 
In  order  to  compare  the  gains  of  the  two  antennas,  we  have  to 
have: 

PD  =  PM  : 


ZD2  RD  =  IM2  RM 


(2.8) 
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Using  Equation  2.3  in  Equation  2.8,  we  get 


ID2RD  =  iM2 


2 


Thus   I: 


VT 


and 


ID   _  *_ 
IM    J2 


From  Equations  2.7  and  2.9,  we  conclude  that: 


EM 


■M 


Thus,  the  electric  fields  will  be  equal  only  if  the  currents 
are  equal  for  the  two  antennas. 

In  Table  14,  characteristics  of  the  two  antennas  from 
the  NEC  output  files  are  presented.  From  Figures  40  and  43, 
we  conclude  that  the  shape  of  near  electric  field  is  the  same 
for  the  dipole  and  monopole  antenna.  The  ratio  of  the  input 
maximum  current  for  the  two  antennas  is: 
0.9764  X  10"2 


20  log10 


■7.28  dB 


0.2259  X  10"1 

From  Figures  40  and  43,  we  see  that  the  0  dB  contours  for  the 
dipole  are  at  a  distance  2"X  away  on  the  x-axis  and  for  the 
monopole  at  a  distance  of  3}  on  the  same  axis.  It  appears 
that  the  field  attenuates  more  rapidly  with  radial  distance 
away  from  the  dipole  as  compared  to  the  monopole.  For  both 
antennas,  the  maximum  value  appears  at  a  40°  angle  from  the 
x-axis.  As  the  observation  point  moves  away  from  the  antenna 
along  the  x-axis,  there  is  a  difference  of  4.6  dB  between  the 
electric   field  values   on  the  contour  plots  of  the  two 
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antennas.  Also,  note  that  as  the  observation  point  moves 
away  from  the  box  surface  along  the  z-axis,  the  difference 
appears  to  be  5.6  dB.  Compared  to  the  theoretical  difference 
of  7.28  dB  that  we  expect,  there  is  a  2  dB  offset  which  is  a 
result  of  comparing  a  A/5  monopole  instead  of  a  A/A  one  to 
the  A/2    dipole. 

The  phase  plots  of  the  Ex  and  Ez  components  of  the 
electric  fields  for  the  A/2  dipole  and  0.06  m  monopole  in 
Figures  41,  42,  44,  and  45  show  that  both  antennas  have  a 
smooth  spherical  wavefront  pattern. 
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TABLE  14 
RESULTS  FROM  NEC  OUTPUTS 


DIPOLE  V2 

IN  FREE 

SPACE 

(5-SEGMENTS) 

6  CM  MONOPOLE 
OVER  PERFECT 

GROUND 
(5-SEGMENTS) 

6CM  MONOPOLE 

MOUNTED  ON 

SURFACE  PATCH 

BOX  CENTER 

(5-SEGMENTS) 

LENGTH 
(meters) 

0.15m  (  3/2) 

0.06m  (  3/5) 

0.06m  (  3/5) 

WAVELENGTH 
(meters) 

0.  3m 

0.3m 

0.3m 

EXCITATION 

(voltage 

source) 

1  Volt 

1  Volt 

1  Volt 

FEED 

POSITION 
(#  OF 
SEGMENT) 

3rd 

1st 

1st 

MAXIMUM 
CURRENT 
(magnitude 
in  amps) 

0.98  X  10"2 

0.22  X  10"1 

0.26  X  10"1 

IMPEDANCE 
( ohms ) 

89.61  +  J49.57 

26  -  J35.82 

10.99  -  J37.53 

INPUT 
POWER 
(watts) 

0.43  x   10"2 

0.66  X  10"2 

0.36  X  10"2 
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CONTOUR   E-PIOLD    (OB   REFER   TO    1V/M 
DIPOLE   LENGTH  X/2   ON   Z   RXIS 
FREE   SPRCE-FRE0=1GHZ 


0.0166     0.5966       1.1766       1.7566       2.3366 

X-RXI5  VRRIRTION    (METERS) 


Figure   40.      Total   E-Field   Contours 
Dipole      "XJ2   on    Z   Axis    in    Free    Space 
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CONTOUR   E-FIELD    (DB   REF   TO    1V/M) 

DIPOLE   LENGTH   X/2   ON   Z-nXIS 

FREE   SPRCE-FREQ=1GHZ 


0.04*  0.073  0.102  0.131  0. 

X-RXIS  VRRIRTION    (METERS) 


Figure    40(a).      Total   E-Field   Contours 
Dipole        A/2    in    Free    Space    (0.5  ^    away    from   the   wire) 
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CONTOUR   E-FIELD    (DB   REF   TO    1V/M) 

DIPOLE   LENGTH  A/2   ON   Z-RXIS 

FREE   SPRCE-FRE0=1GHZ 


0.0150     0.0614       0.1078       0.1542       0.2006       0.2470 

X-RXI5  VARIATION    (METERS) 


Figure   40(b).    Total   E-Field   Contours 
Dipole       A/2    in    Free    Space    (0.8  3    away    from   the   wire) 


85 


FHRSE   OF   X   COMPONENT   OF   E-FIELD 
DIPOLE  X/2   ON   Z   RXIS 
„     IN   FREE   SPRCE-FRE0=1GHZ 


0.01S6     0.5908       1.1650       1.7392       2.3134       2.8875 

X-flXIS  VRRIflTION    (METERS) 


Figure   41.       X-Component ,    E-Field   Phase   Contours 
Dipole      /\/2    in   Free    Space 
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DUQC 


OF   Z   COMPONENT   OF   E~F 
CIPOLE  \/2   ON   Z   RXIS 
IN   FPEE   SPflCE-FREO=iGHZ 


n 


D.0166     0.5908 


1.1650       1.7392       2.3134 

X-flXIS  VRRIflTION    (METERS) 


Figure    42.       Z-Component ,    E-Field   Phase    Contours 
Dipole      A/2    in    Free    Space 
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CONTOUR   E-FIELQ  (DB   REFER   TO    1V/M) 

MONOPOLE   6   CM   ON   Z   RXIS 

OVER   PERFECT   GR0UND-FREQ=1GHZ 

in 

n 

j\\  /    /      &           /                        X''' 
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CI 

0.0166            0.5966                 1.1766                 1.7566                 2.3366                 2.S 

3166 

X-RXI5  VARIATION    (METERS) 

Figure  43.   Total  E-Field  Contours 
Monopole  6cm  Over  Perfect  Ground 
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CONTOUR   E-FIELD    (DS   REP   TO    1V/M) 

MONOPOLE   6   CM   ON   Z-flXIS 

O^ER   PERFECT   GR0UND-FREQ=1GHZ 

d 

03 
rg- 

O 

CO 
LJ    CM 

5  2" 

O 

E— 
CE 

cr 

>    CM 
01 

in  "t" 

—  o 

X 

CE 
M 

d 

§ 

o 

^^^                   /         /                %'                         /^ 

mv     ^\   X 

0.0166            0.0530                 0.1094                 0.1558                 0.2022                 0. 

X-RXI5  VflRIRTION    (METERS) 

H86 

Figure  43(a).   Total  E-Field  Contours 
Monopole  6cm  Over  Perfect  Ground  (13  away  from  the  wire) 
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CONTOUR   L-FIELD    (DB   REF   TO    1V/M 
MONOPOLE   6   CM   ON   Z-RXIS 
OVER   PERFECT   GROUND-FREO=lGHZ 
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\      \ 

0.0166     0.1326       0.2486       0.3646       0.4806 

X-RXI5  VRRIflTION   (METERS) 


Figure   43(b).      Total   E-Field  Contours 
lonopole    6cm  Over   Perfect    Ground    (2  a     away    from   the   wire) 
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CONTOUR   E-FIELD    (DB   REF   TO    1V/M) 

MONOPOLE   6   CM   ON   Z-RXIS 

OVER   PERFECT   GR0UND-FREQ=1GHZ 

r\i 

Ci_ 

;'' ;'  /   /     /        ^'                         ..—--11.77- 

^, 

;;'?;     '      /           ,\v 

(T) 

.->  tn  i    , 

DC 

r\j ;      ,         / 

LJ  ^_ 

z:  j 

2 
CD 

■-'//»  /  /       _..--8-« 

i— 

'  i N"  -''    /'         ,''"'                                   ~**> 

cr 

'  /  y     '       /  •          ,'" 

2 

'  '  /  -'   A 

>  £ 

;  /  <    V       / 

^-" 

i  ,'  ,    ^," 

X 

■'//  '**<£         ..-""" ~"^*n 

M 

-■!//    # 

81 

-'!//  *.* --... 

!^ 

IT  _ 

'~'it/                               \                         \ 

d 

~  'f-      ^                  \                    ^ 

i/i 

^X N              *\                    \                              '.                                         CD 

d_ 

N> .        N            \                  \                          \ 
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X-RXI5  VRRIRTION    (METERS) 

Figure  43(c).   Total  E-Field  Contours 
Monopole  6cm  Over  Perfect  Ground  (83   away  from  the  wire) 
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PHRSE   OF   X   COMPONENT   OF   E-FIELD 

NONOPOLE   6CM   ON   Z   AXIS 

OYER   PERFECT   GR0UNQ-FREQ=1GHZ 

S  


0.0166     0.5908       1.1650       1.7392       2.3134 

X-RXIS  VRRIRTION    (METERS) 


Figure   44.      X-Component ,    E-Field  Phase   Contours 
Monopole   6cm  Over   Perfect    Ground 
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PHPSE   Of 


r    n 


COMPONENT   OP   E-PIELD 


MONOPOLE   6CM   ON   Z   RXIS 
ER   PERPECT   GROOND-PREQ=lGHZ 


0.0166     0.5908       1.1650       1.7392       2.3134 

X-RXI5  VARIATION  (METERS) 


2. 8876 


Figure  45.   Z-Component ,  E-Field  Phase  Contours 
.Monopole  6cm  Over  Perfect  Goound 
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2 .   Monopole  at  Center  of  Patch  Box  vs  Monopole  Over 
Ground  Plane 

Input  data   files   for  NEC  solutions  of  the  near 

electric   field   of   the  monopole   on   the   box   center   are 

contained  in  Appendix  B.   In  Figures  46  to  50,  the  magnitude 

of  the  total  electric  field  and  the  magnitude  and  phase  of 

the  Ex  and  Ez  components  are  shown.    In  Figure  46,  the  E- 

field  zero  dB  contour  is  observed  to  be  closer  to  the  antenna 

(it  is  at  2?     vs.  3^)  than  for  the  monopole  over  a  ground 

plane.    Using  the  ratio  of  the  maximum  input  currents  alone, 

we  calculate  that  the  theoretical  difference  in  the  electric 

fields  of  the  two  antennas  should  be  -1.08  dB.   The  contour 

plots  of  Figures  46  and  43  show  that  there  is  a  difference  of 

-3  dB  along  the  box   surface   direction   and  a   difference  of 

+8  dB  up  the  z-axis,  perpendicular  to  the  surface,  compared 

to  the  isolated  monopole  reference  case.   These  differences 

in  near  field  densities  are  attributed  to  box  radiation  and 

diffraction.   In  this  near  field  region  we  can  observe  beams 

and  nulls  forming.   A  maximum  occurs  at  60°  up  from  the  box 

surface  vs.  4  0°  for  the  monopole  over  the  ground  plane.  The 

main   lobe   (Figure  46)   starts  to  develop   at  a   distance  2} 

(0.6  m)  from  the  antenna,  along  the  x-axis.    Beyond  this 

point,   it  generally  has  the  same  pattern  independent  of 

distance.   Within  the  2j\   distance,  the  pattern  shape  has  not 

yet  fully  developed. 
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In  Figures  47  and  49,  Ex  and  Ez  contour  plots 
demonstrate  the  behavior  of  the  magnitude  (the  Ey  component 
is  negligible  and  is  not  presented) . 

Recalling  the  Ep^^  field  variations  and  noting  that 
Epeak  consists  of  Ez  and  Ex  components,  we  can  scan  Figures 
47  and  49  and  see  how  these  components  contribute  to  the 
total  field  of  Figure  46.  Different  variations  with  distance 
along  the  x-  and  z-axes  and  with  take-off  (elevation)  angle 
are  noted.  The  z-component  is  predominant,  as  it  should  be, 
for  a  z-directed  monopole.  Null  effects  are  most  apparent 
for  Ez  at  moderate  elevation  angles  and,  as  expected,  Ex 
nulls  at  the  zenith. 

Figures  48  and  50  show  the  phase  of  the  x  and  z 
components.  The  electric  field  phase  has  the  same  spherical 
wavefront  pattern  as  the  monopole  over  the  ground  plane 
(Figures  44  and  45),  except  for  the  region  of  the  null, 
caused  by  the  box.  The  disturbance  in  the  electric  field 
phase  appears  at  the  same  27°  location  in  both  phase  plots 
(Figures  48  and  50)  .  It  is  guite  severe  in  the  z  component 
phase  (Figure  50)  and  is  barely  observed  in  the  x  component 
phase  (Figure  48) . 

In  order  to  get  better  insight  into  the  electric  near 
field  variations  and  where  the  maxima  and  minima  occur, 
3-dimensional  plots  are  presented  in  Figures  51  and  52.  These 
plots  were  produced  from  Fortran  program  III,  described  in 
Appendix  B  and  use  the  same  view  plane  as  previous  plots. 
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Each  plot  displays  a  surface  whose  elevation  points 
represents  field  strength  in  the  upper  portion,  and  a  normal 
2-D  contour  plot  "projection"  in  the  lower  portion.  Two 
different  view  points  were  selected,  one  "looking  away  from" 
the  monopole  and  one  towards  it.  In  each  figure  the  vertex 
that  corresponds  to  the  pointed  "spike-like"  area  of  the 
surface  is  the  origin  where  the  antenna  is  mounted.  The 
delay  of  the  field  as  the  observation  point  moves  away  from 
the  monopole  source  is  expected.  The  null  "trough"  is  very 
striking  in  this  type  of  presentation. 
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CONTOUR   E-FIELD    (DB   REF   TO    1V/M) 

MONOPOLE   6   CM   AT   CENTER   OF   SPBOX 

7X7   PRTCHES   ON   TOP-FREQ-IGHZ 
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0.0166     0.5966       1.1766       1.7566       2.3366       2. 915 

X-RXIS  VRRIRTION   (METERS) 


Figure   46.      Total   E-Field   Contours 
Monopole    at   Patch   Box   Center 
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MRSNITUDE    (EX   PERK)    OF   E-PIELD 

MONOPOLE   6CN   RT   CENTER   OP   SPBOX 

ZX7   PATCHES   ON   T0P-PREQ=1GHZ 
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Figure  47.   X-Component ,  E-Field  Contours 
Monopole  at  Patch  Box  Center 
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PHRSE   OP   X   COMPONENT   OP   E-PIELD 

MONCPOLE   6CM   RT   CENTER   OP   SP30X 

1X7   PRTCHES   ON   T0R-PRE0=1GHZ 


0.0166 


1.16S0       1.7392 

•AXIS  VARIATION  (METERS! 


2.3134 


2.8876 


Figure  48.   X-Component ,  E-Field  Phase  Contours 
Monopole  at  Patch  Box  Center 
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MRGNITUEE    (EZ   PERK)    OF   E-riELO 
MONOPOLE   6CM   PT   CENTER   OF   SPBOX   | 
7X7   PRTCHES   ON   T0P-FREQ=1GPZ 


/  /    # 

°&* 

i           !             ^ 

'           / 

b- 

*?                                  & 

$        / 

'       $          / 

- -' 

0)' 

^-' 

r* 

/ 

/ 

<bw 

/ 

/      «$>'  ,-"'>: 

/                 ,''                 \ 

/     ,^'  x'V;:^ 

/    $ 

'''   /  A-^&Pv 

"----'' 

'■^.,y 

'.-■'' 

;...-i4.62- 

~\             /'      / ,','.i\'.''  '.-''"' 

.. 

)        '    /''*'. ''»»'                



h)/ 

^V5:!l"/.'.'."."" 

■*-s 

0.0165     0.5908       1.1650       1.7392       2.3134 

X-RXIS  VRRIRTION    (METERS) 


Figure   49.       Z-Component ,    E-Field   Contours 
Monopole    at   Patch   Box   Center 
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m  ,i^i   u~    Z   oOMrul\LNi    U\     i    I-  i  l L U 

MONQPOLE   6CN   RT   CENTER   Or   SPBGX 

7X7   PnTCHES   ON   TOD-FREQ=i GHZ 


1.1650  1.7392 

5  VRRIRTION    (METERS! 


2.8876 


Figure  50.   Z-Component ,  E-Field  Phase  Contours 
Monopole  at  Patch  Box  Center 
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Figure  51.  Total  E-Field  3-D  Plot,  View  Toward  Monopole 
(Monopole  at  Patch  Box  Center) 
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Figure  52.   Total  E-Field  3-D  Plot,  Viewed  from  Monopole 
(Monopole  at  Patch  Box  Center) 
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3.  Monopole  at  Edge  of  Patch  Box  vs.  Monopole  at  Center 
Figures  53  to  59,  for  the  edge-mounted  monopole 
follow  the  same  format  and  trend  as  the  center-mounted  case. 
The  following  comments  will  amplify  differences  from  the 
center  case  and  interesting  features  of  the  field 
distributions  (NEC  input  data  sets  are  in  Appendix  B) . 

The  edge-mounted  geometry  provides  a  larger  planar 
surface  (box  top)  in  front  of  the  monopole  view  plane  and 
results  in  the  following  differences  with  respect  to  the 
center-mount  geometry: 

The  z-axis   peak   field   is  somewhat   greater   at   given 

distances   from  the  box   surface.   This   is  traced  to  a 

larger  Ex  component. 

The  elevation  plane  null  is  not  as  deep. 

The  phase  contours  for  the  z  component  shown  evidence  of 

two   close-in   nulls,  but  very   little  "phase  wrinkles" 

beyond  the  near  zone. 

The   contour  plots  depict  a  more  uniform  field  overall, 

with  a  less  severe  null. 
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CONTOUR   l-FIlLQ    (DB   RT~   TO    1V/M) 

MO.\0?OlE  6CM  RT  EDGE (3. 63  CM  FROM  CENTER!  CR  SP5CX 
11X11  PfiTCHES  ON  TOP-FREQ-iGHZ 
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Figure    53.      Total   E-Field   Contours 
Monopole    at    Patch   Box   Edge 
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Figure    54.      X-Component ,    E-Field   Contours 
Monopole    at    Patch   Box   Edge 
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Figure    55.      X-Component ,    E-Field  Phase   Contours 
Monopole    at   Patch   Box   Edge 


107 


P 0 L£   6 C M   RTEDGe'oP~SPB OX 

11X11    PnTCHPS   ON   IOP-PREG= 1GHZ 


&2T 


oV 


^ 


.  -\1 


■14.19 


-0.31S75    0.5=444      1.12854      1.70284      2.27703 

X-RXIS  VflRIfiTION    (METERS! 


Figure    56.       Z-Component ,    E-Field   Contours 
Monopole    at   Patch   Box   Edge 
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Figure  57.   Z-Component ,  E-Field  Phase  Contours 


Monopole  at  Patch  Box  Edge 
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Figure  58.   Total  E-Field  3-D  Plot,  View  Toward  Monopole 
(Monopole  at  Patch  Box  Edge) 
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Figure  59.   Total  E-Field  3-D  Plot,  Viewed  From  Monopole 
(Monopole  at  Patch  Box  Edge) 
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4 .  Monopole  at  Corner  of  Patch  Box  vs  Monopole  at  Center 
The  corner-mounted  monopole  near  field  plots  are  in 
Figures  60  to  66.  The  NEC  input  data  is  in  Appendix  B.  The 
fields  for  this  case  fall  in  between  the  center  and  edge- 
mounted  field  configurations.   That  is: 

The   elevation  plane   null   is   between   the  other   two 

configurations'  nulls. 

Phase  wrinkles  show  one  null,  as  does  the  center-mounted 

case,  but  it  is  not  as  severe. 

Contour   plots  indicate  a  "flatter"  field  pattern  with  a 

fairly  uniform  distribution  away  from  the  monopole. 
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Figure    60.      Total   E-Field   Contours 
Monopole    at    Patch   Box   Corner 
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Figure  61.   X-Component ,  E-Field  Contours 
Monopole  at  Patch  Box  Corner 
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Figure    62.       X-Component ,    E-Field   Phase   Contours 
Monopole    at    Patch   Box   Corner 
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Figure   63.       Z-Component ,    E-Field   Contours 
Monopole    at   Patch   Box   Corner 
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Figure   64.       Z-Component ,    E-Field  Phase   Contours 
Monopole    at    Patch   Box   Corner 
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Figure  65 


Total  E-Field  3-D  Plot,  View  Toward  Monopole 
(Monopole  at  Patch  Box  Corner) 
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Figure  66.   Total  E-Field  3-D  Plot,  Viewed  From  Monopole 
(Monopole  at  Patch  Box  Corner) 
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G.   RADIATION  PATTERNS 

Radiation  patterns  were  calculated  and  are  presented  in 
Figures  67  and  68  (Monopole  at  center)  ,  69  -  71,  (Monopole 
at  edge) ,  and  72  -  74  (Monopole  at  corner) . 

1.  Monopole  at  Center 

In  Figure  67 ,  the  vertical  pattern  shows  that  the 
maximum  gain  is  very  close  to  5.15  dBi,  the  theoretical  value 
for  a  monopole  over  an  infinite  perfectly  conducting  ground 
plane.  The  metallic  surface  of  the  box  which  re-radiates 
affects  the  pattern  causing  four  lobes  to  appear  instead  of 
the  usual  two.  In  Figure  68,  the  horizontal  pattern  shows 
omnidirectional  radiation  from  the  box-monopole,  which  is 
fairly  small  electrically,  and  not  expected  to  contribute 
much  directionality  in  azimuth. 

2 .  Monopole  at  Edge 

The  vertical  pattern  (x-axis  cut)  in  Figure  69  shows 
an  unsymmetrical  pattern  with  the  maximum  lobe  of  6.6  dBi 
skewed  away  from  the  monopole.  The  y-axis  cut  in  Figure  7  0 
shows  an  expected  symmetric  vertical  plane  pattern.  The 
horizontal  pattern  of  Figure  71  shows  a  skewed  pattern,  with 
the  maximum  field  opposite  to  the  monopole  location. 

3 .  Monopole  at  Corner 

The  results  of  vertical  and  horizontal  patterns  for 
a  corner-mounted  monopole  (Figures  72  -  74)  are  similar  to 
the  ones  for  the  edge-mounted  case  with  unsymmetrical 
patterns  occurring. 
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Figure  67.   Vertical  Pattern,  Monopole  at  Patch  Box  Center 
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SPBOX  7X7  PATCHES    ON  THE  TOP  (MONOPOLE  AT  CENTER) 
HORIZONTAL  PATTERN 
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Figure  68.   Horizontal  Pattern,  Monopole  at  Patch  Box  Center 
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Figure  69.   Vertical  Pattern  (X-Axis  Cut) 
Mono-pole  at  Patch  Box  Edge 
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Figure  70.   Vertical  Pattern  (Y-Axis  Cut) 
Monopole  at  Patch  Box  Edge 
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SPBOX  11X11  PATCHES    ON  THE  TOP  (MONOPOLE  AT  EDGE  3.63  CM) 
HORIZONTAL  PATTERN 
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Figure  71.   Horizontal  Pattern 
Monopole  at  Patch  Box  Edge 
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SPBOX  11X11  PATCHES  (MONOPOLE  AT  CORNER-5.UCM  ON  DIAGONAL) 
VERTICAL  PATTERN  (X-AXIS  CUT) 
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Figure  72.   Vertical  Pattern  (X-Axis  Cut) 
Monopole  at  Patch  Box  Corner 
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SPBOX  11X11  PATCHES  (MONOPOLE  AT  CORNER-5.UCM  ON  DIAGONAL) 
VERTICAL  PATTERN  (45  DEGREES  CUT) 
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Figure  73.   Vertical  Pattern  (45   Cut) 
Monopole  at  Patch  Box  Corner 
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SPBOX  11X11  PATCHES  (MONOPOLE  AT  CORNER-5.UCM  ON  DIAGONAL) 
HORIZONTAL  PATTERN 
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Figure  74.   Horizontal  Pattern 
Monopole  at  Patch  Box  Corner 
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III.   WIRE  GRID  MODELING  FOR  NEAR  ELECTRIC  FIELD  PREDICTIONS 

A.  GUIDELINES  FOR  WIRE  GRID  MODELING  USING  NEC 

Solid  surfaces  can  be  modeled  in  NEC  with  a  grid  of 
wires,  with  the  restriction  that  the  grid  cells  are  to  be 
small  in  terms  of  a  wavelength.  Wire  grid  modeling 
guidelines  are  given  in  Reference  1. 

B.  WIRE-GRID  BOX  MODEL  FOR  NEAR  ELECTRIC  FIELD  INVESTIGATION 
For  the  wire  grid  modeling  technique,  the  NPG1000  version 

of  NEC  allows  up  to  1000  wire  segments.  Typical  run-times 
for  wire  grid  models  of  the  box-monopole  are  six  times  those 
of  the  surface  patch  models.  The  box  of  Figure  3  is  modeled 
as  a  five-sided  wire-grid  box  of  0.1  m  (  -y3  at  1  GHz)  per 
side  and  cells  of  0.0125  by  0.0125  meters.  The  0.06  m 
monopole  antenna  was  divided  into  5  segments  and  placed  on 
top  of  the  wire  grid  box  at  the  center,  edge  (3.75  cm  from 
center),  and  corner  (5.3  cm  on  the  diagonal).  The  geometry 
of  the  wire  grid  box  is  shown  in  Figures  75  -  77.  Appendix  C 
contains  NEC  input  data  files  for  these  geometries.  The 
antenna  was  fed  at  the  base  segment  for  all  cases.  The  wire 
grid  geometry  of  Reference  8,  was  used  for  calculations  of 
near  electric  field  in  the  present  study.  This  geometry 
produced  good  results  for  admittance  and  average  gain 
compared  with  experimental  data  of  Reference  5.  The  Fortran 
programs  of  Chapter  II  are  also  used  for  magnitude  and  phase 
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contour  plots,  and  3-D  plots  of  near  electric  field  for  the 
wire  grid  box  case  with  the  same  field  point  locations  used 
in  the  surface  patch  model. 
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MONOPOLE  6  CM  AT  CENTER  OF  WIRE  GRID  BOX 


THETA  =  60.00  PHI  =  30.00  ETA  =  90.00 


Figure  75.   Monopole  at  Center  of  Wire  Grid  Box 
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MONOPOLE  AT  EDGE  (3.75  CM)  OF  WIRE  GRID  BOX 


THETA  =  60.00    PHI  =  30.00    ETA  =  90.00 
Figure    76.      Monopole    at    Edge   of   Wire   Grid   Box 
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MONOPOLE  6  CM  AT  CORNER  (5.3  CM)  OF  WIRE  GRID  BOX 


THETA  =  60.00    PHI  =  30.00    ETA  =  90.00 
Figure    77.      Monopole    at    Corner   of   Wire   Grid   Box 
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1.   Monopole  at  Center  of  Wire  Grid  Box  vs.  Monopole  at 
Center  of  Surface  Patch  Box 

The  excellent  agreement  of  near  fields  (Figures  78  to 

84)   for  the  wire  grid  box  with  center-mounted  monopole  to 

those  of  the  surface  patch/center  monopole  model  (Figures  4  6 

to  52)   attests  to  the  equivalence  of  the  two  numerical 

models.    Differences  are  less  than  1  dB,  a  value  which  is 

difficult  to  measure. 
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Figure    78.      Total   E-Field   Contours 
Monopole    at   Wire   Grid   Box   Center 
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Figure    79.      X-Component ,    E-Field   Contours 
Mono.pole    at   Wire   Grid   Box   Center 
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Figure  80.   X-Component ,  E-Field  Phase  Contours 
Monopole  at  Wire  Grid  Box  Center 
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Figure    81.       Z-Component ,    E-Field   Contours 
Monopole    at   Wire   Grid   Box   Center 
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Figure  82.   Z-Component ,  E-Field  Phase  Contours 
Monopole  at  Wire  Grid  Box  Center 
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Figure  83.   Total  E-Field  3-D  Plot,  View  Toward  Monopole 
(Monopole  at  Wire  Grid  Box  Center) 
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Figure  84.   Total  E-Field  3-D  Plot,  Viewed  From  Monopole 
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2 .   Monopole  at  Edge  of  Wire  Grid  Box  vs.  Monopole  at 
Edge  of  Surface  Patch  Box 

The  edge-mounted  monopole  for  the  wire  grid  model 

compares  well  to  the  surface  patch  version.   The  differences 

are  expected  to  be  attributed  to  a  more  accurate  rendition  of 

edge  effects  for  the  grid  model  versus  the  patch  model. 

Observable  differences  between  Figures  85  to  91  and  53  to  59 

relate  to  deeper  nulls  for  the  wire  grid  model. 
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Figure    85.      Total    E-Field   Contours 
Monopole    at    Wire   Grid  Box   Edge 
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Figure   86.      X-Component ,    E-Field   Contours 
Monopole    at   Wire   Grid   Box   Edge 
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Figure    87.      X-Component ,    E-Field  Phase   Contours 
Monopole    at    Wire   Grid   Box   Edge 
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Figure    88.       Z-Component ,    E-Field    Contours 
Monopole    at    Wire   Grid   Box   Edge 
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Figure    89.       Z-Component ,    E-Field   Phase   Contours 
Monopole    at   Wire   Grid   Box   Edge 
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Figure   90 


Total  E-Field  3-D  Plot,  View  toward  Monopole 
(Monopole  at  Wire  Grid  Box  Edge) 
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Total  E-Field  3-D  Plot,  Viewed  From  Monopole 
(Monopole  at  Wire  Grid  Box  Edge) 
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3 .   Monopole  at  Corner  of  Wire  Grid  Box  vs.  Monopole  at 
Corner  of  Surface  Patch  Box 

The  previous  two  correlations  between  surface  patch 

box  and  wire  grid  models  continue  for  the  corner-located 

monopole.     Values   of   1   to   1.5   dB   in   field   strength 

differences,   with  deeper  nulls,   summarize  the  results  of 

comparing  Figures  92  to  98  with  Figures  60  to  66. 
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Figure    92.      Total   E-Field   Contours 
Monopole    at    Wire   Grid   Box   Corner 
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Figure  93.   X-Component ,  E-Field  Contours 
Monopole  at  Wire  Grid  Box  Corner 
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Figure    94.      X-Component ,    E-Field  Phase   Contours 
Monopole    at   Wire   Grid   Box   Corner 
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Figure   95.       Z-Component ,    E-Field   Contours 
Monopole    at   Wire   Grid   Box   Corner 
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Figure   96.       Z-Component ,    E-Field  Phase   Contours 
Monopole    at    Wire   Grid   Box   Corner 
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Figure  97.   Total  E-Field  3-D  Plot,  View  Toward  Monopole 
(Monopole  at  Wire  Grid  Box  Corner) 
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Figure  98.   Total  E-Field  3-D  Plot,  Viewed  From  Monopole 
(Monopole  at  Wire  Grid  Box  Corner) 
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IV.   CONCLUSIONS  AND  RECOMMENDATIONS 

The  goal  of  this  investigation  was  to  accurately  predict 
near  fields  using  the  Numerical  Electromagnetics  Code  (NEC) 
on  ship-like  structures.  The  near  fields  of  Navy  HF  antennas 
are  primary  contributors  to  HF  radiation  hazards  (RADHAZ) . 
Providing  reliable  interference-free  communications  is  a 
formidable  task  because  of  the  high  density  of 
electromagnetic  sources  on  the  topside  of  a  modern  warship. 

Since  no  validation  benchmark  results  for  near  fields 

were  available  for  ship  antenna  installations,  a  modeling 

exercise  was  undertaken  using  NEC  in  order  to  select  the 

optimum  models  which  would  be  used  to  calculate  near  fields. 

Box-like  structures  were  analyzed  using  surface  patch  and 

wire  grid  modeling  techniques.   The  simulation  models  consist 

of  a  /i/5  monopole  mounted  on  the  top  of  a  A/ 3    box  at  three 

different  locations:  center,  edge  and  corner.   Optimum  models 

were  selected  by  varying  the  patch  density  on  the  top  surface 

of  the  box  and  observing  the  convergence  of  the  solution  for 

input  admittance  and  average  power  gain,  in  order  to  ensure 

the  validity  of  the  models  and  improve  the  confidence  in  near 

field   predictions.     Optimum   NEC  models   for   the   three 

different  mounting  geometries  are: 

2 

-  CENTER:  7x7  =  49    patches   on   top   (.0016^  patch  area) 

-  EDGE  AND  CORNER:  11x11  =  121    patches   on   top  ( . 0009^2 
patch  area) 
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NEC  results  for  input  admittance  correlate  very  well  with 
measurements.  Average  power  gain  varied  between  1.80  to  2.02 
(within  the  10%  engineering  tolerance  of  the  theoretical 
value) .  Single  Precision  (NPGNEC)  calculations  gave  the  same 
results  as  Double  Precision  (DNPGNEC) . 

Even  though  edges  are  not  modeled  in  the  surface  patch 
technique  [Ref .  1]  ,  this  study  proves  that,  for  positions 
very  close  to  an  edge,  good  results  can  be  obtained  by 
careful  subdividing  (no  special  subdivision  of  smaller 
patches  in  the  vicinity  of  the  edge  or  corner  was  required) . 

Fortran  algorithms  were  developed  to  produce  near 
electric  field  (magnitude  and  phase)  contours  and  3-D  plots 
using  DISSPLA  [Ref.  9].  The  near  field  for  the  monopole  on 
the  box  has  the  same  characteristics  in  magnitude  and  phase 
as  the  monopole  over  a  ground  plane  except  in  a  region  where 
nulls  occur  from  box  radiation  and  diffraction  effects.  The 
edge/corner-mounted  geometries  produce  slightly  different 
near  electric  field  maps  compared  to  center-mounted  geometry. 

Surface  patch  and  wire  grid  models  for  NEC  gave  similar 
results  for  near  fields. 

Previously,  generalized  guidelines  for  near  field 
modeling  had  not  been  developed  for  NEC  and  the  use  of  wire 
grid  and  surface  patch  modeling  for  near  field  parameters  was 
approached  with  caution.  Guidelines  developed  in  this  study, 
as  well  as  the  results  of  the  near  field  behavior  of  the 
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monopole  antenna  on  the  conducting  box,   can  be  used  for 
future  investigations  on  "ship-like"  structures. 

Recommendations  for  continuation  of  the  present  study 
are: 

Extend  the  model  to  larger  and  more  complex  structures 
which  more  closely  approximate  the  ship's  topside  super- 
structure. 

Develop  additional  Fortran  algorithms  for  use  in  near 
field  investigations. 

Exercise  a   new  version   of  the   PATCH  Code,  "JUNCTION" 

[Ref.  11]  in  order  to  check  the  validity  of  the  code 

and  the  computer/time  cost  in  comparison  to  NEC. 

The  present  study  is  an  important  step  in  the  direction 

of  modeling  the  effects  of  shipboard  environments  on  the  near 

field  of  antenna  structures.   It  demonstrates  how  NEC  surface 

patch  and  wire  grid  modeling  techniques  can  be  applied  for 

accurate   near   field  predictions   of   antennas   in  complex 

environments. 
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APPENDIX  A 
NEC  INPUT  DATA  FILES 


PART  A:   NEC  INPUT  DATA  FILES  FOR 
SURFACE  PATCH  MODELING 


FILE:  SP3OX100  DATA      Bl 


CM  SP30X   5X5  PATCHES  ON  TOP  SURFACE 

Cz    FREQUENCY  1.0  0  GHZ 

SMS, 3, 0.5,-0.5,0,0.5,-0.5,1 

SCO, 0,0, -0.5,1 

SMS,  5, 0.5,-0.5,0,0.5,0,0 

SCO, 0,0. 5, 0,1 

GX  0,110 

SM5, 5, -0.5, -0.5, 1,0. 5, -0.5,1 

SZ3, 0,0.5,0.5,1 

GS  0,0,0.1 

GE  1 

FR    0,0,0,0,1000 

GN  1 

NG 

XC  0 

EN 


1.   Data  File  Surface-Patch  Box  Structure 
(5x5  patches  =  25  patches  on  top 
surface ) . 


FILE:  SP3CX1    DATA     Ml 


CM  SP30X   7X7  PATCHES  ON  TOP  SURFACE 

CE    FREQUENCY  1.00  GHZ 

SMS, 3, 0.5, -0.5, 0,0. 5, -0.5,1 

SCO, 0,0, -0.5,1 

SMS, 3, 0.5, -0.5, 0,0. 5,  0,0 

SCO, 0,0. 5, 0,1 

GX  0,110 

SM7, 7, -0.5, -0.5, 1,0. 5, -0.5,1 

SCO, 0,0. 5, 0.5,1 

GS  0,0,0.1 

CE    1 

FR  0,0,0,0,1000 

GN  1 


WG 

XO  0 
EN 


Data  File  Surface-Patch  Box  Structure 
(7x7  patches  =  49  patches  on  top 
surface ) . 
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FILE:  SPB1      DATA     Ml 

CM  SP30X   9X9  PATCHES  ON  TOP  SURFACE 

CE  FREQUENCY  1.00  GHZ 

SM3, 3, 0.5, -0.5, 0,0. 5, -0.5,1 

SCO, 0,0, -0.5,1 

SM3, 5, 0.5, -0.5, 0,0. 5, 0,0 

SCO, 0,0. 5, 0,1 

GX  0,110 

SM9, 9, -0.5, -0.5, 1,0. 5, -0.5,1 

SCO, 0,0. 5, 0.5,1 

GS  0,0,0.1 

GE  1 

FR  0,0, 0, 0,1000 

GN  1 

UG 

XQ  0 

EN 

3.   Data  File  Surface-Patch  Box  Structure 
(9  x  9  =  81  patches  on  top  surface). 


FILE:  SPBOX11   DATA 


Ml 


CM  6  CM  MONOPOLE  AT  CENTER  OF  SURFACE  PATCH  FREQ .  1.00  GHZ 
g  KEcULATION  70F7AVERAGESGAIN  AND  INPUT  IMPEDANCE 


G\ 

GUI, 5, 0,0, 1,0, 0,1. 6, 0.016 

GS  0,0,0.1 

GE  1 

EX  0,1,1,0,1,0 

RP  0,31,4,1002,0,0,3,15 

XQ  0 

EN 


4.  Data  File  6cm  Monopole  at  Center  of 
SPBOX.   Input  Admittance  and  Average 
Power  Gain  Calculations 
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FILE:  SPEDGE1   DATA      Ml 

CM  SPEDGE1  11X11  PATCHES  ON  TOP  SURFACE  (6CM  AT  EDGE  3.63CM  FROM  CENTER) 

CE  FREQUENCY  1.00  GHZ 

SM3, 3, 0.5, -0.5, 0,0. 5, -0.5,1 

SCO, 0,0, -0.5,1 

SM3, 3, 0.5, -0.5, 0,0. 5, 0,0 

SCO, 0,0. 5, 0,1 

GX  0,110 

SMI 1,1 1,-0. 5, -0.5, 1,0. 5, -0.5,1 

SCO, 0,C.  5, 0.5,1 

GS  0,0,0.1 

GE  1 

FR  0,0,0,0,1000 

GN  1 

WG 

XQ  0 

EN 

5.   Data  File  Surface-Patch  Box  Structure 
(11  x  11  =  121  patches  on  top  surface). 


FILE:  SPEDGE11  DATA      Ml 

CM  6  CM  MONOPOLE  AT  EDGE  3.63CM  FROM  CENTER  FREQUENCY   1.00   GHZ 

CM  TOP  SURFACE  11X11  PATCHES 

CE  CALCULATION  OF  AVERAGE  GAIN  AND  INPUT  IMPEDANCE 

GF 

GUI, 5, -0.36 36 36, 0,1, -0.36 36 36, 0,1. 6, 0.016 

GS  0,0,0.1 

GE  1 

EX  0,1,1,0,1,0 

RP  0,31,13,1002,0,0,3,15 

XQ  0 

EN 

6.   Data  File  6cm  Monopole  at  Edge  (3.63cm  from 

Center)  of  SPBOX.   Input  Admittance  and  Average 
Power  Gain  Calculations. 
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FILE:  SPCNTR1   DATA     Ml 

CM  SPCNTR   11X11  PATCHES  ON  TOP  SURFACE  (6CM  MONOPOLE  AT  CORNER  5.14  CM 

CM  ON  THE  DIAGONAL  FROM  CENTER 

CE  FREQUENCY  1.00   GHZ 

SM3, 3, 0.5, -0.5, 0,0. 5, -0.5,1 

SCO, 0,0, -0.5,1 

SM3, 3, 0.5, -0.5, 0,0. 5, 0,0 

SCO, 0,0. 5, 0,1 

GX  0,110 

SMI 1,1 1,-0. 5, -0.5, 1,0. 5, -0.5,1 

SCO, 0,0. 5, 0.5,1 

GS  0,0,0.1 

GE  1 

FR  0,0,0,0,1000 

GN  1 

WG 

XQ  0 

EN 

7.   Data  File  Surface-Patch  Box  Structure  for 
Monopole  at  Corner  of  SPBOX. 


FILE:  SPCNTR11  DATA      Ml 

CM  6  CM  MONOPOLE  AT  CORNER  5.14  CM  ON  THE  DIAGONAL  FROM  CENTER 

CM  TOP  SURFACE  11X11  PATCHES 

CM  FREQUENCY   1.0  0   GHZ 

CE  CALCULATION  OF  AVERAGE  GAIN  AND  INPUT  IMPEDANCE 

GF 

GW1, 5, 0.36 36 36, 0.36 3636, 1,0. 36 36 36, 0.36 36 36, 1.6, 0.0 16 

GS  0,0,0.1 

GE  1 

EX  0,1,1,0,1,0 

RP  0,31,13,1002,0,45,3,15 

XQ  0 

EN 

8.   Data  File  6cm  Monopole  at  Corner  (5.14cm  on  the 
Diagonal)  of  SPBOX.   Input  Admittance  and 
Average  Power  Gain  Calculations. 
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PART  B:   NEC  INPUT  DATA  FILES  FOR 
NEAR  ELECTRIC  FIELD  IN  2-D 


FILE:  SNC1      DATA 


CM  SPBOX   7X7  PATCHES  ON  TOP  SURFACE 

CE    FREQUENCY  1.00  GHZ 

SM3, 5, 0.5, -0.5, 0,0. 5, -0.5,1 

SCO, 0,0, -0.5,1 

SM3, 3, 0.5, -0.5, 0,0. 5, 0,0 

SCO, 0,0. 5, 0,1 

GX  0,110 

SM7, 7, -0.5, -0.5, 1,0. 5, -0.5,1 

SCO, 0,0.  5, 0.5,1 

GS  0,0,0.1 

GE  1 

FR  0,0,0,0,1000 

G.'l  1 

WG 

NX 

CM  6  CM  MONOPOLE  AT  CENTER  OF  PATCH  BOX 

CE  NEAR  E-FIELD  ON  X-AXIS  AT  0.0166,0,0.115 

GF 

GW1, 5, 0,0, 1,0, 0,1. 6, 0.016 

G3  0,0,0.1 

GE  1 

EX  0,1,1,0,1,0 

PL  2,2,5,1 

NE  0,200, 1,1, 0.0166, 0,0. 115, 0.01, 0,0 


XQ  0 

EN 


1.   Data  File  Near  Electric  Field  Along 

x-Axis  for  Monopole  at  Center  of  SPBOX 


FILE:  SNC111 


CM  SPBOX   7X7  PATCHES  ON  TOP  SURFACE 

CE  FREQUENCY  1.00  GHZ 

SM3, 3, 0.5, -0.5, 0,0. 5, -0.5,1 

SCO, 0,0, -0.5,1 

SM3, 3, 0.5, -0.5, 0,0. 5, 0,0 

SCO, 0,0. 5, 0,1 

GX  0,110 

SM7, 7, -0.5, -0.5, 1,0. 5, -0.5,1 

SCO, 0,0. 5, 0.5,1 

GS  0,0,0.1 

GE  1 

FR  0,0,0,0,1000 

GN  1 

WG 

NX 

CM  6  CM  MONOPOLE  AT  CENTER  OF  PATCH  BOX 

CE  NEAR  E-FIELD  ON  Z-AXIS  AT  0,0,0.175 

GF 

GW1, 5, 0,0, 1,0, 0,1. 6, 0.016 

GS  0,0,0.1 

GE  1 

EX  0,1,1,0,1,0 

PL  2,2,5,3 

NE  0,1,1,400,0,0,0.175,0,0,0.01 

XQ  0 

EN 

2.   Data  File  Near  Electric  Field  Along 

z-Axis  for  Monopole  at  Center  of  SPBOX 
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FILE:  NEDGEX    DATA      Bl 

CM  SPEDGE1  11X11  PATCHES  ON  TOP  SURFACE  (6CM  AT  EDGE  3.63CM  FROM  CENTER) 

CE  FREQUENCY  1.00  GHZ 

SM3, 3, 0.5, -0.5, 0,0. 5, -0.5,1 

SCO, 0,0, -0.5,1 

SM3, 3, 0.5, -0.5, 0,0.  5, 0,0 

SCO, 0,0. 5, 0,1 

GX  0,110 

SMI 1,1 1,-0. 5, -0.5, 1,0. 5, -0.5,1 

SCO, 0,0. 5, 0.5,1 

GS  0,0,0.1 

GE  1 

FR  0,0,0,0,1000 

GN  1 

WG 

CM  6  CM  MONOPOLE  AT  EDGE  3.63  CM  FROM  CENTER  (11X11  PATCHES  ON  TOP) 

CE  NEAR  E-FIELD  ON  X-AXIS 

GF 

GH  1,5,-0.36  36  36, 0,1,-0.36  3636,0,1.6,0.016 

GS  0,0,0.1 

GE  1 

EX  0,1,1,0,1,0 

PL  2,2,5,1 

HE  0,20  0,1,1,-0.01976  36,0,0.115,0.01,0,0 

XQ  0 

EN 

3.   Data  File  Near  Electric  Field  Along  x-Axis  for 
Monopole  at  Edge  of  SPBOX. 


FILE:  NEDGEY    DATA      Bl 

CM  SPEDGE1  11X11  PATCHES  ON  TOP  SURFACE  (6CM  AT  EDGE  3.63CM  FROM  CENTER) 

CE  FREQUENCY  1.00  GHZ 

SM3, 3, 0.5, -0.5, 0,0. 5, -0.5,1 

SCO, 0,0, -0.5,1 

SM3, 3, 0.5, -0.5, 0,0. 5, 0,0 

SCO, 0,0. 5, 0,1 

GX  0,110 

SMI 1,1 1,-0. 5, -0.5, 1,0. 5, -0.5,1 

SCO, 0,0. 5, 0.5,1 

GS  0,0,0.1 

GE  1 

FR  0,0,0,0,1000 

GN  1 

WG 

NX 

CM  6  CM  MONOPOLE  AT  EDGE  3.63  CM  FROM  CENTER  (11X11  PATCHES  ON  TOP) 

CE  NEAR  E-FIELD  ON  Y-AXIS 

GF 

GW  1,5,-0.36  36  36,0,1,-0.36  3636,0,1.6,0.016 

GS  0,0,0.1 

GE  1 

EX  0,1,1,0,1,0 

PL  2,2,5,2 

NE  6,1,20  0,1,-0.036  36  36,0.0166,0.115,0,0.01,0 

XQ  0 

EN 

4.   Data  File  Near  Electric  Field  Along  y-Axis  for 
Monopole  at  Edge  of  SPBOX. 
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FILE:  NEDGEZ    DATA      Bl 

CM  SPEDGE1  11X11  PATCHES  ON  TOP  SURFACE  (6CM  AT  EDGE  3.63CM  FROM  CENTER) 

CE  FREQUENCY  1.00  GHZ 

SM3, 3, 0.5, -0.5, 0,0. 5, -0.5,1 

SCO, 0,0, -0.5,1 

SM3, 3, 0.5, -0.5, 0,0. 5, 0,0 

SCO, 0,0. 5, 0,1 

GX  0,110 

SMI 1,1 1,-0. 5, -0.5, 1,0. 5, -0.5,1 

SCO, 0,0. 5, 0.5,1 

GS  0,0,0.1 

GE  1 

FR  0,0,0,0,1000 

GN  1 

WG 

NX 

CM  6  CM  MONOPOLE  AT  EDGE  3.63  CM  FROM  CENTER  (11X11  PATCHES  ON  TOP) 

CE  NEAR  E-FIELD  ON  Z-AXIS 

GF 

GW  1,5,-0.36  36  36,0,1,-0.36  36  36,0,1.6,0.016 

GS  0,0,0.1 

GE  1 

EX  0,1,1,0,1,0 

PL  2,2,5,3 

NE  0,1,1,20  0,-0.036  36  36,0,0.175,0,0,0.01 

XQ  0 

EN 

5.   Data  File  Near  Electric  Field  Along  z-Axis  for 
Monopole  at  Edge  of  SPBOX. 


FILE:  NCORNX    DATA      Bl 

CM  SPCNTR   11X11  PATCHES  ON  TOP  SURFACE  (6CM  MONOPOLE  AT  CORNER  5.14  CM 
CM  ON  THE  DIAGONAL  FROM  CENTER 
CE  FREQUENCY  1.00   GHZ 

,0.5,-0.5,1 

,0.5,0,0 
5,1,0.5,-0.5,1 


SM3, 

3,0. 

,5,-0.5,0 

SCO, 

0,0, 

-0.5,1 

SM3, 

3,0. 

5,-0.5,0 

SCO, 

0,0. 

5,0,1 

GX  0 

,  11C 

1 

SM11 

,11, 

-0.5,-0.. 

SCO, 

0,0. 

5,0.5,1 

GS  0 

,0,C 

1.1 

GE  1 

FR  0 

,o,c 

1,0,1000 

GN  1 

HG 

NX 

CM  6 

CM 

MONOPOLE 

CM  MONOPOLE  AT  CORNER  5.14  CM  ON  THE  DIAGONAL  FROM  CENTER 
CM  TOP  SURFACE  11X11  PATCHES 
CM  FREQUENCY   1.00   GHZ 
CE  NEAR  E-FIELD  ON  X-AXIS 
GF 

GW1 ,  5 ,  0  .  36  36  36  ,  0  .  36 36  36  , 1 ,  0  .  36 36  36  ,  0  .  3636 36  , 1  .  6  ,  0  .  016 
GS  0,0,0.1 
GE  1 

EX  0,1,1,0,1,0 
PC  2,2,5,1 

NE  0,20  0,1,1,0.0529636,0.0363636,0.115,0.01,0,0 
XQ  0 
EN 

6.   Data  File  Near  Electric  Field  Along  x-Axis  for 
Monopole  at  Corner  of  SPBOX. 
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FILE:  NCORNY    DATA 


Bl 


CM  SPCNTR   11X11  PATCHES  ON  TOP  SURFACE  C6CM  MONOPOLE  AT  CORNER  5.14  CM 

CM  ON  THE  DIAGONAL  FROM  CENTER 

CE  FREQUENCY  1.00   GHZ 

SM3, 3, 0.5, -0.5, 0,0.  5, -0.5,1 

SCO, 0,0, -0.5,1 

SM3, 3, 0.5, -0.5, 0,0. 5, 0,0 

SCO, 0,0. 5, 0,1 

GX  0,110 

SMI 1,1 1,-0. 5, -0.5, 1,0. 5, -0.5,1 

SCO, 0,0. 5, 0.5,1 

G3  0,0,0.1 

GE  1 

FR  0,0,0,0,1000 

GN  1 

WG 

NX 

CM  6  CM  MONOPOLE  AT  CORNER  5.14  CM  ON  THE  DIAGONAL  FROM  CENTER 

CM  TOP  SURFACE  11X11  PATCHES 

CM  FREQUENCY   1.0  0   GHZ 

CE  NEAR  E-FIELD  ON  Y-AXIS 

GF 

GW1,5,0. 36 36 36, 0.3636 36, 1,0. 36 36 36, 0.36 3636, 1.6, 0.0 16 

GS  0,0,0.1 

1 

0,1,1,0,1,0 

2,2,5,2 

0  , 1 ,  20  0, 1 , 0 . 036  36  36 , 0  .  0529636 ,0.115,0,0.01,0 

0 


GE 
EX 
PL 
NE 
XQ 
EM 


Data  File  Near  Electric  Field  Alorij 
Monopole  at  Corner  of  SPBOX. 


y-Axis  for 


FILE:  NCORNZ    DATA 


Bl 


CM  SPCNTR   11X11  PATCHES  ON  TOP  SURFACE  (6CM  MONOPOLE  AT  CORNER  5.14  CM 
CM  ON  THE  DIAGONAL  FROM  CENTER 
CE  FREQUENCY  1.00   GHZ 
SM3, 3, 0.5, -0.5, 0,0. 5, -0.5,1 
SCO, 0,0,-0  .5,1 
SM3, 3, 0.5, -0.5, 0,0. 5, 0,0 
SCO, 0,0. 5, 0,1 
GX  0,110 

SMI 1,1 1,-0. 5, -0.5, 1,0. 5, -0.5,1 
SCO, 0,0. 5, 0.5,1 
GS  0,0,0.1 
GE  1 

FR  0,0,0,0,1000 
GN  1 
WG 
NX 

CM  6  CM  MONOPOLE  AT  CORNER  5.14  CM  ON  THE  DIAGONAL  FROM  CENTER 
CM  TOP  SURFACE  11X11  PATCHES 
CM  FREQUENCY   1.00   GHZ 
CE  NEAR  E-FIELD  ON  Z-AXIS 
GF 

GW1, 5, 0.36 3636, 0.36 36 36, 1,0. 3636 36, 0.363636, 1.6,  0.016 
GS  0,0,0.1 
_  GE  1 


EX  0,1,1,0,1,0 
PL  2,2,5,3 

NE  0,1,1,20  0,0.0363636,0.036  36  36,0.175,0,0,0 
XQ  0 
EN 
8. 


01 


Data  File  Near  Electric  Field  Along  z-Axis  for 
Monopole  at  Corner  of  SPBOX. 
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APPENDIX  B 
FORTRAN  PROGRAMS  AND  NEC  INPUT i DATA  FILES 

PART  A:   DESCRIPTION  OF  FORTRAN  PROGRAMS  FOR  INTERFACING 
NEC  OUTPUT  WITH  DISSPLA  GRAPHICS  PACKAGE 

PROGRAM  I:   CON1 

This  Fortran  program  is  constructed  in  order  to  generate 
contour  plots  of  magnitude  Epea]<;  (Etotal)  of  the  near 
electric  field  using  NEC  output  in  conjunction  with  the 
DISSPLA  graphics  package  [Ref.  9],  Part  of  NEC's  typical 
output  which  shows  the  calculated  Ep^^,  is  given  in  the 
example  1  data  file,  and  is  produced  using  the  PL  card  with 
options  2,  2,  5,  1  and  an  NE  card  [Ref.  1]  .  In  example  1, 
the  2nd  column  is  the  x-variation,  with  points  spaced  as  is 
defined  by  the  NE  card,  while  the  third  column  is  the 
magnitude  of  Epea]<:  (Etotal)  which  corresponds  to  each  point 
in  space. 

The  CON1  program  reads  the  values  from  the  3rd  column  of 
NEC  near  field  output,  i.e.,  example  1  data  file.  It  then 
arranges  these  values  in  an  NX  by  NY  array,  for  input  to  the 
subroutine  CONTOR,  which  in  turn  creates  the  contour  plots. 
C0N1  also  selects  the  maximum  (Emax)  and  minimum  (Emin) 
values  of  the  electric  field  (3rd  column  of  example  1  data 
file) .  Using  these  values  in  combination  with  the  NLEVEL 
value  (the  number  of  levels  between  the  contours,  which  in 
the  present  study  is  twenty) ,  it  calculates  the  variable  FINC 
which  is  the  contouring  interval,  which  is  needed  by  CONTOR. 
The  values  of  Emin,  Emax,  and  FINC  are  displayed  on  the 

screen.   The  calling  statement  is: 
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CALL  CONTOR  (A,  NX,  NY,  FINC) 

where:  A   -  NX  by  NY  array  containing  the  values  to 
be  contoured 

NX  -  Number  of  columns  in  the  array  (in  the  present 
study  it  is  30) 

NY  -  Number  of  rows  in  the  array  (in  the  present 
study  it  is  30) 

The  values  in  the  array  A  have  been  arranged  by  CON1  as 
follows : 

Z   A(1,NY)       A(2,NY)...  ...A(NX,NY) 

A(1,NY-1)  .  .  . 


A(l,2)  .  .  . 

A(l, 1)        A(2, 1)  .  .  .  .  .  .A(NX,1) 


►  X 

In  the  contour  plots  of  study,  the  magnitude  of  Epg^, 
the  electric  field,  is  presented  in  a  form  of  a  section  of  a 
plane  in  3-dimensional  space  composed  of  a  (30  x  30)  900- 
point-window  in  space.  The  unit  of  magnitude  used  to  express 
the  electric  near  field  in  all  the  plots  is  dB  with  reference 
to  one  volt  per  meter: 


1  =  1 

dB(lV/m)  =  201og10( ) 

lV/m 


When  the  magnitude  of  Ep^^  (E-^otal)  ^s  Plotted  using  NEC 
output  in  the  form  of  example  1,  the  C0N1  FORMAT  is  defined 
as  (17x,  D10.3).  However,  when  Epeak  of  the  inividual 
components  (Ex  or  Ez)  is  plotted,  the  PL  card  with  options, 
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2,2,4,1  is  used  in  conjunction  with  an  NE  card  [Ref.  1].  In 
this  case,  NEC  output  has  the  form  of  example  2.  Columns  in 
the  example  2  data  file  are:  x-variation,  Ex,  x-component 
phase,  Ey,  y-component  phase,  Ez,  z-component  phase.  In  this 
case,  when  the  magnitude  of  individual  components,  x- 
component  (Ex)  or  z-component  (Ez),  is  plotted,  the  FORMAT  in 
C0N1  must  be:  (15x,  D10.3)  and  (55x,  D10.3)  correspondingly. 

CONTOR:  This  subroutine  is  constructed  using  the  DISSPLA 
graphics  package  [Ref.  9]  and  is  modified  for  this  thesis  for 
near  electric  field  contours.  The  subroutine  has  comments  in 
each  step  explaining  its  function.  The  variables  XORIG, 
XSTP,  XMAX,  YORIG,  YSTP,  and  YMAX,  have  to  be  defined 
manually  by  the  user  for  each  case.  For  each  plot,  the 
program  must  be  compiled  by  the  "FORTVS"  command  before 
execution.  In  the  present  study,  the  parameters  labeled 
y-axis  are  actually  the  z-axis  parameters. 

The  subroutine  MYCON  represents  negative  dB  contours  by 
dashed  curves  while  the  positive  ones  are  solid  curves.  The 
number  in  parentheses  of  the  statements,  COMMON  WORK  and  CALL 
BCOMON  [Ref.  9],  must  be  5,000  for  magnitude  plots  while  for 
phase  plots  it  must  be  20,000  because  DISSPLA  needs 
additional  workspace  for  smoothing  the  phase  values. 

PROGRAM  II:   CON3 

This  Fortran  program  generates  contour  plots  of  the 
phase  of  individual  components  (x  or  z  components)  of  the 
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electric  field  from  NEC  output  and  the  DISSPLA  graphics 
package  [Ref.  9].  Part  of  a  typical  NEC  output  is  presented 
in  the  example  2  data  file.  This  output  is  produced  from  the 
the  PL  card  with  options  2,2,4,1  and  an  NE  card  [Ref.  1]  . 
The  columns  presented  in  example  2  are:  x-variation,  Ex,  x- 
component  phase,  Ey,  y-component  phase,  Ez,  z-component 
phase. 

The  program  reads  the  values  from  the  3rd  and  7th 
columns  (Phase  of  x  and  z  components)  and  arranges  them 
in  the  Array  A  consisting  of  Nx  by  Ny  points  using  the 
FORMAT  (25x,  D10.3)  and  the  FORMAT  (65x,  D10.3)  for  the  phase 
of  x  and  z-components  correspondingly.  For  the  case  of  the 
phase  plot  in  the  present  study,  the  array  is  selected  as  100 
x  100  =  10,000  points.  This  is  necessary  because  the  near 
field  region  phase  is  very  complicated,  having  abrupt 
variations  and  the  DISSPLA  program  needs  more  points  for  the 
smoothing  and  accuracy.  The  number  of  levels  between 
contours  in  the  phase  plots  is  selected  to  be:  NLEVEL  =  10 
versus  2  0  (which  was  used  for  the  magnitude  plot)  .  The 
number  in  the  statement:  COMMON  WORK  and  CALL  BCOMON  must  be 
20,000,  because  in  the  case  of  phase  plots,  DISSPLA  needs 
additional  workspace  for  smoothing. 

PROGRAM  III:   CON2 

This  Fortran  program  generates  3-D  plots  of  magnitude 
Epeak  (Etotal)  of  the  near  electric  field.  In  this  type  of 
plot,  a  section  of  a  plane  in  3-dimensions  consisting  of 
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(30  x  30  )  900  points  is  presented.  Electric  field  values 
are  presented  as  a  3rd  dimension  (vertical  axis)  .  The 
contour  plot  appears  at  the  bottom  of  the  graph  (exactly  the 
same  plot  as  the  one  produced  by  CON1)  while  on  the  top,  the 
electric  field  variation  as  a  3-D  surface  is  presented.  The 
program  is  the  same  as  C0N1  except  that  the  user  must  input 
the  number  of  levels  (NLEVEL)  manually  from  the  keyboard  and 
also  has  the  option  to  select  coordinates  of  the  viewpoint: 
VX,  VZ,  and  VE.  CON2  uses  the  subroutine  DMH002  of  DISSPLA 
[Ref.  9]  which  is  listed  at  the  end  of  the  program. 
Parameters  in  the  statements  CALL  GRAF3D,  CALL  GRAF  and  CALL 
RLVEC3  have  to  be  input  by  the  user  and  the  program  has  to  be 
compiled  by  "FORTVS"  before  execution.  All  other 
information,  a  user  may  need,  appears  in  comments  in  the 
program. 

In  order  for  DISSPLA  to  be  automatically  executed  for 
each  of  the  Programs  I,  II,  or  III,  an  EXEC  file  has  been 
constructed  correspondingly  (CON1,  2,  or  3  EXEC).  The  user 
has  to  input  the  filename  of  the  NEC  output  data  file  (which 
has  the  form:  FILENAME  PLOTDATA)  in  the  EXEC  file  and  then 
typing  CON1,  CON2  OR  CON3  in  CMS,  causes  the  execution  of 
DISSPLA  automatically.  For  all  the  programs  the  output  is: 
DISSPLA  METAFILE  T,  which  can  be  copied  to  the  A  or  B-disk 
and  then,  using  the  DISSPOP  option,  the  user  can  select  the 
type  of  output  needed  from  a  menu  (Sherpa,  IBM7  9,  etc.). 
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FILE: 

DIP4      PLOTDATA  Bl 

N8  09 

0 

.242E+01 

0 

.924E-01 

N£09 

0 

.252E+01 

0 

.955E-01 

NS09 

0 

.262E+01 

0 

.945E-01 

N309 

0 

.272E-I-01 

0 

.955E-01 

NS09 

0 

.232E+01 

0 

.960E-01 

N3  0  9 

0 

.292E+01 

0 

.965E-Q1 

NS0  9 

0 

.166E-01 

0 

.511E-02 

N309 

0 

.117E+00 

0 

.624E-02 

N3  09 

0 

.217E+00 

0 

.116E-01 

NS09 

0 

.517E+00 

0 

.168E-01 

N3  0  9 

0 

.417E+00 

0 

.220E-01 

N809 

0 

.517E+00 

0 

.271E-01 

N8  0  9 

0 

.617E+00 

0 

.321E-01 

NS09 

0 

.717E+00 

0 

.369E-01 

MS  09 

0 

.817E+G0 

0 

.416E-01 

N309 

0 

.917E+00 

0 

.461E-01 

N8  09 

0 

.102E+01 

0 

.504E-01 

N8  0  9 

0 

.112E+01 

0 

.545E-01 

N8  09 

0 

.  122E+01 

0 

.534E-0I 

N3  0  9 

0 

.152E+01 

0 

.621E-01 

NS09 

0 

.142E+01 

0 

.655E-01 

NS0  9 

0 

,152E+01 

c 

.688E-01 

N309 

0. 

.162E+01 

0 

.718E-01 

N809 

0. 

172E+01 

0 

746E-01 

MS  09 

0. 

132E+01 

0. 

771E-01 

N8  0  9 

0. 

192E+01 

0. 

795E-01 

NS09 

0. 

202E+C1 

0. 

316E-01 

N809 

0. 

212E+01 

0. 

836E-01 

N309 

0. 

222E+01 

0  . 

853E-01 

N809 

0. 

232E+01 

0. 

869E-01 

N8  0  9 

0. 

242E+01 

0. 

8S3E-01 

N809 

0. 

252E+01 

0. 

895E-01 

NS0  9 

0. 

Z62E+01 

0. 

905E-01 

N3  0  9 

0. 

272E+01 

0. 

914E-01 

NS09 

0. 

2S2E+01 

0. 

922E-01 

N309 

0. 

292E+01 

0. 

92SE-01 

Example    1:      NEC  Output 
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FILE:  DIPH44    PLOTDATA  Bl 

N804  0.106E+01  0 . 501 E-01-0 . 150E+03  0.000E+0C  O.OOOE+OO  0.183E-01  0.460E+02 

NS04  0.109E+01  0.511E-01-0.161L+05  O.OOOE+UO  O.OOOE+00  0.191E-01  0.334E+02 

N804  0.112E+01  0.520E-01-0.174E+03  C.OOOE+00  O.OOOE+OO  0.200E-01  0.207E+02 

NS04  0.115E+C1  0.529E-01  0.174E+03  O.OOOE  +  OO  O.OOOE+OO  0.208E-01  0.769E+01 

N804  C.118E+01  0.538E-01  0.161E+03  O.OOOE+OO  O.OOOE+OO  0 . 216E-01-0 . 553E+01 

NS04  0.121E+01  0.547E-01  0.149E+03  O.OOOE+00  O.OOOE+00  0  .  225E-01-0 . 190E+02 

N504  0.123E+01  0.555E-01  0.135E+03  O.OOOE+OO  O.OOOE+00  0 . 233E-01-0 . 327 E+02 

N804  C.126E+01  0.563E-01  0.122E+03  O.OOOE+OO  O.OOOE+OO  0 .  242E-01-0 .  466E+02 

NS04  0.129E+01  0.571E-01  0.103E+03  O.OOOE+OO  O.OOOE+OO  0 . 251 E-01-0 . 6 OSE+02 

NS04  0.132E+01  0.579E-01  0.942E+02  O.OOOE+OO  O.OOOE+00  0  .  26 OE-01-0 . 751E+02 

N80<  0.135E+01  0.586E-01  0.800E+02  O.OOOE+OO  O.OOCE+00  0 .  268E-01-0 . 898E  +  02 

N804  0.138E+01  0.593E-01  0.655E+02  O.OOOE+OO  O.OOOE+00  0 . 277 E-01-0 . 105E+03 

N804  0.141E+01  0.600E-01  0.508E+02  O.OOOE+OO  O.OOOE+00  0 . 286E-01-0 . 120E+03 

N804  0 . l^^E+01  0.607E-01  0.358E+02  O.OOOE+00  O.OOOE+00  0 . 295E-01-0 . 135E+03 

N804  0.147E+01  0.613E-01  0.205E+02  O.OOOE+OO  O.OOOE+00  0 . 304E-01-0 . 151 E+03 

N804  0.150E+01  0.619E-01  0.505E+01  O.OOOE+OO  O.OOOE+OO  0 .  313E-01-0 .  166E+03 

N804  0.152E+01  0 . 625E-01-0 . 107E+02  O.OOOE+OO  O.OOOE+00  0.321E-01  0.178E+03 

N804  0.155E+01  0.630E-01-0.266E+02  O.OOOE+OO  O.OOOE+00  0.330E-01  0.161E+03 

N804  0.153E+01  0.636E-01-0.428E+02  O.OOOE+OO  O.OOOE+00  0.339E-01  0.145E+03 

N804  0.161E+01  0.641E-01-0.593E+02  O.OOOE+OO  O.OOOE+OO  0.348E-01  0.128E+03 

N804  0.164E+01  0.645E-01-0.759E+02  O.OOOE+OO  O.OOOE+OO  0.357E-01  0.111E+03 

N804  0.167E+01  0  .  650E-01-0 .  923E+02  O.OOOE+OO  O.OOOE+OO  0.365E-01  0.943E+C2 

N804  0.170E+01  0.654E-01-0.110E+03  O.OOOE+OO  O.OOOE+OO  0.374E-01  0.770E+02 

N304  0.173E+01  0.658E-01-0.127E+03  O.OOOE+OO  O.OOOE+OO  0.383E-01  0.595E+02 

N804  0.176E+01  0 . 662E-01-0 . 145E+03  O.OOOE+OO  O.OCOE+OO  0.391E-01  0.417E+02 

N804  0.179E+01  0  .  666 E-01-0 . 163E+03  O.OOOE+00  O.OOOE+OO  0.400E-01  0.238E+02 

N804  0.1S1E+01  0.669E-01  0.179E+03  O.OOOE+OO  O.COOE+00  0.408E-01  0.565E+01 

N804  0.1S4E+01  0.672E-01  0.161E+03  O.OOOE+00  O.OOOE+00  0  .  416E-01-0 . 127E+02 

N804  0.137E+01  0.675E-01  0.143E+03  O.OOOE+OO  O.OOOE+OO  0 . 425E-01-0 . 313E+02 

N804  0.190E+01  0.678E-01  0.124E+03  O.OOOE+00  O.OOOE+00  0 . 433E-01-0 . 500E+02 

N804  0.193E+01  0.6S0E-01  0.105E+03  O.OOOE+00  O.OOOE+00  0 . 441E-01-0 . 690E+C2 

N804  0.196E+01  0.683E-01  0.864E+02  O.COOE+00  O.OOOE+00  0 . 449E-01-0 .881 E+02 

N804  0.199E+01  0.685E-01  0.671E+02  O.OOOE+OO  O.OOOE+OO  0 . 457E-01-0 . 107 E+03 

N804  0.202E+01  0.687E-01  0.477E+02  O.OOOE+OO  O.OCOE+OO  0  .  465E-01-0 . 1 27 E+03 

N804  0.205E+01  0.688E-01  0.281E+02  O.OOOE+OO  O.OOOE+00  0 . 473E-01-0 . 147 E+03 

NS04  0.208E+01  0.690E-01  0.836E+01  O.OOOE+OO  O.OOOE+OO  0  .  481 E-01-0 . 167E+03 

NSC4  0.210E+01  0.691E-01-0.116E+02  O.OOOE+OO  O.OOOE+OO  0.488E-01  0.173E+03 

NS04  0.213E+01  0 . 692E-01-0 . 318E+C2  O.OOOE+OO  O.COOE+00  0.496E-01  0.153E+03 

N804  0.216E+01  0 . 693E-01-0 . 521 E+02  O.OOOE+OO  O.OOOE+OO  0.503E-01  0.133E+03 

N804  0.219E+01  0.694E-01-0.726E+02  O.OOOE+OO  O.OOOE+OO  0.511E-01  0.112E+03 

N804  0.222E+01  0 . 695E-01-0 . 933E+02  O.OOOE+OO  O.OOOE+OO  0.51SE-01  0.912E+02 

N804  0.225E+01  0  .  695E-01-0 . 1 14E+03  O.OOOE+OO  O.OOOE+OO  0.525E-01  0.702E+02 

N804  0.228E+01  0 . 696E-01-0 . 135E+03  O.OOOE+OO  O.OOOE+00  0.532E-01  0.491E+02 

N804  0.231E+01  0 .696E-01-0 . 156E+03  O.OOOE+OO  O.OOOE+00  0.539E-01  0.278E+02 

N804  0.234E+01  0.696E-01-0.173E+03  O.OOOE+OO  O.OOOE+OO  0.545E-01  0.635E+01 

N804  0.237E-r01  0.696E-01  0.161E+03  O.OOOE+OO  O.OOOE+00  0 .  552E-01-0 .  153E+02 

N804  0.239E+01  0.696E-01  0.139E+03  O.OOOE+OO  O.OOOE+OO  0  .  559E-01-0 . 371E+02 

N804  0.242E+01  0.695E-01  0.117E+03  O.OOCE+00  O.OOOE+OO  0 .  565E-01-0 . 590E+02 

NS04  0.245E+01  0.695E-01  0.950E+02  O.OOOE+00  O.OOOE+OO  0 . 571E-01-0 . 811 E+02 

N804  0.243E+01  0.694E-01  0.723E+02  O.OOOE+OO  O.OOOE+00  0 . 573E-01-0 . 103E+03 

N804  0.251E+01  0.693E-01  0.505E+02  O.OOOE+00  O.OOOE+00  0  .  584E-01-0 . 126E+03 

N804  0.254E+01  0.692E-01  0.280E+02  O.OOOE+00  O.OOOE+OO  0  .  590E-01-0 . 148E+03 

N804  0.257E+01  0.691E-01  0.534E+01  O.OOOE+OO  O.OOOE+OO  0 . 595E-01-0 . 171E+03 

N804  0.260E+01  0 . 690E-01-0 . 174E+02  O.OOOE+OO  O.OOOE+OO  0.601E-01  0.166E+03 

N804  0.263E+01  0 .689E-01-0 . 404E+02  O.OOOE+OO  O.OOOE+OO  0.607E-01  0.143E+03 

N804  0.266E+01  0 .688E-01-0 .  635E+02  O.OOOE+OO  O.OOOE+OO  0.612E-01  0.120E+03 

N804  0.268E+01  0 . 686E-01-0 .867E+02  O.OOOE+OO  O.OOOE+OO  0.618E-01  0.968E+02 

N804  0.271E+01  0.685E-01-0.110E+03  O.OOOE+OO  O.OOOE+OO  0.623E-01  0.734E+02 

N804  0.274E+01  0 .683E-01-0 . 134E+03  O.OOOE+OO  O.OOOE+OO  0.628E-01  0.498E+02 

N804  0.277E+01  0 . 681E-01-0 . 157E+03  O.OOOE+OO  O.OOOE+OO  0.633E-01  0.261E+02 

NS04  0.280E+01  0.680E-01  0.179E+03  O.OOOE+OO  O.OOOE+OO  0.638E-01  0.234E+01 

N804  0.283E+01  0.678E-01  0.155E+03  O.OOOE+OO  O.OOOE+OO  0  .  643E-01-0 . 216E+02 

N804  D.286E+01  0.676E-01  0.131E+03  O.OOOE+OO  O.OOOE+00  0 . 647E-01-0 . 457 E+02 

N804  0.289E+01  0.6~74E-01  0.107E+03  O.OOOE+OO  O.OOOE+OO  0 .  652E-01-0 .  699E+02 


2.  Example  2:   NEC  Output 
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FILE:  com 


FORTRAN   Bl 


REAL*8  RO 

DEFINE  THE  ARRAY  ACNX,NY)  WHICH  CONTAINS 

THE  VALUES  TO  BE  CONTOURED . NX=NUMBER  OF  COLUMNS 

AND  NY=NUMBER  OF  ROWS  OF  THE  ARRAY. 

DIMENSION  A(30,30) 
EMIN=1E20 
EMAX=-1E20 
DO  50  1=1,30 

DO  40  J=l,30 
READ  (12,20)RO 
USE  THIS  FORMAT  WHEN  THE  MAGNITUDE  EPEAKC ETOTAL ) 
IS  TO  BE  PLOTTED. NEC  OUTPUT  IN  THIS  CASE  HAS  TO 
BE  PRODUCED  BY  PL  2,2,5,1  AND  NE  CARDS. 


20  FORMAT  (17X,D10.3) 

C  USE  THE  FORMATS  BELOW  WHEN  THE  MAGNITUDE  EPEAK 
C  OF  X  AND  Z  COMPONENTS  CORRESPODINGLY  IS  TO  BE 
C  PLOTTED. NEC  OUTPUT  IN  THESE  CASES  HAS  TO  BE 
C  PRODUCED  BY  PL  2,2,4,1  AND  NE  CARDS. 
C 

C20  FORMAT  (15X.D10.3) 

C20  FORMAT  (55X,D10.3) 

A(J,I)=20XLOG10(RO) 
C  IN  THIS  STEP  THE  PROGRAM  PREDICTS  THE  MINIMUM  AND 
C  AND  MAXIMUM  VALUES  OF  THE  ELECTRIC  FIELD. 
C 

EMIN  =  MIN(EMIN,A(J,D) 
EMAX=MAX(EMAX,A(J,I)) 
40  CONTINUE 

50     CONTINUE 
NX  =  30 
NY  =  30 
C  DEFINE  THE  NUMBER  OF  LEVELS  BETWEEN  THE  CONTOURS. 

NLEVEL=20 
C  THE  PROGRAM  CALCULATES  THE  CONTOURING  INTERVAL. 

FINC  =  (EMAX-EMIfn/(NLEVEL-l) 
C  EMIN,EMAX,FINC  ARE  DISSPLAYED  ON  THE  SCREEN. 

WRITE  (6,*)  EMIN,EMAX,FINC 
C  THE  FOLLOWING  STATEMENT  IS  THE  CALLING  STATEMENT 
C  OF  THE  SUBROUTINE  CONTOR. 

CALL  CONTOR  (A, NX, NY, FINC) 

STOP 

END 


3.       PROGRAM    I:       CON1    FORTRAN.       Contour   plots    of 
Magnitude    for   Near   Electric   Field. 


FILE:    com 


EXEC 


Bl 


FILEDEF    12    DISK    FILENAME    PLOTDATA    A 
EXEC    DISSPLA    CON1 


CON1  EXEC:   Fortran  File  for  Automatic  Execution 
of  DISSPLA. 
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FILE:  CONTOR    FORTRAN   Bl 


SUBROUTINE  CONTOR ( A , NX, NY, FINC) 
THIS  SUBROUTINE  CONTOURS  AN  NX  BY  NY  ARRAY  OF  REGULARLY  SPACED  POINTS 

A  =  AN  NX  BY  NY  ARRAY  OF  REGULARLY  SPACED  POINTS 
NX:  NUMBER  OF  POINTS  IN  THE  X-DIRECTION 
NY:  NUMBER  OF  POINTS  IN  THE  Y-DIRECTION 
FINC:  CONTOUR  INTERVAL 

DIMENSION  A(NX,NY) 
THE  NUMBER  IN  THE  PARENTHESIS  OF  THE  FOLLOWING 
COMMAND  DEFINES  THE  WORKSPACE  IN  THE  DISSPLA 
GRAPHICS  PACKAGE. THIS  NUMBER  HAS  TO  BE  5000  WHEN 
THE  MAGNITUDE  EPEAK  IS  TO  BE  PLOTTED  AND  20000 
WHEN  THE  PHASE  OF  THE  INDIVIDUAL  COMPONENTS  IS 
TO  BE  PLOTTED. IN  THE  CASE  OF  PHASE  PLOT  THE  ARRAY 
IS  LARGER  AND  DISSPLA  NEEDS  ADDITIONAL  WORKSPACE 
FOR  SMOOTHING. SAME  COMMENTS  CORRESPONDS  TO  THE 
"CALL  BCOMON"  STATEMENT. 

COMMON  WORKC5000) 


C  SET  PARAMETERS  FOR  AXES: 
C  DEFINE  X  AXIS  PARAMETERS 
C 

XORIG=0.05<U 

XSTP=0.53 

XMAX=2.9541 
C    DEFINE   Y   AXIS    PARAMETERS    (CORRESPOND   TO    Z   AXIS) 

YORIG=0.115 

YSTP=0.55 

YMAX=5.015 
C  CALL  COMPRS  TO  USE  THE  DISSPOP  OPTION. 

CALL  COMPRS 

CALL  SETCLRC 'CYAN') 


C  SET  PAGE  AND  PLOT  SIZES,  SET  UP  AXES  FOR  PLOT: 

CALL  PAGE(3.5,11.0) 

CALL  BCOMONC5000) 

CALL  AREA2D(6. 0,3.0) 
C 
C  LABEL  AXES: 

CALL  XNAMEC 'X-AXIS  VARIATION  (METERS)* ' , 100 ) 

CALL  YNAMECZ-AXIS  VARIATION  (METERS)*' ,  100) 
C 

CALL  GRAF(XORIG,XSTP,XMAX,YORIG,YSTP,YMAX) 

CALL  FRAME 
C 
C  TITLE  OF  THE  PLOT: 

CALL  HEADINC 'CONTOUR  E-FIELD  (DB  REF  TO  1V/M)$ » , 100 , 2 . , 3) 

CALL  HEADINCMONOPOLE  6CM  AT  CORNER  OF  WIRE  GRID  BOX$', 100,1 

CALL  HEADIN( 'FREQ=lGHZ-(5.3  CM  ON  DIAGONAL )$', 100 , 1 ., 3) 
C 
C  MAKE  CONTOURS  AND  DRAW: 

CALL  SETCLR('RED') 

CALL  CONMIN(3.0) 

CALL  CONANG(60.) 

CALL  CONLIN(0, 'MYCON', » LABELS ',  2, 10) 

CALL  CONMAK(A,NX,NY,FINC) 

CALL  CONTURd,  'LABELS',  'DRAW') 


c 

c 

END 

PLOT: 

CALL  ENDPL(O) 

CALL  DONEPL 

RETURN 

END 

3) 


C**x****xx*x**x*xx*xxxxx*****xx*xxx**x*xxx*xx*x*xxxx**xxx**x*x***x*** 
C 

SUBROUTINE  MYCON( RARAY, IARAY) 

DIMENSION  RARAY(2),IARAY(1) 
C 
C  THIS  ROUTINE  MAKES  NEGATIVE  CONTOURS  DASHED  ,  ZERO  LINES  HEAVIER 


CON00010 
.CON00020 
CON00030 
CON00040 
CON00050 
CON00060 
CON00070 
CON00080 
CON00090 
CON00100 
CON00110 
C0N00120 
CON00130 
CON00140 
CON00150 
CON00160 
CON00170 
CON00180 
CON00190 
CON00200 
CON00210 
CON00220 
CON00230 
CON00240 
CON00250 
CON00260 
CON00270 
CON00220 
CON00290 
CON00300 
CON00310 
CON00320 
CON00350 
CON00340 
CON00350 
CON00360 
CON00370 
CON00330 
CON00390 
CONQ0400 
CON00410 
CONQ0420 
CON00430 
CON00440 
C0N00450 
CON00460 
CON00470 
CON00430 
CONQ0490 
CON00500 
CON00510 
CON00520 
CON00530 
CON00540 
C0N00550 
CON00560 
CON00570 
CON00580 
CON00590 
CON00600 
CON00610 
CON00620 
CON00630 
CON00640 
CON00650 
CON00660 
CON00670 
CON00630 
CON00690 
CON00700 
CON00710 
CON00720 
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FILE:  CONTOR    FORTRAN   Bl 

C  WHILE  THE  POSITIVE  CONTOURS  ARE  PRESENTED  BY  SOLID  LINES.  CON00730 

C  CON00740 

CALL  RESET('DASH')  CON00750 

IF  (RARAYC1)  .GE.  0.)  GO  TO  10  CON00760 

CALL  DASH  CON00770 

10  RARAYC2)  =  1.  CON00780 

IARAY(l)  =  1  CON00790 

IF  (RARAY(l)  .EQ.  0.)  IARAY(l)  =  2  CON00800 

RETURN  CON00810 

END  CON00820 

5.   CONTOR  FORTRAN:   Subroutine  Which  is  Used  to 
Generate  the  Contour  Plots. 
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FILE:  C0N3      FORTRAN   Bl 

C  THIS  PROGRAM  READS  THE  VALUES  OF  THE  PHASE  OF  EX  OR  EZ  COMPONENT  CON00010 

C  FROM  A  NEC  OUTPUT  FILE  (NEC  OUTPUT  IN  THIS  CASE  HAS  TO  CON00020 

C  BE  PRODUCED  BY  A  PL  2,2,4,1  AND  NE  CARDS. THE  OUTPUT  CON00050 

C  IN  THIS  CASE  IS  LIKE  THE  ONE  IN  EXAMPLE  2  DATA  FILE)  CON0004O 

C  AND  USING  THE  SUBROUTINE  CONTCR  FORTRAN  PLOTS  THE  CON00050 

C  VARIATION  OF  THE  PHASE  USING  DISSPLA.  CON0006C 

C  CON00070 

REALX8  RO  CON00080 

C  DEFINE  THE  ARRAY  A(NX,NY)  WHICH  IN  THIS  CASE  IS  COMPOSED  CON00090 

C  OF  100X100=10000  ELEMENTS.  NX=NUMBER  OF  COLUMNS  AND  CON00100 

C  NY=NUMBER  OF  RONS  OF  THE  ARRAY.  CON00110 

DIMENSION  A(100,100)  CON00120 

EMIN=1E20  COH00130 

EMAX=-1E20  CON00140 

DO  50  1=1,100  CON00150 

DO  40  J=l,100  CON00160 

READ  (12,20)RO  CON00170 

C  USE  THIS  FORMAT  WHEN  THE  PHASE  OF  X  COMPONENT  C0N001S0 

C  OF  THE  ELECTRIC  FIELD  IS  TO  BE  PLOTTED  CON0019O 

C                 *  CON00200 

20            FORMAT  (25X,D10.3)  CON00210 

C  USE  THIS  FORMAT  WHEN  THE  PHASE  OF  Z  COMPONENT  CON00220 

C  OF  THE  ELECTRIC  IS  TO  BE  PLOTTED  C0N00230 

C  CON00240 

C20            FORMAT  (65X,D10.3)  CON00250 

ACJ,I)=RO  C0N00260 

C  IN  THIS  STEP  THE  PROGRAM  PREDICTS  THE  MINIMUM  CON00270 

C  AND  MAXIMUM  OF  THE  PHASE  VALUES.  CON00280 

C  CON00290 

EMIN=MIN(EMIN,A(J,I))  CON0  03  0  0 

EMAX=MAX(EMAX,A(J,I))  CON00310 

40          CONTINUE  CON00320 

50     CONTINUE  CON00330 

NX=100  CONO0340 

NY=100  CON00350 

C  DEFINE  THE  NUMBER  OF  LEVELS  BETWEEN  CONTOURS  CON00360 

NLEVEL=10  CON00370 

C  IN  THIS  STEP  THE  PROGRAM  CALCULATES  THE  CONTOURING  CON00330 

C  INTERVAL  CON00390 

C  CON00400 

FINC=(EMAX-EMIN)/(NLEVEL-1)  CON0  0410 

C  THE  VALUES  OF  EMIN, EMAX, FINC  ARE  DISPLAYED  ON  THE  CON00420 

C  SCREEN  CON00430 

C  CON00440 

WRITE  (6,X)  EMIN, EMAX, FINC  CON00450 

C  THE  CALLING  STATMENT  OF  THE  SUBROUTINE  CONTOR  FORTRAN  CON00460 

CALL  CONTOR  (A, NX, NY, FINC)  CON00470 

STOP  CON00480 

END  CON00490 

6.       PROGRAM    II:       CON3    FORTRAN.       Contour  Plots    of   Phase 
for   Near   Electric    Field. 


FILE:    C0N3  EXEC  Bl 


FILEDEF    12    DISK    FILENAME    PLOTDATA    A 
EXEC    DISSPLA    C0N3 


7.   CON3  EXEC:   Fortran  File  for  Automatic  Execution 
of  DISSPLA. 
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FILE:  COM2      FORTRAN   Bl 

C  THIS  PROGRAM  READS  THE  VALUES  OF  THE  MAGNITUDE  EPEAK  OF  E-FIELD  CON00010 

C  FROM  NEC  OUTPUT  FILE  AND  GENERATES  A  PARTICULAR  PLOT  WHICH  HAS  CON00020 

C  THE  3-D  SURFACE  AT  THE  TOP  WHILE  AT  THE  BOTTOM  THE  CONTOUR  CON00030 

C  PLOT  IS  PRESENTED. THE  PROGRAM  USES  THE  SUBROUTINE  "DMH002"  OF  CON00040 

C  DISSPLA  GRAPHICS  PACKAGE  WHICH  IS  ATTACHED  HERE  AT  THE  END.  C0N00050 

C  CON00060 

REALMS  RO  CON00070 

DIMENSION  A(30,30)  CON00080 

CALL  COMPRS  CON00090 

EMIN=1E20  CON00100 

EMAX=-1E20  CON00110 

DO  50  1=1,30  CON00120 

DO  40  J=l,30  CON00130 

READ  (12,20)R0  CON00140 

20            FORMAT  (17X,D10.3)  CON00150 

A(J,I)=20*LOG10(RO)  CON00160 

C  IN  THIS  STEP  THE  PROGRAM  PREDICTS  THE  MAXIMUM  AND  MINIMUM  VALUES  CON00170 

C  OF  THE  MAGNITUDE  OF  THE  ELECTRIC  FIELD.  CON00180 

C  CON00190 

EMIN=MIN(EMIN,A(J,I))  CON0  020  0 

EMAX=MAX(EMAX,A(J,I))  CON0  0210 

40          CONTINUE  CON00220 

50     CONTINUE  CON00230 

NX=30  CON00240 

NY=30  CON00250 

C  DEFINE  THE  NUMBER  OF  LEVELS  (INPUT  FROM  KEYBOARD)  CON00260 

C  CON00270 

WRITE  (6,70)'  ENTER  NLEVEL  >  «  CON00280 

READ  (5,*)  NLEVEL  CON00290 

C      NLEVEL=20  CON00300 

C  IN  THIS  STEP  THE  PROGRAM  CALCULATES  THE  CONTOURING  INTERVAL.  CON00310 

FINC=(EMAX-EMIN)/(NLEVEL-1)  CON0  0320 

C  THE  VALUES  OF  EMIN, EMAX, FINC  ARE  DISSPLAYED  ON  THE  SCREEN.  CON00330 

WRITE  (6,*)  EMIN, EMAX, FINC  CON00340 

C  IN  THIS  STEP  THE  USER  HAS  TO  INPUT  BY  THE  KEYBOARD  THE  COORDINATES     CON00350 

C  OF  THE  VIEW  POINT.  CON00360 

C  CON00370 

WRITE  (6,70)  '  ENTER  VX,VZ,VE,>  »  CON00380 

READ  (5,*)  VX,VZ,VE  CON00390 

C  THIS  IS  THE  CALLING  STATEMENT  OF  THE  SUBROUTINE  "DMH002"  CON00400 

C  CON00410 

CALL  DMH002  (A, NX, NY, FINC, EMIN-30 , EMAX, NLEVEL , VX, VZ, VE)  CON00420 

CALL  DONEPL  CON00430 

70     FORMAT  (A)  CON00440 

STOP  CON00450 

END  CON00460 

C**********************************************************************   CON 00470 

SUBROUTINE  DMHO 02 (ZDAT, NX, NY, FINC, EMIN, EMAX, NLEVEL , VX, VZ, VE)  CON0  048  0 

DIMENSION  ZDAT(30,30)  CON00490 

COMMON  WORK(5000)  CON00500 

CALL  RESET  (3HALL)  CON00510 

CALL  PAGE  (11. ,8.5)  CON00520 

C   SET  AXIS  PARAMETERS  AND  ALPHABETS  CON00530 

CALL  INTAXS  CON00540 

CALL  ZAXANG  (90.)  CON00550 

CALL  SIMPLX  CON00560 

CALL  MX1ALF  (5HL/CST, • * f )  CON00570 

CALL  MX2ALF  (  5HSTAND,  f  V.  •  )  CON00580 

C   DEFINE  3-D  WORK  AREA  AND  AXIS  CON00590 

CALL  AREA2D  (9.5,6.5)  CON00600 

C   DEFINE  THE  LABEL  OF  THE  PLOT  CON00610 

CALL  HEADIN  ('%3-D  ELECTRIC  FIELD  (DB  REF  TO  1V/M)$ ' , 36 , 2 . 0, 3)     CON00620 

CALL  HEADIN('%MONOPOLE  6CM  AT  CORNER  OF  WIRE  GRID  BOX$ ' , 40, 2 . 0,    CON00630 

13)  CON00640 

CALL  HEADIN(»%(5.3  CM  ON  DIAGONAL  FROM  CENTER)$ » , 33, 1 . 0, 3)  CON00650 

CALL  RESET  C6HHEIGHT)  CON00660 

C   DEFINE  THE  LABELS  FOR  THE  THREE  AXES  CON00670 

CALL  X3NAME  (.'XX   VARIATION  (METERS  )$  f ,  21 )  CON00680 

CALL  Y3NAME  ('%Z  VARIATION  (METERS)$ » , 21 )  CON00690 

CALL  Z3NAME  ('^ELECTRIC  FIELDS', 15)  CON00700 

CALL  VOLM3D  (10. ,10. ,8.)  CON00710 

C   DEFINE  3-D  VIEW  POINT  CON00720 
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FILE:  C0N2      FORTRAN   Bl 

CALL  VIEW  (VX,VZ,VE)  CON00730 

C  CALL  GRAF3D(X0RIG,XSTEP,XMAX,Y0RIG,YSTEP,YMAX,EMIN,FINC,EMAX)  CON00740 

C  CON00750 

CALL  GRAF3D  CO . 0541, .58,2. 9541, 0 . 115, .58,3 . 015, EMIN,  FINC,  CON00760 

1EMAX)  CON00770 

C   DEFINE  THE  SURFACE  VIA  A  MATRIX  C0N00730 

CALL  BLSUR  CON00790 

CALL  SURMAT  (ZDAT, 2, 30, 2, 30, 0)  CON00800 

CALL  DASH  CON00S10 

XDEL=0.  CON00820 

DO  231  1=1,3  CON00830 

CALL  BSHIFT  (XDEL,0.)  CON00840 

C   DRAW  BEDPOST  EFFECT  CON00850 

CALL  RESETC 'BASALF')  C0N00860 

CALL  HEIGHTC .2)  CON00870 

CALL  COMPLX  CON00880 

C   DEFINE  A  MESSAGE  ON  THE  PLOT  IF  IT  IS  DESIRED  CON00890 

C     CALL  MESSAGC 'SUBROUTINE  DMH002$ '  ,  100 , 3 . 5,  0 . )  CON00900 

CALL  RESETC 'HEIGHT')  CON00910 

CALL  RESETC 'COMPLX')  CON00920 

C  DEFINE  THE  4  DASHED  LINES  WHICH  CONNECT  THE  4  CORNERS  OF  THE  CONTOUR   CON00930 

C  GRAPH  WITH  THE  3-D  GRAPH. THESE  4  DASHED  LINES  ARE  AN  OPTIONAL  FEATURE . C0N00940 

C  CON00950 

CALL  RLVEC3  ( 0 . 0541 , 0 . 115 , ZDATC 1 , 1 ), 0 . 0541 , 0 . 115 , EMIN, 0 )  CON00960 

CALL  RLVEC3  ( 2 . 9541 , 0 . 115, ZDATC 30 , 1 ) , 2 . 9541 , 0 . 115, EMIN, 0 )  CON00970 

CALL  RLVEC3  C 2 . 9541 , 3 . 015, ZDATC 30 , 30 ) , 2 . 9541 , 3 . 015, EMIN, 0 )  CON00980 

CALL  RLVEC3  C 0 . 0541 , 3 . 015, ZDATC 1 , 30) , 0 . 0541 , 3 . 015, EMIN, 0 )  CON00990 

XDEL=XDEL+.01  CON01000 

231   CONTINUE  CON01010 

CALL  RESET  (3HD0T)  CON01020 

CALL  RESET  (6HBSHIFT)  CON01030 

C   ENTER  GRFITI  LOOP  AND  DEFINE  2-D  PLOT  CON01040 

CALL  GRFITI  ( 0 . , 0  .  , 0  .  0, 1  .  , 0  .  , 0  .  0,  0  .  ,1  .  , 0  .  )  CON01050 

CALL  INTAXS  CON01060 

CALL  AREA2D  (10. ,10.)  CON01070 

C  CALL  GRAF  CXORIG, XSTEP, XMAX, YORIG, YSTEP, YMAX)  CON01080 

C  CON01090 

CALL  GRAF  ( 0 . 0541 ,  . 53 , 2 . 9541 , 0  .  115,  .  58 , 3 . 015)  CON01100 

C   CONTOURING  IS  OPTIONAL  FEATURE  AND  CORRESPONDS  TO  THE  BOTTOM  PLOT     CON01110 

C   SET  CONTUR  PARAMETERS  CON01120 

CALL  BCOMON  (5000)  CON01130 

CALL  CONMAK  (ZDAT, 30 , 30 , FINC)  CON01140 

C     CALL  CONLIN  ( 4, 5HS0L ID, 6HLABELS, 2, 10 )  CON01150 

CALL  CONLIN  ( 0 , 'MYCON ' , » LABELS » , 2, 10)  CON01160 

CALL  CONANG  (90.)  CON01170 

CALL  FRAME  CON01180 

CALL  HEIGHT  (.3)  CON01190 

C     CALL  CONTUR  ( 5, 6HL ABELS, 4HDRAW)  CON01200 

CALL  CONTUR  ( 1 ,' LABELS ',' DRAW' )  CON01210 

C   END  GRFITI  LOOP  CON01220 

CALL  END3GR(0)  CON01230 

CALL  ENDPL  (0)  CON01240 

RETURN  CON01250 

END  CON01260 

C  CON01270 

CXXXXXXXXXXXXXXXXX*XXXXXX**X*XXXXX*X*X**XXXXXXX**XXXXXXXXXXXXXXXXXXXX    CON 01 28  0 

C  CON01290 

SUBROUTINE  MYCON( RARAY, IARAY)  CON01300 

DIMENSION  RARAY(2),IARAY(1)  CON01310 

C  CON01320 

C  THIS  ROUTINE  MAKES  NEGATIVE  CONTOURS  DASHED  AND  THE  ZERO  LINE  HEAVIER .CON01330 

C  THE  POSITIVE  CONTOURS  ARE  PRESENTED  BY  SOLID  LINES  CON01340 

CALL  RESET( 'DASH')  CON01350 

IF  (RARAY(l)- .GE.  0.)  GO  TO  10  CON01360 

CALL  DASH  CON01370 

10  RARAY(2)  =  1.  CON01380 

IARAY(l)  =  1  CON01390 

IF  (RARAY(l)  .EQ.  0.)  IARAY(l)  =  2  CON01400 

RETURN  CON01410 

END  CON01420 

8.   PROGRAM  III:   CON2  FORTRAN.   3-Dimensional  Plots  of 
Magnitude  for  Near  Electric  Field. 
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FILE:  COM2      EXEC     El 


FILEDEF  12  DISK  FILENAME  PLOTDATA  A 
EXEC  DISSPLA  C0.N2 


9.   CON2  EXEC:   Fortran  File  for  Automatic 
Execution  of  DISSPLA. 


PART  B:   NEC  INPUT  DATA  FOR  NEAR  FIELDS 
SURFACE  PATCH  MODELING. 


FILE:  DIP4      DATA      Bl 


CM  DIPOLE  L/2  ON  Z  AXIS . FREQUENCY  1  GHZ 

CE    IN  FREE  SPACE.  NEAR  E-FIELD  (50X30  POINTS) 

GW  1,5, 0,0,-0.  075, 0,0,0. 075,0. 0016 

GE 

FR  0,0,0,0,1000 

EX  0,1,3,0,1,0 

PL  2,2,5,1 

NE  0,30,1,30,0.0166,0,0.115,0.1,0,0.1 

XQ  0 

E.N 


1.   Data  File  Dipole  "\/2    in  Free  Space.   Near 

Electric  Field  (for  Contour  Plot  of  Magnitude) 


FILE:  DIPH44    DATA     Bl 

CM  DIPOLE  L/2  ON  Z  AXIS . FREQUENCY   1  GHZ. IN  FREE  SPACE. 

CE  NEAR  E-FIELD. 100X100  POINTS. 

GW  1,5,0,0,-0.075,0,0,0.075,0.0016 

GE 

FR  0,0,0,0,1000 

EX  0,1,3,0,1,0 

PL  2,2,4,1 

NE  0,10  0,1,10  0,0.0166,0,0.115,0.029,0,0.029 

XQ  0 

EN 

2.   Data  File  Dipole  3/2  in  Free  Space.   Near 
Electric  Field  (for  Contour  Plot  of  Phase) 
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FILE:  M0N111    DATA 


CM  MONOPOLE  6CM  ON  Z  AXIS . FREQUENCY  1  GHZ 

CE  OVER  PERFECT  GROUND.  NEAR  E-FIELD  (30X30  POINTS) 

GN  1,5,0,0,0,0,0,0.06,0.0016 

GE  1 


0.115,0.1,0,0.1 


3.  Data  File  6cm  Monopole  Over  Ground  Plane 
Near  Electric  Field  (for  Contour  Plot  of 
Magnitude) . 


GN 

1 

FR 

0 

.0, 

0, 

.0 

,1000 

EX 

0, 

.1, 

1, 

,0 

,1, 

,0 

PL 

2 

■  Z, 

5, 

1 

HE 

0, 

,3C 

1,] 

.,30, 

,0. 

0166, 

.0 

XQ 

0 

EN 

FILE:  M0NPH1    DATA      Bl 


CM  MONOPOLE  6CM  ON  Z  AXIS . FREQUENCY  1  GHZ 

CE  OVER  PERFECT  GROUND.  NEAR  E-FIELD .( 1 00X100  POINTS) 

GW  1,5,0,0,0,0,0,0.06,0.0  016 

GE  1 

GN  1 

FR  0,0,0,0,1000 

EX  0,1,1,0,1,0 

PL  2,2,4,1 

NE  0,100, 1,100, 0.0166, 0,0. 115, 0.029, 0,0. 029 

XQ  0 

EN 


Data  File  6cm  Monopole  Over  Ground  Plane 
Near  Electric  Field  (for  Contour  Plot  of 
Phase) . 
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FILE:  Nl        DATA      Bl 

CM  SPBOX   7X7  PATCHES  ON  TOP  SURFACE 

CE  FREQUENCY  1.00  GHZ 

SM3, 3, 0.5, -0.5, 0,0. 5, -0.5,1 

SCO, 0,0,-0.5,1 

SM3, 3, 0.5, -0.5, 0,0. 5, 0,0 

SCO, 0,0. 5, 0,1 

GX  0,110 

SM7, 7, -0.5, -0.5, 1,0. 5, -0.5,1 

SCO, 0,0. 5, 0.5,1 

GS  0,0,0.1 

GE  1 

FR  0,0,0,0,1000 

GN  1 

HG 

NX 

CM  6  CM  MONOPOLE  AT  CENTER  OF  PATCH  BOX  .30X30  POINTS 

CE  NEAR  E-FIELD  ON  X-AXIS  AT  0  .  0166 , 0, 0  .  115CC0NT0URED  FIELD) 

GF 

GN1, 5, 0,0, 1,0, 0,1. 6, 0.016 

GS  0,0,0.1 

GE  1 

EX  0,1,1,0,1,0 

PL  2,2,5,1 

ME  0,30,1,30,0.0166,0,0.115,0.1,0,0.1 

XQ  0 

EN 

5.   Data  File  6cm  Monopole  at  Center  of  Surface 
Patch  Box.   Near  Electric  Field  (for  Contour 
Plot  of  Magnitude). 


FILE:  NPH1      DATA      Bl 

CM  SPBOX   7X7  PATCHES  ON  TOP  SURFACE 

CE  FREQUENCY  1.00  GHZ 

SM3, 3, 0.5, -0.5, 0,0. 5, -0.5,1 

SCO, 0,0, -0.5,1 

SM3, 3, 0.5, -0.5, 0,0. 5, 0,0 

SCO, 0,0. 5, 0,1 

GX  0,110 

SM7, 7, -0.5, -0.5, 1,0. 5, -0.5,1 

SCO, 0,0. 5, 0.5,1 

GS  0,0,0.1 

GE  1 

FR  0,0,0,0,1000 

GN  1 

HG 

NX 

CM  6  CM  MONOPOLE  AT  CENTER  OF  PATCH  BOX 

CE  NEAR  E-FIELD  ON  X-AXIS  AT  0  .  0166 , 0 , 0  .  115  .  POINTS  100X100 

GF 

GW1, 5, 0,0, 1,0, 0,1. 6, 0.016 

GS  0,0,0.1 

GE  1 

EX  0,1,1,0,1,0 

PL  2,2,4,1 

NE  0,10  0,1,10  0,0.0166,0,0.115,0.029,0,0.029 

XQ  0 


6.   Data  File  6cm  Monopole  at  Center  of  Surface  Patch 
Box.   Near  Electric  Field  (for  Contour  Plot  of 
Phase).  184 


FILE:  NED1      DATA      Bl 

CM   11X11  PATCHES  ON  TOP  SURFACE  (6CM  AT  EDGE  3.63CM  FROM  CENTER) 

CE  FREQUENCY  1.00  GHZ 

SM3, 3, 0.5, -0.5, 0,0. 5, -0.5,1 

SCO, 0,0, -0.5,1 

SM3, 3, 0.5, -0.5, 0,0. 5, 0,0 

SCO, 0,0. 5, 0,1 

GX  0,110 

SMI  1,1 1,-0. 5, -0.5, 1,0. 5, -0.5,1 

SCO, 0,0. 5, 0.5,1 

GS  0,0,0.1 

GE  1 

FR  0,0,0,0,1000 

GN  1 

WG 

NX 

CM  6  CM  MONOPOLE  AT  EDGE  3.63  CM  FROM  CENTER  (11X11  PATCHES  ON  TOP) 

CE  NEAR  E-FIELD  ON  X-AXIS  (CONTOURS  30X30  POINTS) 

GF 

GN  1,5,-0.36  36  36,0,1,-0.3636  36,0,1.6,0.016 

GS  0,0,0.1 

GE  1 

EX  0,1,1,0,1,0 

PL  2,2,5,1 

NE  0,30,1,30,-0.01976  36,0,0.115,0.1,0,0.1 

XQ  0 

EN 

7.   Data  File  6cm  Monopole  at  Edge  of  Surface  Patch 
Box.   Near  Electric  Field  (for  Contour  Plot  of 
Magnitude) . 


FILE:  NEDPH1    DATA      Bl 

CM   11X11  PATCHES  ON  TOP  SURFACE  (6CM  AT  EDGE  3.63CM  FROM  CENTER) 

CE  FREQUENCY  1.00  GHZ 

SM3, 3, 0.5, -0.5, 0,0. 5, -0.5,1 

SCO, 0,0, -0.5,1 

SM3, 3, 0.5, -0.5, 0,0. 5, 0,0 

SCO, 0,0. 5, 0,1 

GX  0,110 

SMI 1,11, -0.5, -0.5, 1,0. 5, -0.5,1 

SCO, 0,0. 5, 0.5,1 

GS  0,0,0.1 

GE  1 

FR  0,0,0,0,1000 

GN  1 

WG 

NX 

CM  6  CM  MONOPOLE  AT  EDGE  3.63  CM  FROM  CENTER  (11X11  PATCHES  ON  TOP) 

CE  NEAR  E-FIELD  ON  X-AXIS  (PHASE  100  POINTS-ONE  NE  CARD) 

GF 

GW  1,5,-0.36  36  36,0,1,-0.36  36  36,0,1.6,0.016 

GS  0,0,0.1 

GE  1 

EX  0,1,1,0,1,0 

PL  2,2,4,1 

NE  0,10  0,1,10  0,-0.01976  36,0,0.115,0.029,0,0.029 

XQ  0 

EN 

8.   Data  File  6cm  Monopole  at  Edge  of  Surface  Patch 
Box.   Near  Electric  Field  (for  Contour  Plot  of 
Phase).  185 


FILE:  NC01      DATA      Bl 

CM  11X11  PATCHES  ON  TOP  SURFACE  (6CM  MONOPOLE  AT  CORNER  5.14  CM) 

CM  ON  THE  DIAGONAL  FROM  CENTER 

CE  FREQUENCY  1.00   GHZ 

SM3, 5, 0.5, -0.5, 0,0. 5, -0.5,1 

SCO, 0,0,-0 .5, 1 

SM3, 3, 0.5, -0.5, 0,0. 5, 0,0 

SCO, 0,0. 5, 0,1 

GX  0,110 

SMI 1,1 1,-0. 5, -0.5, 1,0. 5, -0.5,1 

SCO, 0,0. 5, 0.5,1 

GS  0,0,0.1 

GE  1 

FR  0,0,0,0,1000 

GN  1 

WG 

NX 

CM  6  CM  MONOPOLE  AT  CORNER  5.14  CM  ON  THE  DIAGONAL  FROM  CENTER 

CM  TOP  SURFACE  11X11  PATCHES 

CM  FREQUENCY   1.00   GHZ 

CE  NEAR  E-FIELD  ON  X-AXIS . (CONTOURS  30X30  POINTS) 

GF 

GW1, 5, 0.36 36 36, 0.36 3636, 1,0.  36 36 36, 0.36 36 36, 1.6, 0.0 16 


GS 

0 

0 

0 

1 

GE 

1 

EX 

0 

1, 

1 

0 

,1 

0 

PL 

2 

2 

5 

1 

NE 

0 

3C 

,1, 

30 

0 

05296 

XQ 

0 

36,0.0363636,0.115,0.1,0,0.1 

0 
EN 

9.  Data  File  6cm  Monopole  at  Corner  of  Surface  Patch 
Box.   Near  Electric  Field  (for  Contour  Plot  of 
Magnitude ) . 

FILE:  NCOPH1    DATA      Bl 

CM  11X11  PATCHES  ON  TOP  SURFACE  (6CM  MONOPOLE  AT  CORNER  5.14  CM) 

CM  ON  THE  DIAGONAL  FROM  CENTER 

CE  FREQUENCY  1.00   GHZ 

SM3, 3, 0.5, -0.5, 0,0. 5, -0.5,1 

SCO, 0,0, -0.5,1 

SM3, 3, 0.5, -0.5, 0,0. 5, 0,0 

SCO, 0,0. 5, 0,1 

GX  0,110 

SMI 1,1 1,-0. 5, -0.5, 1,0. 5, -0.5,1 

SCO, 0,0. 5, 0.5,1 

GS  0,0,0.1 

GE  1 

FR  0, 0,0,0,1000 

GN  1 

WG 

NX 

CM  6  CM  MONOPOLE  AT  CORNER  5.14  CM  ON  THE  DIAGONAL  FROM  CENTER 

CM  TOP  SURFACE  11X11  PATCHES 

CM  FREQUENCY   1.00   GHZ 

CE  NEAR  E-FIELD  ON  X-AXIS . (PHASE  100  POINTS-ONE  NE  CARD) 

GF 

GUI, 5, 0.36 36 36, 0.36 36 56, 1,0.  36 36 36, 0.36 36 36, 1.6, 0.016 

GS  0,0,0.1 

GE  1 

EX  0,1,1,0,1,0 

PL  2,2,4,1 

NE  0,10  0,1,10  0,0.05296  36,0.036  36  36,0.115,0.029,0,0.029 

XQ  0 

EN 

10.  Data  File  6cm  Monopole  at  Corner  of  Surface  Patch 
Box.   Near  Electric  Field  (for  Contour  Plot  of 
Phase).  186 


PART  C:   RADIATION  PATTERNS  FOR  MONOPOLE 
ON  THE  SURFACE  PATCH  BOX. 


FILE:  SPB0XV1   DATA      Bl 

CM  SPBOX   7X7  PATCHES  ON  TOP  SURFACE 

CE  FREQUENCY  1.00  GHZ  (VERTICAL  PATTERN) 

SM3, 5, 0.5, -0.5, 0,0. 5, -0.5,1 

SCO, 0,0, -0.5,1 

SM3, 3, 0.5, -0.5,  0,0. 5,  0,0 

SCO, 0,0.5,0,1 

GX  0,110 

SM7, 7, -0.5, -0.5, 1,0. 5, -0.5,1 

SCO, 0,0. 5, 0.5,1 

GS  0,0,0.1 

GE  1 

FR  0,0,0,0,1000 

GN  1 

NG 

NX 

CM  6  CM  MONOPOLE  AT  CENTER  OF  SURFACE  PATCH  FREQ.  1.00  GHZ 

CM  TOP  SURFACE  7X7  PATCHES 

CE  RADIATION  PATTERNS  (VERTICAL) 

GF 

GUI,  5,  0,0, 1,0,  0,1.  6,  0.0 16 

GS  0,0,0.1 

GE 

EX  0,1,1,0,1,0 

PL  3,1,0,4 

RP  0,131,1,1002,-90,0,1,0 

XQ  0 

EN 


Data  File  6cm  Monopole  at  Center  of  Surface  Patch 
Box.   Vertical  Pattern. 


FILE:  SPB0XH1   DATA      Bl 

CM  SPBOX   7X7  PATCHES  ON  TOP  SURFACE 

CE  FREQUENCY  1.00  GHZ  (HORIZONTAL  PATTERN) 

SM3, 3, 0.5, -0.5, 0,0. 5, -0.5,1 

SCO, 0,0, -0.5,1 

SM3, 3, 0.5, -0.5, 0,0. 5, 0,0 

SCO, 0,0. 5, 0,1 

GX  0,110 

SM7, 7, -0.5, -0.5, 1,0. 5, -0.5,1 

SCO, 0,0. 5, 0.5,1 

GS  0,0,0.1 

GE  1 

FR  0,0,0,0,1000 

GN  1 

WG 

NX 

CM  6  CM  MONOPOLE  AT  CENTER  OF  SURFACE  PATCH  FREQ.  1.00  GHZ 

CM  TOP  SURFACE  7X7  PATCHES 

CE  RADIATION  PATTERNS  (HORIZONTAL) 

GF 

GW1, 5, 0,0, 1,0, 0,1. 6, 0.016 

GS  0,0,0.1 

GE 

EX  0,1,1,0,1,0 

PL  3,2,0,4 

RP  0,1,361,1002,90,0,0,1 

XQ  0 

EN 

2.   Data  File  6cm  Monopole  at  Center  of  Surface  Patch 
Box.   Horizontal  Pattern. 
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FILE:  SPEDGEV1  DATA      Bl 

CM  SPEDGE1  11X11  PATCHES  ON  TOP  SURFACE  (6CM  AT  EDGE  3.63CM  FROM  CENTER) 

CE  FREQUENCY  1.00  GHZ  (VERTICAL  PATTERN-X  AXIS  CUT) 

SMS, 3, 0.5, -0.5, 0,0. 5, -0.5,1 

SCO, 0,0, -0.5, 1 

SM3, 3, 0.5, -0.5, 0,0. 5, 0,0 

SCO, 0,0. 5, 0,1 

GX  0,110 

SMI 1,1 1,-0. 5, -0.5, 1,0. 5, -0.5,1 

SCO, 0,0. 5, 0.5,1 

GS  0,0,0.1 

GE  1 

FR  0,0,0,0,1000 

GN  1 

UG 

NX 

CM  6  CM  MONOPOLE  AT  EDGE  3.63CM  FROM  CENTER  FREQUENCY   1.00   GHZ 

CM  TOP  SURFACE  11X11  PATCHES 

CE  RADIATION  PATTERNS  (VERTICAL-X  AXIS  CUT) 

GF 

,0,1.6,0.016 


GW] 

,5, 

-0 

363636 

0,1 

-0. 

36 

GG 

0 

0 

0 

1 

GE 

EX 

0 

1 

1 

0 

1,0 

PL 

3 

1 

0 

u 

RP 

0 

1S1 

1 

1002 

-90 

0,1 

,0 

XQ 

0 

EN 

3.   Data  File  6cm  Monopole  at  Edge  of  Surface  Patch  Box. 
Vertical  Pattern  (x-Axis  Cut). 


FILE:  SPEDGEV2  DATA      Bl 

CM  SPEDGE1  11X11  PATCHES  ON  TOP  SURFACE  (6CM  AT  EDGE  3.63CM  FROM  CENTER) 

CE  FREQUENCY  1.00  GHZ  (VERTICAL  PATTERN-Y  AXIS  CUT) 

SM3, 3, 0.5, -0.5, 0,0. 5, -0.5,1 

SCO, 0,0, -0.5,1 

SM3, 3, 0.5, -0.5, 0,0. 5, 0,0 

SCO, 0,0. 5, 0,1 

GX  0,110 

SMI 1,1 1,-0. 5, -0.5, 1,0. 5, -0.5,1 

SCO, 0,0. 5, 0.5,1 

GS  0,0,0.1 

GE  1 

FR  0,0,0,0,1000 

GN  1 

WG 

NX 

CM  6  CM  MONOPOLE  AT  EDGE  3.63CM  FROM  CENTER  FREQUENCY   1.00   GHZ 

CM  TOP  SURFACE  11X11  PATCHES 

CE  RADIATION  PATTERNS  (VERTICAL-Y  AXIS  CUT) 

GF 

GW1 ,5,-0.363636,0,1,-0.363636,0,1.6,0.016 

GS  0,0,0.1 

GE 

EX  0,1,1,0,1,0 

PL  3,1,0,4 

RP  "0,131,1,1002,-90,90,1,0 

XQ  0 

EN 

4.   Data  File  6cm  Monopole  at  Edge  of  Surface  Patch  Box. 
Vertical  Pattern  (y-Axis  Cut). 
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FILE:  SPcDGEHl  DATA      Bl 

CM  SPEDGE1  11X11  PATCHES  ON  TOP  SURFACE  C6CM  AT  EDGE  3.63CM  FROM  CENTER) 

CE  FREQUENCY  1.00  GHZ  (HORIZONTAL  PATTERN) 

SMS, 3, 0.5, -0.5, 0,0. 5, -0.5,1 

SCO, 0,0, -0.5,1 

SMS, 3, 0.5,-0.5,0,0.5,0,0 

SCO, 0,0.5,0,1 

GX  0,110 

SMI 1,1 1,-0. 5, -0.5, 1,0. 5, -0.5,1 

SCO, 0,0. 5, 0.5,1 

GS  0,0,0.1 

GE  1 

FR  0,0,0,0,1000 

GN  1 

WG 

NX 

CM  6  CM  MONOPOLE  AT  EDGE  3.63CM  FROM  CENTER  FREQUENCY   1.00   GHZ 

CM  TOP  SURFACE  11X11  PATCHES 

CE  RADIATION  PATTERNS  (HORIZONTAL) 

Gr 

GW1,5, -0. 363636, 0,1,-0. 3636 36, 0,1. 6, 0.016 

GS    0,0,0.1 

Gc 

EX    0,1,1,0,1,0 

PL    3,2,0,4 

RP    0,1,361,1002,90,0,0,1 

XQ    0 

EN 

5.   Data  File  6cm  Monopole  at  Edge  of  Surface  Patch  Box. 
Horizontal  Pattern. 


FILE:  SPCV1     DATA      31 

CM  SPCNTR   11X11  PATCHES  ON  TOP  SURFACE  (6CM  MONOPOLE  AT  CORNER  5.14  CM 

CM  ON  THE  DIAGONAL  FROM  CENTER 

CE    FREQUENCY  1.00   GHZ.  VERTICAL  PATTERN  (X-AXIS  CUT) 

SMS, 3, 0.5, -0.5, 0,0. 5, -0.5,1 

SCO, 0,0, -0.5,1 

SM3, 3, 0.5, -0.5, 0,0. 5, 0,0 

SCO, 0,0. 5, 0,1 

GX  0,110 

SMI  1,1 1,-0.  5, -0.5, 1,0.  5, -0.5,1- 

SCO, 0,0. 5, 0.5,1 

GS  0,0,0.1 

GE  1 

FR  0,0,0,0,1000 

GN  1 

HG 

NX 

CM  6  CM  MONOPOLE  AT  CORNER  5.1<+  CM  ON  THE  DIAGONAL  FROM  CENTER 

CM  TOP  SURFACE  11X11  PATCHES 

CM  FREQUENCY   1.00  GHZ 

CE  RADIATION  PATTERNS  (VERTICAL-  X  AXIS  CUT) 

GF 

GUI,  5, 0.36 36 36, 0.36 36 36, 1,0.  36 36 36, 0.36 36 36, 1.6, 0.016 

GS  0,0,0.1 

GE  1 

EX  0,1,1,0,1,0 

PL  3,1,0,4 

RP  4,181,1,1002,-90,0,1,0 

XQ  0 

EN 


6.   Data  File  6cm  Monopole  at  Corner  of  Surface  Patch  Box 
Vertical  Pattern  (x-Axis  Cut). 
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FILE:  SPCV2 


DATA 


CM  SPCNTR   11X11  PATCHES  ON  TOP  SURFACE  C6CM  MONOPOLE  AT  CORNER  5.14  CM 

CM  ON  THE  DIAGONAL  FROM  CENTER) 

CE  FREQUENCY  1.00   GHZ.  VERTICAL  PATTERN  (45  DEGREES  CUT) 

SM3, 3, 0.5, -0.5, 0,0. 5, -0.5,1 

SCO, 0,0, -0.5,1 

SMS, 3, 0.5, -0.5, 0,0. 5, 0,0 

SCO, 0,0. 5, 0,1 

GX  0,110 

SMI 1,1 1,-0. 5, -0.5, 1,0. 5, -0.5,1 

SCO, 0,0. 5, 0.5,1 

GS  0,0,0.1 

GE  1 

FR  0,0,0,0,1000 

GN  1 

HG 

NX 

CM  6  CM  MONOPOLE  AT  CORNER  5.14  CM  ON  THE  DIAGONAL  FROM  CENTER 

CM  TOP  SURFACE  11X11  PATCHES 

CM  FREQUENCY   1.00   GHZ 

CE  RADIATION  PATTERNS.  VERTICAL  (45  DEGREES  CUT) 

GF 

GN1,  5, 0.36 36 36,  0.36 36 36, 1,0. 36 36 36, 0.36 36 36, 1.6, 0.0 16 

GS  0,0,0.1 

GE  1 

EX  0,1,1,0,1,0 

PL  3,1,0,4 

RP  0,131,1,1002,-90,45,1,0 

XQ  0 


7.   Data  File  6cm  Monopole  at  Corner  of  Surface  Patch 
Box.   Vertical  Pattern  (45°  Cut). 


FILE:  SPCH1 


DATA 


CM  SPCNTR   11X11  PATCHES  ON  TOP  SURFACE  (6CM  MONOPOLE  AT  CORNER  5.14  CM 
CM  ON  THE  DIAGONAL  FROM  CENTER 


GHZ. 
5,-0.5,1 


HORIZONTAL  PATTERN 


0.5,-0.5,1 


CE  FREQUENCY  1.00 

SM3, 3, 0.5, -0.5, 0,0 

SCO, 0,0, -0.5,1 

SMS, 3, 0.5, -0.5, 0,0. 5, 0,0 

SCO, 0,0. 5, 0,1 

GX  0,110 

SM11, 11, -0.5, -0.5,1 

SCO, 0,0. 5, 0.5,1 

GS  0,0,0.1 

GE  1 

FR  0,0,0,0,1000 

GN  1 

WG 

NX 

CM  6  CM  MONOPOLE  AT  CORNER  5.14  CM  ON  THE  DIAGONAL  FROM  CENTER 

CM  TOP  SURFACE  11X11  PATCHES 

CM  FREQUENCY   1.00   GHZ 

CE  RADIATION  PATTERNS  (HORIZONTAL) 

GF 

GW1  ,"5,  0.56  36  36,  0.36  36  36, 1,0.  36  36  36,  0.36  36  36, 1.6,  0.0 16 

GS  0,0,0.1 

GE  1 

EX  0,1,1,0,1,0 

PL  3,2,0,4 

RP  0,1,361,1002,90,0,0,1 

XQ  0 

EN 


Data  File  6cm  Monopole  at  Corner  of  Surface  Patch 
Box.   Horizontal  Pattern. 
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APPENDIX  C 
NEC  INPUT  DATA  FILES 


PART  A:   NEC  INPUT  DATA  FILES  FOR  WIRE  GRID 
MODELING.   STRUCTURE  GEOMETRY 


FILE:  WG11      DATA      Bl 

CM  WIRE  GRID  10  CM.  LONG  WITH  6CM  M0N0P0LE  AT  CENTER 

CM  CELLS  .125  X  .125  GROUND 

CE  DATA  FILE  FOR  PLOTTING  USING  PLOTDGLP  EXEC 

GW  100,4, .5, -.5, 1,0, -.5,1, .01 

GM  1,7, 0,0, 0,0,0, -.125, 100. 100 

GW  20  0,3,  .5, -.5,1, .5, -.5,0,  .01 

GM  2,3,0,0,0,0,-125,0,20  0.201 

GM  3, 3, 0,0,0, -.125, 0,0, 200. 201 

GW  300,4, .5, -.5,1, .5,0,1, .01 

GM  4, 7, 0,0, 0,0,0, -.125, 300. 303 

GW  400, 4,. 5, -.375, 1,0, -.375,1,. 01 

GM  5,2,0,0,0,0,-125,0,400.404 

GW  500, 4,  .37 5, -.5,1,  .37 5,  0,1, .01 

GM  6, 2, 0,0,0, -.125, 0,0, 500. 505 

GX  0,110 

GW  99,3,  .5,0,0,  .5,0,1,  .01 

GM  7,3,0,0,90,0,0,0,099.099 

GW  2, 3,0, -.5, 1,0,.  5,1,.  01 

GW  3, S,.5, 0,1, -.5, 0,1, . 01 

GW  1,5, 0.005, 0,1, 0.005, 0,1. 6,. 016 

GM  0,35,0,0,10,0,0,0,001.001 

GE 

EN 

1.   Data  File  Geometry  of  6cm  Monopole  at  Center 
of  Wire  Grid  Box. 


FILE:  WEDG11    DATA      Bl 

CM  WIRE  GRID  10  CM.  LONG  WITH  6CM  MONOPOLE  AT  EDGEC3.75CM) 

CM  CELLS  .125  X  .125  GROUND 

CE  FILE  FOR  PLOTTING  USING  PLOTDGLP  EXEC 

GW  100,4,  .5, -.5, 1,0, -.5,1,  .01 

GM  1,7, 0,0, 0,0,0, -.125, 100. 100 

GW  200,3,  .5, -.5,1,  .5, -.5,0,  .01 

GM  2,3,0,0,0,0,.125,0,200.201 

GM  3, 3, 0,0,0, -.125, 0,0, 200. 201 

GW  300,4,  .5, -.5,1,  .5,0,1, .01 

GM  4, 7, 0,0, 0,0,0, -.125, 300. 303 

GW  40  0,4, .5, -.375, 1,0, -.375,1, .01 

GM  5,2,0,0,0,0,  .125,0,400.404 

GW  50  0,4, .375, -.5,1,  .375,0,1,  .01 

GM  6, 2, 0,0,0, -.125, 0,0, 500. 505 

GX  0,110 

GW  99,3, .5,0,0,  .5,0,1,  .01 

GM  7,3,0,0,90,0,0,0,099.099 

GW  2, 8,0, -.5, 1,0,  .5,1, .01 

GW  3, 3, .5, 0,1, -.5, 0,1, .01 

GW  1,5, 0.005, 0,1, 0.0 05, 0,1. 6,. 016 

GM  0,35,0,0,10,0,0,0,001.001 

GM  0,0, 0,0,0, -.375, 0,0, 001. 001 

GE 

EN 

2.   Data  File  Geometry  of  6cm  Monopole  at  Edge 
of  Wire  Grid  Box. 
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FILE:  WC0G11    DATA      Bl 

CM  WIRE  GRID  10  CM.  LONG  WITH  MONOPOLE  6  CM  AT  C0RNERC5.3CM  ON  DIAGONAL) 

CM  CELLS  .125  X  .125  GROUND 

CE  FILE  FOR  PLOTTING  USING  PLOTDGLP  EXEC 

GW  100,4,  .5, -.5, 1,0, -.5,1,  .01 

GM  1,7, 0,0, 0,0,0, -.125, 100. 100 

GU  200,8, .5, -.5,1,  .5, -.5,0,  .01 

GM  2,3,0,0,0,0, .125,0,200.201 

GM  3, 3, 0,0,0, -.125, 0,0, 200. 201 

GW  300,4,  .5, -.5,1,  .5,0,1,  .01 

GM  4, 7, 0,0, 0,0,0, -.125, 300. 303 

GW  400, 4,. 5, -.375, 1,0, -.37 5,1,. 01 

GM  5,2,0,0,0,0,.125,0,400.404 

GW  500, 4,. 37 5, -.5,1,. 37 5, 0,1,. 01 

GM  6, 2, 0,0,0, -.125, 0,0, 500. 505 

GX  0,110 

GW  99,8, .5,0,0, .5,0,1,  .01 

GM  7,3,0,0,90,0,0,0,099.099 

GW  2, 8,0, -.5, 1,0, .5,1, .01 

GW  3, 8,. 5, 0,1, -.5, 0,1,. 01 

GW  1,5,0.0  05,0,1,0.005,0,1.6, .016 

GM  0,35,0,0,10,0,0,0,001.001 

GM  0,0,0,0,0, .375, .375,0,001.001 

GE 

EN 

3.   Data  File  Geometry  of  6cm  Monopole  at  Corner 
of  Wire  Grid  Box. 


192 


PART  B:   NEC  INPUT  DATA  FILES  FOR  NEAR  ELECTRIC  FIELDS. 
WIRE  GRID  MODELING 

File:   Wl       DATA    Bl 

CM  WIRE  GRID  10  CM.  LONG 

CM  CELLS  .125  X  .125  GRCUND 

CE. 

GH  10  0,4,  .  5, -.5, 1,0, -.5,1,  . 01 

GM  1,7, 0,0, 0,0,0, -.125, 100. 100 

GH  20 0,3,.  5, -.5,1,  .  5, -.5,0,. 01 

GM  2,3,0,0,0,0,  .125,0,200.201 

GM  5, 3, 0,0,0, -.125, 0,0, 200. 201 

GH  300, A,. 5, -.5,1,. 5, 0,1,. 01 

GM  4, 7, 0,0, 0,0,0, -.125, 300. 303 

GH  400, 4,. 5, -.37 5, 1,0, -.37 5,1,. 01 

GM  5,2,0,0,0..  0,  .125,0,400.404 

GH  500, 4,.  37 5, -.5,1,. 375, 0,1,. 01 

GM  6, 2,  0,0,0, -.125, 0,0, 500. 5 05 

GX  0,110 

GS  0,0,  .1 

GE  1 

FR  0,0,0,0,1000 

GU    1 

WG 

NX 

CM  6  CM  M0N0P0LE  AT  CENTER  OF  HIRE  GRID  FREQ=1  GHZ 

CE  NEAR  E-FIELD  ON  X-AXIS   (CONTOURS  30  POINTS) 

GF 

GH  2, 3,0. -.5, 1,0,. 5,1,. 01 

GH  3, 3,. 5, 0,1, -.5, 0,1,. 01 

Gyl    1,5,  0,0, 1,0, 0,1. 6,. 016 

GS  0, 0,  .1 

GE 

EX  0,1,1,0,1,0 

PL  2,2,5,1 

NE  0,30,1,30,0.0166,0,0.115,0.1,0,0.1 

XO  0 

E.N 

1.   Data  File  6cm  Monopole  at  Center  of  Wire 
Grid  Box  (for  Contour  Plot  of  Magnitude) 

FILE:  HPH1      DATA      Bl 

CM  WIRE  GRID  10  CM.  LONG 

CM  CELLS  .125  X  .125  GROUND 

CE 

GH  100, 4,.  5, -.5, 1,0, -.5,1,  . 01 

GM  1,7, 0,0, 0,0,0, -.125, 100. 100 

GH  20 0,3,.  5, -.5,1,.  5, -.5,0,. 01 

GM  2,3,0, 0, 0,0,  .125, 0,200.201 

GM  3, 3, 0,0,0, -.125, 0,0, 200. 201 

GH    300,4,  .5, -.5,1,  .5,0,1,  .01 

GM  4, 7, 0,0, 0,0,0, -.125, 300. 303 

GH  400,4,  .5, -.37 5, 1,0, -.37 5,1,. 01 

GM  5,2, 0, 0,  0, 0,  .125, 0,400.404 

GH  500, 4,. 375, -.5,1,. 375, 0,1,. 01 

GM  6, 2, 0,0,0, -.125, 0,0, 500. 505 

GX  0,110 

GS  0,0,  .1 

GE  1 

FR  0,0,0,0,1000 

GN  1 

WG 

NX 

CM  6  CM  MONOPOLE  AT  CENTER  OF  WIRE  GRID  FREQ=1  GHZ 

CE  NEAR  E-FIELD  ON  X-AXIS   (PHASE  100  POINTS) 

GF 

GH    2, 3,0, -.5, 1,0,  .5,1,  .01 

GH    3, 3,. 5, 0,1, -.5, 0,1,. 01 

GH  1,5, 0,0, 1,0, 0,1. 6,.  016 

GS  0,0,  .1 

GE 

EX  0,1,1,0,1,0 

PL  2,2,4,1 

NE  0,100,1,100,0.0166,0,0.115,0.029,0,0.029 

XQ  0 

EN 

2.   Data  File  6cm  Monopole  at  Center  of  Wire 
Grid  (for  Contour  Plot  of  Phase). 
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FILE:  WED1      DATA      Bl 

CM  WIRE  GRID  10  CM.  LONG 

CM  CELLS  .125  X  .125  GROUND 

CE 

GW  100,4, .5, -.5, 1,0, -.5,1, .01 

GM  1,7, 0,0, 0,0,0, -.125, 100. 100 

GW  200,3, .5, -.5,1, .5, -.5,0, .01 

GM  2,3,0,0,0,0, .125,0,200.201 

GM  3, 3, 0,0,0, -.125, 0,0, 200. 201 

GW  300,4,  .5, -.5,1,  .5,0,1,  .01 

GM  4, 7, 0,0, 0,0,0, -.125, 300. 303 

GW  40 0,4,. 5, -.37 5, 1,0, -.37 5,1,. 01 

GM  5,2,0,0,0,0,  .125,0,400.404 

GW  500, 4,. 37 5, -.5, 1,  .37  5,0,1,. 01 

GM  6, 2, 0,0,0, -.125, 0,0, 500. 505 

GX  0,110 

GS  0,0,  .1 

CE    1 

FR  0,0,0,0,1000 

GN  1 

WG 

NX 

CM  6  CM  MONOPOLE  AT  EDGEC3.75  CM  FROM  CENTER)  OF  WIRE  GRID  FREQ=1  GHZ 

CE  NEAR  E-FIELD  ON  X-AXIS   (CONTOURS  30  POINTS) 

GF 

GW  2, 3,0, -.5, 1,0,  .5,1,  .01 

GW  3, 3,. 5, 0,1, -.5, 0,1,. 01 

GW  1,5, -.37 5, 0,1, -.375, 0,1. 6, .016 

GS  0,0,  .1 

GE 

EX  0,1,1,0,1,0 

PL  2,2,5,1 

NE  0,30,1,30,-0.0209,0,0.115,0.1,0,0.1 

XQ  0 

EN 

3.      Data   File    6cm  Monopole    at    Edge   of   Wire   Grid   Box 
(for   Contour   Plot    of   Magnitude). 

FILE:    WEDPH1        DATA        Bl 

CM   WIRE   GRID    10    CM.     LONG 

CM   CELLS    .125    X    .125    GROUND 

CE 

GW    100,4, .5,-.5,l,0,-.5,l, .01 

GM   1,7, 0,0, 0,0,0, -.125, 100. 100 

GW   200,8, .5, -.5,1, .5, -.5,0, .01 

GM   2,3,0,0,0,0, .125,0,200.201 

GM   3, 3, 0,0,0, -.125, 0,0, 200. 201 

GW    300,4,  .5, -.5,1,  .5,0,1,  .01 

GM   4, 7, 0,0, 0,0,0, -.125, 300. 303 

GW    400, 4,. 5, -.375, 1,0, -.375,1,. 01 

GM   5,2, 0, 0,0,0,  .125,0,400.404 

GW   500,4, .375, -.5,1, .375,0,1. .01 

GM   6, 2, 0,0,0, -.125, 0,0, 500. 505 

GX    0,110 

GS    0,0,  .1 

GE    1 

FR    0,0,0,0,1000 

GN    1 

WG 

NX 

CM   6    CM   MONOPOLE   AT    EDGEC3.75    CM    FROM   CENTER)    OF   WIRE   GRID    FREQ=1    GHZ 

CE   NEAR    E-FIELD    ON    X-AXIS      (PHASE    100    POINTS) 

GF 

GW   2, 8,0, -.5, 1,0, .5,1, .01 

GW   3,8, .5, 0,1, -.5, 0,1, .01 

GW    1,5, -.375, 0,1, -.375, 0,1. 6, .016 


GS    0,0,  .1 

GE 

EX    0,1,1,0, 

,1,0 

PL    2,2,4,1 

NE   0,100,1; 

100, 

-0 

XQ    0 

EN 

0209,0,0.115,0.029,0,0.029 


4.   Data  File  6cm  Monopole  at  Edge  of  Wire  Grid  Box 
(for  Contour  Plot  of  Phase). 
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FILE:  WC01      DATA 


CM  WIRE  GRID  10  CM.  LONG 
CM  CELLS  .125  X  .125  GROUND 
CE 

GW  100,4,  .5, -.5, 1,0, -.5,1,  .01 
GM  1,7,  0,0, 0,0,0, -.125, 10 0.1 00 
GW  200,8,  .5, -.5,1,  .5, -.5,0,  .01 
GM  2,3,0,0,0,0, .125,0,200.201 
GM  3, 3, 0,0,0, -.125, 0,0, 200. 201 
GW  300,  4,  .  5,  -.5,1,.  5,  0,1,.  01 
GM  4, 7, 0,0, 0,0,0, -.125, 300. 303 
GW  4  0  0,4,  .5, -.37 5, 1,0, -.375,1, .01 
GM  5,2,0,0,0,0, .125,0,400.404 

5,1, .375,0,1, .01 
125,0,0,500.505 


GW 

500,4, 

.  .375,- 

GM 

6, 

,2,0. 

,0,0,-. 

GX 

0. 

,110 

G5 

0. 

,0, .: 

L 

GE 

1 

FR 

0, 

,0,0. 

,0,1000 

GN 

1 

WG 

NX 

CM 

6 

CM  MONOPOL 

GW  1. 

,5, .375, 

GS  0, 

,0,  .1 

GE 

EX    0; 

,1,1,0,1 

PL  2, 

,2,5,1 

NE  0, 

,30,1,30 

XQ  0 

EN 

CM  MONOPOLE  AT  C0RNERC5.3CM  ON  DIAGONAL)  OF  WIRE  GRID  FREQ=1  GHZ 
CE  NEAR  E-FIELD  ON  X-AXIS   (CONTOURS  30  POINTS) 
GF 

GW  2, 3,0, -.5, 1,0,. 5,1,. 01 
GW  3, 8,  .5, 0,1, -.5, 0,1,  .01 

375,1,  .375,  .375,1.6, .016 


0.0541,0.0  37  5,0.115,0.1,0,0.1 


5.   Data  File  6cm  Monopole  at  Corner  of  Wire  Grid  Box 
(for  Contour  Plot  of  Magnitude). 

FILE:  WC0PH1    DATA      Bl 

CM  WIRE  GRID  10  CM.  LONG 

CM  CELLS  .125  X  .125  GROUND 

CE 

GW  100,4, .5, -.5, 1,0, -.5,1, .01 

GM  1,7, 0,0, 0,0,0, -.125, 100. 100 

GW  200,8, .5, -.5,1, .5, -.5,0, .01 

GM  2,3,0,0,0,0, .125,0,200.201 

GM  3, 3, 0,0,0, -.125, 0,0, 200. 201 

GW  300,4, .5, -.5,1, .5,0,1, .01 

GM  4, 7, 0,0, 0,0,0, -.125, 300. 303 

GW  400, 4,. 5, -.37 5, 1,0, -.375,1,. 01 

GM  5,2,0,0,0,0, .125,0,400.404 

GW  500, 4,. 375, -.5,1,. 37 5, 0,1,. 01 

GM  6, 2, 0,0,0, -.125, 0,0, 500. 505 

GX  0,110 

GS  0,0,  .1 

GE  1 

FR  0,0,0,0,1000 

GN  1 

WG 

NX 

CM  6  CM  MONOPOLE  AT  C0RNER(5.3CM  ON  DIAGONAL)  OF  WIRE  GRID  FREQ=1  GHZ 

CE  NEAR  E-FIELD  ON  X-AXIS   (PHASE  100  POINTS) 

GF 

GW  2,8,0, -.5, 1,0, .5,1, .01 

GW  3,8, .5, 0,1, -.5, 0,1, .01 

GW  1,5,  .375,  .375,1, .375,  .375,1.6, .016 

GS  0,0, .1 


GE 

EX 

0 

,1, 

.1, 

,0, 

,1,0 

PL 

2, 

■  2, 

.4; 

,1 

NE 

0, 

,100. 

,li 

,100, 

0 

.0541 

XQ 

0 

EN 

0.0375,0.115,0.029,0,0.029 


6.   Data  File  6cm  Monopole  at  Corner  of  Wire  Grid  Box 
(for  Contour  Plot  of  Phase). 
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